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\ 

X 

Stanton  number  measurements  liave  been  made  for  a transpired  turbu- 
lent boundary  layer  ou  a roug'n  surface.  Teats  were  conducted  at  uniform 
blowing  tractions  trom  0 to  .QOd,  with  uniform  surface  temperature  and 
uniform  tree-stream  velocity.  live  x-Reynolds  number  range  of  these  tests 

and  the  roughness  Reynolds  number  range  from  2U 


a»  -«■—  -y 

was  from  10  to  ^10 

-.*»  '■ 

to  200.  The'  data  are  believed  to  be  accurate  to  within  i .0001  Stanton 


number  units  over  most  of  c2:<s  range.  At  each  test  condition,  several 
velocity  profiles  were  taken  to  measure  the  boundary  layer  growth.  The 
boundary  layer  momentum  thickness  variation  along  the  test  surface  has 
been  used  to  estimate  rough-plate  skin  friction. 

The  data  indicate  th*»  expected  increase  in  both  skin  friction  and 
heat  transfer  due  to  rougliness.  The  data  display  some  unusual  features 
when  plotted  against  boundary  layer  size^Vor  any  given  value  of  F, 
for  all  free-stream  velocities,  the  Stantlm.  number  data  fall  on  a single 
curve  when  plotted  against  enthalpy  thi~Juutsdu  and  the  skin  friction  data 
collapse  (but  less  well)  to  a single  curve  when  plotted  against  momentum 
thickness.  This  behavior  might  have  been  expected  for  a "fully  rough" 
plate  but  was  not  expected  at  the  lower  end  of  the  roughness  Reynolds 
number  range  of  these  tests.  The  blowing  data  indicate  that  within  the 
uncertainty  of  the  experiment,  a form  of  the  Couette  flow  model  used  for 
smooth  surface  boundary  layers  can  also  bt  used  for  rough  surfaces  to 
predict  the  effects  of  blowing  on  Stanton  number. 

Predictions  of  the  experimental  boundary  layers  have  been  carried 
out  using  a finite-difference  boundary  layer  prediction  method  which 
employs  the  Patanker-Spuldlng  finite-difference  formulation  and  a mean 
field  closure  with  a mixing-length  model  employing  van  Driest  type  damp- 
ing. The  effects  of  roughness  liave  been  incorporated  into  this  predic- 
tion program  by  modification  of  tlie  mixing- length  model  used  for  smooth- 
surface  turbulent  boundary  layers. 

The  apparatus  constructed  for  these  tests  is  a closed-loop  wind 
tunnel  using  air  at  essentially  ambient  condition  as  both  the  transpired 
and  free-stream  fluid.  The  rough  surface  consists  of  24  porous  plates 


l v 


forming  an  eight  foot  lcng  test  section.  The  individual  plates  were 
fabricated  by  brazing  together  50-mil  OFHC  copper  balls  stacked  in 
their  most  dense  array.  In  the  test  section  assembly  each  plate  is 
individually  controlled  with  its  own  electric  heater , transpiration  air 
supply  and  instrumentation • The  construction  of  the  test  section  is 
such  that  bloviug,  no  blowing,  and  suction  cases  can  be  tested. 
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NOMENCLATURE 


van  Driest  damping  function  for  a smooth  surface. 

van  Driest  damping  function  for  a rough  surface. 

Blowing  parameter,  F/St. 

Defined  by  Equation  (4-9). 

Specific  heat  of  fluid,  Btu/lb  *F. 

m 

Skin  friction  coefficient,  T0/(pu^/2) 

Pipe  diameter,  ft. 

Parameter  defined  by  the  equation  y+  » ^ in  (Ey+ ) 

Blowing  fraction,  PovQ/pasuw. 

Proportionality  factor  in  Newton's  Second  Law,  lb  ft/lbc  sec^. 

ty  lb  * 

Free-stream  mass  velocity  * p u . lb  / sec  ft  . 

«D  HO  Q 

Roughness  particle  size,  ft. 

Equivalent  sand  grain  roughness  size,  ft. 
mixing  length , ft . 

Mixing  length  at  wall,  ft. 

2 

Mass  flux  through  the  plate  surface,  lb  /sec  ft  . 

m 

Exponent  in  equation  (1-1). 

Ball  radius,  ft. 

Temperature,  *F. 

Boundary  layer  free-stream  velocity,  ft/sec. 

Shear  velocity  ■ ^O8c/P0  . ft/sac. 

Velocity  normal  to  the  test  surface,  ft/sec. 

Average  velocity  in  a pipe  flow,  ft/sec. 

Distance  measured  downstream  xrom  test  section  inlet,  ft. 

Distance  measured  downstream  from  test  section  inlet  to  virtual 
origin,  ft. 

Distance  normal  to  the  surface,  ft. 
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Greek 

2 

a Thermal  diffusivity,  ft  / sec. 

2 

v Kinematic  viscosity,  ft  /sec. 

T Wall  shear  stress,  lb./ft2. 
o r 

6 Boundary  layer  displacement  thickness*  ft. 

99Z  thickness  of  the  momentum  boundary  layer,  ft. 

\ Proportionality  constant  used  in  outer  region  mixing-length, 
formulation  (1  - Ai^). 

tc  Mixing-length  constant,  1 • icy. 

0 Boundary  layer  momentum  thickness,  ft. 

A Boundary  layer  enthalpy  thickness,  ft. 

•X 

p Density,  lb  /ft  . 

y Viscosity,  lb  /sec  ft. 

m 

to  Humidity,  lb  water  vapor/lb.  dry  air 
m m 

Dimensionless  Groups 

Re^  Pipe  Reynolds  number  - dv/v. 

Re^  Roughness  size  Reynolds  number  • ku^/v. 

Re , Roughness  Reynolds  number  " k^u^/v. 

Re^  x- Reynolds  number  - xu^/v. 

Re.  Momentum  thickness  Reynolds  number  “ 6un/v. 

Re^  Enthalpy  thickness  Reynolds  number  “ Au^/v. 

St  Stanton  number  « h/Gc 

ASt  Stanton  number  error,  see  Equation  (2-2). 

StQ  Stanton  number  without  blowing  ■ h/Gc. 

Pr  Prandtl  number  - v/a. 

Prt  Turbulent  Prandtl  number. 
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Subscripts 

® Refers  to  free-etteam  condition*. 

w Refers  to  value  on  fluid  aide  of  wall  - fluid  interface. 

T Refers  to  condition  of  tranapiration  flow  baaeath  poroua  plate. 
H Thermal  (uaed  aa  mining- length  subscript), 
m Momentum  (uwd  aa  mixing- length  subscript), 

a Smooth  surface  value. 

Suparacrlpta 

+ Refers  to  variables  non-diatnsicnalized  in  veil  coordinates; 
y - yuT/v. 
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CHAPTER  1 


INTRODUCTION  AND  BACKGROUND 

Heat  transfer  between  a surface  and  a fluid  stream  is  affected  by 
the  surface  condition,  the  fluid  properties,  the  velocity  field,  and  the 
thermal  boundary  conditions.  Studies  of  this  complex  interaction  are 
generally  conducted  by  fixing  all  but  one  of  these  descriptors  and  vary- 
ing the  remaining  one  to  determine  its  effect.  Host  prior  studies  of 
boundary  layer  heat  transfer  have  been  concerned  with  the  effects  of 
fluid  properties,  temperature  boundary  conditions,  or  main  stream  veloc- 
ity distribution.  Current  interest  in  protecting  surfaces  by  transpira- 
tion or  ablation  1ms  introduced  a class  of  problems  in  which  surface 

1 

roughness  is  an  unavoidable  feature. 

A.  General  Background 

Much  of  what  we  ’ now  about  rough  surface  hydrodynamics  is  based  on 

2 

ideas  and  results  of  the  pipe  flow  experiments  by  Nikuradse  [9].  To 
avoid  the  difficulties  of  defining  an  arbitrary  rough  surface  with  an 
uncertain  distribution  of  roughness  elements,  he  investigated  the  flow 
through  sand-grain  roughened  tubes.  Roughness  elements  for  these  exper- 
iments consisted  of  selected  sand,  carefully  sieved  and  attached  in 
maximum  density  to  the  tube  walls.  The  roughness  In  these  experiments 

was  described  by  a single  parameter,  k , the  size  of  the  sand-grain 

s 

elements.  This  ’sand-grain'  measure  of  roughness  has  become  a standard 
in  skiu  friction  studies.  It  is  still  common  practice  to  express  the 
effect  of  an  arbitrary  roughness  in  terms  of  an  'equivalent  sand-grain 
roughness' , k^. 

This  introductory  paragraph  was  first  used  by  Professor  Moffat  in  the 
original  proposal  to  the  Navy  for  the  rough  surface  research  contract 
that  supported  these  tests.  This  same  introduction  was  repeated  again 
in  each  of  the  12  quarterly  reports  describing  our  progress  on  the  project. 
It  seems  appropriate  that  it  be  used  again  to  Introduce  the  final  report- 
ing of  this  phase  of  the  roughness  work. 

2 

Numbers  in  brackets  refer  to  references  listed  at  the  end  of  this  report. 

1 


The  pipe  flow  experiments  shoved  three  domains  of  behavior  in  terms 
of  the  roughness  Reynolds  number,  Re.;  - uTk  /v.  For  values  of  rJUghness 
Reynolds  number  less  than  5 the  flow  b shaved  as  though  the  surface  were 
smooth.  For  values  of  roughness  Reynolds  ntxnber  greater  than  70,  the 
pipe  friction  factor  became  independent  of  Reynolds  number.  Re  , » vd/v. 
This  state  was  described  as  'fully  rough'.  Values  of  roughness  Reynolds 
number  between  5 and  70  defined  a region  of  'transitonal  roughness'. 

Prandtl  [10]  and  von  Kansan  [11]  used  Nikuradse's  results  (in  1934) 
to  predict  the  behavior  of  the  rough  surface  boundary  layer.  These  pre- 
dictions indicated  that  thi  boundary  layer  would  attain  a 'fully  rough* 
state  such  that  the  skin  friction  would  be  a function  only  of  x/ks« 

Early  rough  surface  boundary  layer  experimental  studies  were  carried  out 
by  Moore  [12]  and  [13].  Moore  looked  at  air  flowing  over  a flat 

plate  roughened  by  regular  arrays  of  square  bars  attached  to  its  surface. 
Ham-*' a experiment  was  carried  out  with  air  flowing  over  surfaces  rough- 
ened by  screens  attached  to  them.  Several  other  experimenters  looked  et 
the  rough  surface  boundary  layer,  but  most  notable  among  these  were  two 
studies  by  Perry  et  al.  [14,15]  and  Liu  et  al.  [16],  The  experiments  by 
Perry  and  his  co-workers  were  carried  out  with  lateral  rectangular  bars 
on  a flat  plate  and  with  zero  and  adverse  pressure  gradients.  They  found 
that  pressure  gradient  did  not  alter  the  roughness  effects  on  the  bound- 
ary layer.  They  also  attempted  to  simulate  a change  in  effective  sand- 
grain  roughness  by  changing  the  spacing  between  the  array  of  bars  used  to 
roughen  the  plate.  They  found  that  changing  the  bar  spacing  did  not 
change  the  effective  sand-grain  roughness  of  the  surface  as  had  been  ex- 
pected. This  led  them  to  define  a second  type  of  roughness  whose  beha- 
vior, unlike  sand-grain  roughness,  could  not  be  correlated  in  terms  of 
an  external  length  scale,  but  scaled  on  the  'logarithmic  asymptote*. 

This  is  the  distance  between  the  top  of  the  roughness  elements  and  the 
apparent  surface  of  the  plate,  the  plane  where  the  velocity  profiles 
extrapolate  to  zero. 

The  experiment  by  Liu  was  conducted  in  water  using  an  array  of  rec- 
tangular bars  to  roughen  the  otherwise  smooth  surface.  Liu  also  varied 
his  bar  spacing  to  change  the  roughness  of  his  surface  and  determined 
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equivalent  sand-grain  roughnesses  for  all  surface  roughness  that  he 
tested.  He  was  not  able,  however,  to  correlate  his  skin  friction  data 
from  differen:  surtacc  roughnesses  with  •»  single  x-Reynolds  number  ex- 
pression. Liu’s  report  also  contains  a very  complete  table  of  both  the 
experimental  and  analytical  work  tliat  had  been  done  till  then  (1966)  for 
both  tough  pipe  and  boundary  layer  type  flowu. 

More  recent  rougli  surface  hydrodynamic  studies  iiave  been  made  by 
Crass  [17],  who  used  a tiydrogen  bubble  technique  to  measure  instantaneous 
velocity  distributions  in  a water  tunnel  above  a sand-grain  roughness. 

His  study  documented  details  of  the  turbulent  structure  near  the  wall,. 

Uu  [18]  used  a floating  element  balance  to  measure  skin  friction  in  an 
air  tunnel  with  a sand-grain  roughened  surface.  These  tests  showed 
'fully  rough'  behavior  and  agreed  well  with  the  Prandtl-Schlichting  pre- 
diction. Tsuji  and  lida  [19]  examined  velocity  profiles  over  rough 
surfaces  and  showed  that  they  could  be  predicted  using  a modified  mixing- 
length  approach,  maintaining  a non-zero  value  of  the  mixing  length  at 
the  wail.  In  a related  study,  Antonia  and  Luxton  [20]  investigated  the 
effect  of  au  abrupt  change  from  a smooth  to  a rough  surface.  In  this 
study  the  rough  surface  was  constructed  from  parallel  square  bars  like 
many  earlier  studies,  but  oue  of  the  bars  was  instrumented  with  static 
pressure  taps  to  provide  a direct  measure  of  form  drag.  The  method  of 
measuring  form  drag  on  individual  roughness  elements  was  also  used  by 
Perry  [14],  Townes  ct  al.  [21,22]  have  studied  the  structure  of  turbu- 
lent flow  in  sand-grain  roughened  pipes. 

B.  Previous  Heat-Tranofer  Experimental  Work 

Much  less  has  been  done  in  the  field  of  heat  transfer.  One  of  the 
first  systematic  rough  surface  experimental  studies  was  carried  out  by 
Nunner  [23].  These  experiments  used  air  flowing  through  rough  pipes. 

The  results  were  used  to  establish  a simple  enpirical  relationship  between 
the  increase  in  Nusselt  number  due  to  roughness  and  the  increase  in  the 
skin  friction.  Several  important  heat  transfer  studies  followed,  notably 
by  Dipprey  and  Sabersky  [24],  Owen  and  Thompson  [25]  and  Gcwen  and  Smith 
126]. 
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Dipprey  and  Sabersky  looked  at  the  flow  of  four  fluids  of  different 
Prandtl  numbers  through  one  smooth  pipe  and  three  rough  pipes.  These 
studies  demonstrated  that  rough  wall  heat  transfer  varied  with  the 
Prandtl  number,  even  in  the  fully  rough  regime  where  molecular  viscosity 
effects  seem  not  important. 

Owen  and  Thompson  proposed  a flow  model  near  the  rough  surface  which 
explains  why  increasing  surface  roughness  can  increase  the  heat  transfer. 
Based  cn  their  ideas,  they  developed  an  expression  for  a sublayer  Stanton 
number  which  was  able  to  correlate  a large  block  of  the  rough  pipe  heat 
transfer  data. 

Gowen  and  Smith  studied  flow  in  several  rough  pipes  using  fluids  of 
three  different  Prandtl  numbers.  In  one  of  their  pipes  they  were  able  to 
make  temperature  profile  measurements  in  the  flow  over  the  rough  wall. 

Like  earlier  investigators,  they  developed  an  expression  for  rough  wall 
Stanton  uumber  based  on  the  rough  wall  skin  friction,  but  in  this  case 
they  were  able  to  include  information  obtained  fron  their  temperature 
profiles. 

There  were,  of  course,  other  rough  wall  heat  transfer  studies,  mostly 
confined  to  internal  flow  applications.  Two  papers  which  sn'smarlze  this 
work  are  by  Sood  and  Johnson  [29J  and  more  recently  by  Norris  [30].  The 
survey  by  Norris  points  out  that  of  the  many  forms  of  rough  wall  heat 
transfer  correlations  that  have  been  proposed,  none  is  simpler  or  has  had 
a great  deal  more  success  <r.  correlating  the  data  than  the  simple  expres- 
sion by  Nunner, 


where  the  exponent  n can  be  expressed  as  a function  of  Prandtl  number. 
The  general  form  of  this  relationship  is  shown  in  Fig.  1.1,  and  Norris's 
recommended  Prandtl  variation  is  shown  in  Fig.  1.2.  An  approximate  ex- 
pression for  the  exponent  is: 

n - .68(Pr)0,215  . (1-2) 


4 


Aii  interesting  feature  of  all  the  pipe  experiments  is  that  there  is  a 
'ceiling*  which  corresponds  to  a rough-to-smooth  skin  friction  ratio  of 
about  4.0.  This  ceiling  seems  to  vary  with  the  type  of  roughness,  with 
saud -grain  roughness  producing  lower  heat  transfer  than  surfaces  rough- 
ened with  wires  and  square  ribs.  This  range  in  the  heat  transfer  ceiling 
is  shaded  in  Fig.  1.1.  For  increases  In  skin  friction  above  this  ceiling, 
there  is  not  a corresponding  increase  in  the  heat  transfer. 

There  continues  to  be  current  interest  in  rough  wail  heat  transfer. 
Much  of  it  is  specifically  related  to  heat  transfer  to  nose  cones  and 
other  re-entry  vehicles.  Unfortunately,  much  of  the  work  done  Is  not 
in  the  open  literature.  Two  related  studies  by  keshotko  [2?]  and  Boldman 
and  Graham  [28]  examine  heat  transfer  to  a nozzle  with  a rough  surface. 

C.  Hough-Surface  Heat-Transfer  Models 

Several  authors  have  presented  prediction  methods  for  rough  wall 
boundary  layer  heat  transfer.  Integral  boundary  layer  prediction 
schemes  which  include  roughness  effects  have  been  described  by  Dvorak 
[■*2,33]  and  by  Chen  131].  Nestler  [34]  has  proposed  a scheme  using  a 
correlation  by  Owen  and  Thompson  to  relate  the  increase  in  Stanton  num- 
ber due  to  roughness  to  the  increase  in  the  skin  friction  and  other 
boundary  layer  parameters.  Finite-dif iereuce  turbulent  boundary  layer 
prediction  schemes  which  Include  roughness  effects  have  been  described 
by  Lumsdaine  et  al.  135]  and  by  McDonald  and  Fish  [36],  bach  of  these 
studies  has  modified  the  mixing-length  distribution  to  Introduce  rough- 
ness effects.  The  Lumsdaine  paper  discusses  a prediction  program  for 
skin  friction  only.  It  employs  a van  Driest  type  damped  mixing-length 
model  near  the  wall.  To  account  for  roughness,  they  fellow  the  recom- 
mendation of  van  DriesL  [39]  and  add  an  additional  term  to  the  mixing- 
length  damping  expression  which  includes  a roughness  parameter. 

The  McDonald  and  Fish  method  includes  both  skin  friction  and  heat 
transfer  calculations.  Roughness  effects  are  entered  through  the  mixing- 
length  damping  by  defining  an  incremental  damping  due  to  roughness  which 
is  added  to  the  uormal  boundary  layer  damping.  This  damping  effect  is 
the  sole  means  ol  including  roughness  effects  into  both  the  hydrodynamic 
and  heat-transfer  predictions  in  this  method. 
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D.  Objectives  of  the  Present  Research 

The  first  objective  of  this  study  was  to  find  a suitable  rough  sur- 
face to  test.  Ue  sought  a surface  which  was  at  the  same  time  repeatable, 
describable,  and  also  porous. 

Hie  second  objective  was  to  measure  the  heat  transfer  from  the  rough 
surface  to  a turbulent  boundary  layer  over  a wide  range  of  free-stream 
velocities.  As  a secondary  part  of  this  objective,  the  overall  hydrody- 
namic performance  of  the  rough  surface  was  to  be  determined  in  terms  of 
boundary  layer  growth  and  skin  friction. 

The  third  objective  was  to  study  the  effect  of  blowing  on  the  rough 
surface  heat  transfer  to  determine  whether  or  not  transpiration  changed 
the  effective  roughness. 

A final  objective  was  adopted  after  the  apparatus  had  been  comple- 
ted and  testing  had  begun:  to  study  the  effects  of  roughness  and  blowing 
on  heat  transfer  in  the  transition  region.  At  the  lowest  test  section 
velocity,  it  was  discovered  that  a fully  developed,  laminar  boundary 
layer  existed  over  the  first  quarter  of  the  test  section.  It  was  de- 
cided that  instead  of  tripping  the  boundary  layer  at  the  test  section 
iuiet  to  immediately  establish  s turbulent  layer,  that  we  would  examine 
the  effects  of  blowing  and  test  section  velocity  in  the  transition 
region  as  these  data  would  be  of  considerable  interest  and  could  be  ob- 
tained at  minimal  cost. 

£.  General  Features  of  the  Progr.  m 

Spherics}  roughness  elements  were  selected  for  this  study.  The 
test  surface  was  constructed  from  o. 050  inch  diameter  O.F.H.C. 
copper  balls  arranged  such  that  the  surface  formed  a regular  array  of 
hemispherical  caps.  This  surface  was  mere  regular  than  sand  grains, 
being  geometrically  describable,  but  more  like  sand-grain  roughness  than 
the  transverse-bar  roughness  elements  used  in  other  tests. 

High-conductivity  copper  balls  were  used  in  fabrication  of  the  sur- 
face to  ensure  a high  thermal  conductivity  of  the  porous  plates.  With 
a high-conductivity  surface,  heater?  can  be  embedded  into  the  plate  and 
yield  a uniform  surface  temperature,  Independent  of  whether  the  transpi- 
ration is  "on"  or  "off".  This  allows  us  to  use  the  same  surface  to  test 


6 


both  the  unblown  as  well  as  the  blown  boundary  layers.  This  is  a neces- 
sary requirement  if  we  are  to  determine  if  b Lowing  changes  the  effects 
of  surface  roughness. 

The  apparatus  constructed  to  test  this  rough  surface  has  the  capa- 
bility of  operating  over  a wide  range  of  free-streaa  conditions  selected 
to  optimize  the  usefulness  of  the  data.  Fig.  1..)  shows  an  operating  map 
of  the  rough  surface  heat  transfer  apparatus.  The  coordinates  used  for 
this  map  describe  the  performances  in  terms  of  a roughness  particle-size 
Reynolds  number  and  an  x-Reynoida  number.  The  range  of  test  section 
velocities  of  the  rig  was  selected  such  that  at  a given  particle-size 
Reynolds  number,  the  x-Reyrolds  number  could  be  varied  by  a factor  of 

10  . In  these  coordinates,  it's  possible  to  overlap  operating  conditions 
with  the  existing  smooth  surface  boundary  layer  apparatus.  The  rig 
operating  parameters  were  also  restricted  to  ensure  an  essentially 
constant  property  boundary  layer  to  minimize  the  effects  of  variable 
fluid  properties. 

The  following  chapters  describe  the  heat  transfer  apparatus  and  its 
qualification  tests,  the  data  which  have  been  obtained  in  these  tests, 
and,  finally,  the  attempts  to  predict  the  data  with  a finito-dif ference 
boundary  layer  prediction  program. 
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CHAPTER  II 


THE  EXPERIMENTAL  APPARATUS 

The  apparatus  used  in  these  experiments  was  constructed  especially 
for  this  study  and  in  the  discussion  that  follows,  it  will  be  referred 
to  as  the  Roughness  Rig.  The  Roughness  Rig  is  located  in  the  Thermo- 
sciences Laboratory  on  the  second  floor  of  the  Mechanical  Engineering 
Bull  ing  at  Stanford.  Its  basic  design  was  copied  after  an  existing 
heat  transfer  facility  which  the  Heat  and  Mass  Transfer  Group  at  Stan- 
ford has  used  over  the  past  six  years  to  investigate  the  transpired 
turbulent  boundary  layer  on  a smooth  surface.  This  heat  and  mass  trans- 
fer apparatus  was  described  first  by  Moffat  [1].  References [2  - 8] 
describe  modifications  that  have  been  made  since  Moffat's  original  ex- 
periments. This  existing  HMT  Rig  is  still  very  much  in  use  and  continues 
to  play  an  important  role  in  the  Heat  and  Mass  Transfer  Group's  research 
activities. 

A.  Description  of  the  Apparatus 

The  Roughness  Rig  is  a closed  loop  wind  tunnel  using  air  at  essentially 
ambient  conditions.  Its  test  surface  consists  of  a 24-segment  porous 
plate,  18  inches  wide  and  8 feet  long.  Figure  2.1  shows  a flow  schematic 
of  the  four  main  rig  systems:  main  air  supply,  transpiration  air  supply, 

plate  heater  electrical  power  system,  and  the  heat  exchanger  cooling 
water  system.  A photograph  of  the  Roughness  Rig  is  shown  in  Figure  2.2. 

The  following  is  a description  of  the  four  main  rig  systems. 

A.l  The  Main  Air  Supply  System 

The  flow  path  of  the  main  air  system  is  as  follows:  (1)  main  air 

blower,  (2)  overhead  ducting  to  an  oblique  header,  (3)  main-stream  heat 
exchanger,  (4)  screen  box,  (5)  nozzle  to  test  section  inlet,  (6)  8 foot 
long  test  section  and  (7)  a multistage  diffuser  which  returns  the  main- 
stream air  back  to  the  blower. 
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The  Roughness  Rig  main  air  supply  blower  is  a 445-BL  Class  3 Buffalo 
Blower  with  a 20  horsepower  belt  drive.  At  rated  conditions,  this  blower 
delivets  8300  cfm  and  develops  a head  of  12  inches  of  water.  The  blower 
and  drive  are  mounted  on  a 900  lb  seismic  base  to  minimize  vibration. 
Flexible  boots  connect  the  blower  to  the  remainder  of  the  tunnel.  The 
main  air  stream  velocity  in  the  test  section  is  varied  by  changing  the 
pulleys  and  belts  on  the  blower  and  drive.  The  range  of  blower  speeds 
attainable  with  this  drive  arrangement  is  from  the  blower  rated  speed 
of  2400  rpm  down  to  372  rpm. 

From  the  main  blower  discharge,  air  is  delivered  through  a 24" 
diameter  overhead  duct  to  an  oblique  inlet  header  on  the  main-stream 
heat  exchanger.  Ducting  construction  is  from  galvanized  sheet  metal 
with  gasketed  Joints  and  sealed  internally  at  all  seams  with  a silastic 
gasket  sealer  material.  The  oblique  inlet  header  to  the  heat  exchanger 
was  designed  based  on  recommendations  by  Wolf  [40].  The  header  shape 
is  spec ified  for  uniform  flow  distribution  and  minimum  pressure  loss. 

The  header  supplies  air  to  a 5 row,  33"  x 48"  heat  exchanger  used 
for  main  air  stream  temperature  control.  The  heat  exchanger  cooling 
water  is  continuously  pumped  from  a holding  tank  through  both  the  main- 
stream heat  exchanger  and  the  transpiration  air  heat  exchanger.  In 
operation,  the  cooling  water  temperature  is  adjusted  until  the  main- 
stream air  at  the  test  section  inlet  is  at  eiioient  temperature.  This 
avoids  the  formation  of  a thermal  boundary  layer  prior  to  the  test  section 
inlet. 

The  heat  exchanger  is  followed  by  a screen  box  containing  four 
stainless  steel,  II 40  mesh,  0.0065"  dia.  wire  screens.  The  screens  were 
sized  to  minimize  the  mean  field  velocity  disturbances  as  well  as  reduce 
the  main  air  stream  turbulence  level.  Based  on  the  work  of  Schubauer 
et.  aL  [41],  this  screen  pack  should  reduce  mean  velocity  disturbances 
by  a factor  of  600  and  turbulent  fluctuations  by  a factor  of  10.  The 
screen  pack  combined  with  the  large  area  contraction  nozzle  provides  a 
uniform  and  well  damped  inlet  velocity  field  to  the  test  section. 

Details  of  the  inlet  velocity  field  measurements  are  reported  in  Section 
F. 
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The  nozzle  accelerates  the  flow  from  the  screen  box  to  the  test 
section  inlet,  a 19.8  to  1 area  contraction.  The  nozzle  wall  design  is 
based  on  a polynoainal  shape  as  recommended  by  Rouse  and  Hassan  [42] . 

This  conforms  with  experience  in  designing  similar  nozzles  at  Stanford. 

The  wall  shape,  which  is  a two-dimensional  contraction,  has  been  chosen 
such  that  both  the  fli-st  and  second  derivatives  of  flow  area  are  zero  at 
the  nozzle  exit.  There  is  a slight  acceleration  in  the  inlet  flow  to 
avoid  possible  separation  in  this  regicn.  This  was  accomplished  by 
designing  the  nozzle  length  for  40"  and  cutting  off  the  first  inch  when 
the  nozzle  was  actually  built.  A teledeltos  model  of  each  of  the  nozzle 
walls  was  used  to  insure  that  the  nozzle  smoothly  accelerated  the  flow 
and  there  was  not  separation  at  inlet  or  exit. 

The  test  section  consists  of  the  test  plate  assembly,  two  side  walls 
and  a movable  upper  surface.  At  the  inlet,  the  test  section  is  20  inches 
wide  and  4 inches  high.  The  upper  surface  is  pivoted  at  the  inlet  so  that 
it  can  be  adjusted  to  give  either  an  increasing  or  decreasing  flow  area 
in  the  flow  direction.  The  top  and  side  walls  are  constructed  from  contin- 
uous sheets  of  1/2  inch  thick  plexiglass.  The  top  is  sufficiently  flexi- 
ble so  that  it  can  be  warped  to  accomodate  the  nonlinear  growth  of  the  test 
section  boundary  layers.  The  side  walls  contain  two  sets  of  static  press- 
ure taps,  on  2 and  12  inch  centers  in  the  flow  direction.  The  pressure 
taps  on  12  inch  centers  in  the  flow  direction  are  used  to  position  the 
top  to  achieve  the  desired  free  stream  velocity  distribution.  The  tests 
described  here  were  conducted  with  uniform  free  stream  velocity.  To 
obtain  this  test  condition,  the  top  wall  is  set,  experimentally,  for 
each  run  by  adjusting  it  until  there  is  no  measurable  change  in  static 
pressure  along  the  test  section  centerline.  In  practice,  local  deviations 
of  +.001  inches  of  water  were  accepted.  A probe  sled,  which  spans  the  test 
section,  locks  onto  the  side  walls  in  fixed  positions  over  the  center  of 
each  of  the  24  individual  test  plates.  Probes,  supported  from  this  sled 
extend  down  through  access  holes  in  the  movable  top.  In  addition  to  the 
access  holes  along  the  test  section  centerline  the  second,  eleventh  and 
twenty-third  (of  twenty-four)  of  the  individual  test  plates  are  pro- 
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vided  with  a full  set  of  access  holes  extending  across  the  width  of  the 
test  section  on  one  inch  centers.  These  are  primarily  to  examine  boundary 
layer  uniformity  in  the  transverse  direction. 

Flow  from  the  test  section  exits  into  a multistage  vaned  diffuser. 

The  diffuser  inlet  has  a movable  top  to  allow  alignment  with  the  test 
section  top.  The  inlet  section  is  followed  by  three  separate  vaned  two- 
dimensional  diffusion  sections  finally  emptying  into  a plenum  box.  The 
diffuser  area  ratio  between  test  section  and  plenum  box  is  in  excess  of 
7:1.  The  plenum  box  is  connected  to  the  main  blower  inlet  through  a 
flexible  boot.  £ach  of  the  individual  diffuser  sections  were  designed 
based  on  the  experimental  work  of  Cocharan  and  Kline  [43]  and  the  total 
diffuser  recovers  approximately  40%  of  the  kinetic  energy  head  in  the 
flow  leaving  the  test  section.  A small  charging  blower  is  attached  to 
the  plenum  box  to  make  up  air  leakage  from  the  tunnel.  This  blower  is 
used  to  ensure  that  the  test  section  static  pressure  is  equal  to  the 
ambient  pressure  in  the  laboratory.  This  avoids  any  leakage  of  air  into 
or  out  of  the  test  section  which  might  disturb  the  two-dimensionality 
of  the  flow. 

Appendix  A contains  a brief  description  of  the  construction  details 
of  the  screen  box,  nozzle  and  diffuser. 
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A. 2 The  Transpiration  Air  Supply  System 

The  flow  path  of  the  transpiration  air  system  is  as  follows:  (1) 

inlet  air  filter  box,  (2)  transpiration  blower,  (3)  transpiration  heac 
exchanger,  A)  transpiration  flow  header  box,  (5)  individual  delivery 
tubes  to  each  of  the  porous  plate  sections. 

The  purpose  of  the  inlet  filter  box  is  to  prevent  clogging  and  con- 
tamination of  the  porous  plates.  It  is  constructed  using  5 micron  re- 
tention filter  felt  material  and  is  sized  to  have  60  square  feet  of  filter 
area. 

The  filter  box  is  connected  to  a Buffalo  type  V,  size  25  blower. 

This  blower  is  driven  by  a 15  horsepower,  3600  RPM  motor,  connected 
directly  to  the  blower.  In  this  configuration,  the  blower  runs  at  one 
speed  and  transpiration  flow  is  controlled  by  ball  valves  in  the  indivld- 
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ual  delivery  lines.  This  particular  blower  was  designed  for  a higher 
volumetric  flow  than  is  required  at  low  transpiration  rates.  In  order 
to  avoid  off-design  operation,  the  blower  was  equipped  with  dump  valves 
so  that  a part  of  its  discharge  air  could  be  dumped  back  into  the  lab. 

It  was  found,  however,  that  the  blower  performance  was  stable  over  the 
entire  flow  range  of  interest  in  these  experiments  without  venting  its 
discharge . 

Transpiration  air  from  the  blower  is  delivered  through  a 10  inch 
flexible  duct  to  a large  box  containing  the  transpiration  air  heat  ex- 
changer and  bypass  system  which  allows  partial  or  total  heat  exchanger 
by-pass.  A 5 row,  18  x 24  inch  heat  exchanger  is  ut.ed  in  the  transpira- 
tion system.  The  heat  exchanger  receives  its  cooling  water  from  the 
same  recirculating  water  system  which  supplies  the  main  air  heat  ex- 
changer. From  the  heat  exchanger  enclosure,  the  transpiration  air  is 
passed  through  several  mixing  turns,  then  dumped  into  a header  box.  In 
operation,  the  heat  exchanger  and  its  by-pass  system  are  used  to  insure 
that  the  transpiration  air  is  delivered  to  the  header  at  ambient  temper- 
ature. This  avoids  the  formation  of  temperature  gradients  in  the  de- 
livery system  which  acts  Itself  like  a heat  exchanger  between  the  tran- 
spiration air  and  the  laboratory. 

The  transpiration  header  box  is  physically  located  under  the  test 
section.  Individual  supply  lines,  one  for  each  porous  plate,  are  con- 
nected to  the  side  of  the  header  box.  The  first  section  of  each  supply 
line  is  a 3 foot  long  riser  reaching  from  the  header  box  to  a ball  valve 
used  to  control  the  flow  to  each  plate.  In  each  riser  section,  upstream 
of  the  valve,  is  a hot-wire  type  flowmeter.  Details  concerning  the 
fabrication  and  calibration  of  these  flowmeters  are  described  in  Appendix 
B.  The  hot-wire  flowmeter  was  chosen  because  of  its  wide  sensitivity 
range.  The  range  of  flows  which  must  be  accurately  metered  to  cover  the 
operating  range  of  the  Roughness  Rig  is  from  less  than  1 cfm  to  over  50 
cfm.  To  cover  this  range  with  a loss-based  metering  system  such  as  an 
orifice  or  rotameter  would  require  two  or  more  parallel  meters  in  each 
line.  To  complete  the  supply  line  from  the  control  valve  to  the  test 
plate  assembly,  1 inch  flexible  tubing  was  used.  This  final  link  in 
the  supply  line  maken  the  Roughness  Rig  transpiration  system  compatible 
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with  the  existing  HMT  Rig  plate  assemblies.  It  would  be  possible  to 
install  the  existing  smooth  plate  assemblies  on  the  Roughness  Rig  and 
they  would  be  completely  compatible  with  the  transpiration  system.  To 
minimize  the  into  ’action  between  the  transpiration  air  and  the  surround- 
ings, the  header  box  and  supply  lines  have  been  insulated.  The  header 
box  is  covered  with  a layer  of  aluminum  foil  backed,  rock  wool  insula- 
tion and  the  supply  lines  have  been  encased  in  a cardboard  and  masonite 
zone  box  to  insure  a thermally  uniform  environment  for  the  system. 
Thermocouples  in  the  individual  supply  lines  indicate  temperature  dif- 
ferences on  the  order  of  one  degree  Fahrenheit  from  end  to  end  of  the 
8 foot  header  box,  the  hotter  end  being  nearest  the  transpiration  blower. 

When  a blown  boundary  layer  is  eatablished  in  the  test  section,  the 
transpiration  system  operates  in  an  open  loop  mode.  Transpiration  air 
is  drawn  in  through  the  filter  box  and  delivered  to  the  test  section. 

It  is  then  dumped,  downstream  at  the  diffuser  plenum  box.  The  small 
charging  blower  which  is  used  to  control  test  section  pressure  is  dis- 
connected and  air  is  bled  from  the  diffuser  plenum  box  through  a control 
slide  valve.  Test  section  static  pressure  can  easily  be  balanced  to 
match  the  ambient  pressure  by  controlling  the  flow  of  air  through  the 
control  valve. 

When  boundary  layer  suction  is  required  the  flexible  duct  connection 
is  moved  to  the  suction  side  of  the  blower  and  the  blower  discharge  is 
used  to  charge  the  diffuser  plenum  to  provide  test  section  static  pres- 
sure control. 

A. 3 The  Plate  Header  Electrical  Power  System 

The  plate  heater  power  supply  is  a 750  amp,  24  kilowatt  Lincoln  Arc 
Welder.  The  field  r^-istance  of  the  welder  has  been  modified  to  fix  its 
output  at  22  volts  O.C.  Power  connections  are  made  to  the  ground  and 
'high  point'  taps  on  the  welder  to  minimize  line  voltage  droop  with  in- 
creasing load.  Power  is  delivered  to  a bus  bar  box  mounted  on  the  side 
of  the  Roughness  Rig  through  overhead  copper  bus  bar  system.  Each  in- 
dividual plate  has  its  own  heater  which  consists  of  a single  piece  of 
#26  AWG  stranded  copper  wire  with  Irradiated  PVC  insulation.  This  wire 
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is  laced  back  and  forth  in  eight  grooves  in  the  back  of  the  plate.  Ends 
of  the  heater  wire  are  connected  to  il  12  AWG  stranded  copper  leads  which 
pass  through  the  side  cf  the  aluminum  support  casting.  One  heater  lead 
is  connected  to  a precision  ammeter  shunt,  one  for  each  plate,  and  then 
to  the  ground  bus  bar  in  the  box.  The  other  lead  is  connected  to  a 
power  transistor  mounted  on  the  other  bus  bar  in  the  box.  The  power 
transistor  is  part  of  the  power  control  circuit  for  each  plate  heater. 
Plate  power  is  controlled  by  individual  amplifier  circuits,  one  for  each 
plate,  which  can  adjust  heater  voltage.  All  of  the  control  circuit  ele- 
ments except  the  power  transistors  on  the  bus  bar  hav*  been  reduced  to 
a printed  circuit  assembly, 6 channels  per  board  and  are  mounted  in  a 
card  box  behind  a control  panel  in  the  instrument  console  located  on  the 
left  of  the  rig,  sse  Fig.  2.2.  On  this  control  panel  are  potentiometers 
with  digital  read  out  faces  and  small  edge  voltmeters  which  Indicate 
heater  voltage  for  each  channel.  The  plate  voltage  and  therefore  the 
power  is  proportional  to  the  potentiometer  setting  for  that  channel.  By 
recording  pot  settings  from  a particular  run,  test  conditions  can  be  very 
nearly  repeated  by  resetting  the  pots  to  their  former  netting.  A detailed 
description  of  the  power  control  circuits  is  given  in  Appendix  C.  The 
bus  bar  on  which  the  power  transistors  are  mounted  acts  as  a heat  sink 
for  the  power  transistors.  Under  certain  conditions  these  transistors 
reject  a substantial  heat  load  and  it  is  necessary  to  provide  water 
cooling  of  the  bus  bar.  This  was  accomplished  by  soft  soldering  1/2  x 1 
inch  rectangular  copper  wave  guide  into  two  parallel  grooves  milled 
lengthwise  in  the  surface  of  the  8 foot  long,  3 x 3/8  inch  copper  bus  bar. 
The  wave  guide  is  used  as  a cooling  water  passage  to  reduce  transistor 
temperatures.  A magnetic  valve  wired  into  the  starter  switch  circuit 
for  the  welder  power  supply  automatically  turns  the  bus  bar  cooling 
water  on  when  the  welder  is  started.  A time  delay  relay  set  for  a one 
minute  delay  insures  that  the  bus  bar  cooling  continues  until  after  the 
welder  has  coasted  down  once  it  has  been  stopped. 

Heater  power  measurements  are  made  by  measuring  the  voltage  drop 
across  the  heater  and  across  a precision  shunt  in  the  heater  circuit. 
Heater  voltage  and  shunt  voltage  leads  are  carried  in  shielded  pairs  to 
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selector  switch  read-out  stations.  These  selector  switches  are  'zone 
boxed'  to  prevent  any  thermally  induced,  stray  signals  which  could  intro- 
duce errors  in  power  readings. 

These  read-out  stations  are  also  equipped  with  VIDAR  plug  connections 
which  will  allow  automatic  data  acquisition  at  a late-  date.  The  ammeter 
shunts  used  for  power  read-out  were  individually  calibrated,  in  place, 
against  a high  precision  shunt  from  the  Thermosciences  Measurements  Center. 
In  addition  to  individual  shunt  calibration,  the  shunt  connections  were 
potted  with  RTV  compound  to  insulate  them  against  possible  external  thermal 
effects  and  each  of  the  24  shunts  were  provided  with  an  air  cooling  system 
to  insure  that  there  would  be  no  stray  thermally  Induced  EMF  signals. 

A. 3 The  Heat  Exchanger  Cooling  Water  System 

Cooling  water  for  both  the  main  air  and  transpiration  air  heat  ex- 
changers is  supplied  from  a single  loop  which  continuously  circulates 
cooling  water  through  both  heat  exchangers  from  an  insulated  holding  tank. 
By  maintaining  a high  water  circulation  rate,  in  excess  of  20  GPM,  temper- 
ature gradients  across  the  heat  exchangers  are  minimized.  This  ensures 
uniform  temperature  in  the  air  being  cooled.  The  main  air  heat  exchanger 
is  sized  to  remove  the  heat  load  from  the  main  blower  as  well  as  the 
heater  power  added  in  the  test  section.  The  transpiration  air  heat  ex- 
changer removes  the  heat  load  due  to  the  transpiration  blower.  Temper- 
ature control  of  the  cooling  water  is  achieved  by  dumping  a portion  of 
the  heat  exchanger  return  water  and  making  it  back  up  from  the  water 
supply  main.  The  make-up  water  is  mixed  in  the  holding  tank  to  provide 
a buffer  against  temperature  fluctuations  in  the  supply  water.  It  was 
found  that  the  cooling  water  temperature  dominated  the  exit  air  temper- 
ature from  both  the  main  and  transpiration  air  heat  exchangers.  This 
was  not  unexpected  for  an  air-water  system  with  high  effectiveness  heat 
exchangers  such  as  are  used  here.  By  holding  the  circulating  water 
temperature  just  below  the  ambient  temperature,  the  exit  air  temperature 
from  both  heat  exchangers  can  be  held  at  the  laboratory  ambient  temper- 
ature. 
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B.  The  Test  Plate  Assembly 

The  test  plate  assembly  forms  the  bottom  surface  of  the  Roughness 
Rig  test  section.  Its  design  has  been  copied  after  the  test  plate 
assembly  used  In  the  existing  HMT  Rig.  It  is  22  inches  wide,  2 inches 
wider  than  the  inside  dimension  of  the  test  section,  and  96  inches  long 
in  the  flow  direction.  It  is  made  up  of  four  subassemblies,  bolted  to 
a common  support  structure.  Each  subassembly  consists  of  an  aluminum 
casting  in  which  are  mounted  6 individual  porous  plates,  each  4 inches 
wide  and  18  inches  long.  The  porous  region  of  the  test  surface  thus 
is  96  inches  long  and  18  inches  wide  . A photograph  of  a machined  casting 
is  shown  in  Fig.  2.3.  A cross  section  through  one  of  the  plates  and  the 
compartment  in  the  aluminum  casting  beneath  it  is  shown  in  Fig.  2.4. 

Transpiration  air  enters  through  the  air  delivery  tube  from  the 
flow  control  valve.  Like  the  smooth  plate  HMT  Rig,  a one  inch  delivery 
tube  is  used  for  each  compartment.  The  air  jet  entering  the  lower  com- 
partment is  diverted  by  a baffle  plate  to  avoid  direct  impingement  on 
the  pre-plate.  The  air  flow  in  this  first  inlet  plenum  may  be  poorly 
distributed  and  may  have  significant  eddies.  To  protect  the  working 
plate,  the  upper  surface  of  the  inlet  plenum  is  a porous  bronze  pre- 
plate which  provides  a pressure  drop  type  damping  of  the  maldistribution. 
Air  passing  through  the  pre-plate  is  distributed  no  worse  than  the 
variation  in  permeability  of  the  commercially  available  pre-plates. 

Directly  above  the  pre-plate  is  the  thermocouple  location  for 
measuring  the  temperature  of  the  air  delivered  to  the  porous  plates.  A 
layer  of  honeycomb  material  is  fastened  to  the  bottom  surface  of  the 
porous  plates.  This  honeycomb  has  hexagonal  cells,  3/16  inch  in  diameter 
and  is  3/8  inch  thick.  The  honeycomb  acts  as  a 'flow  straightener'  for 
the  flow  from  the  pre-plates.  It  also  eliminates  eddies  or  jets  which 
could  cause  both  flow  and  temperature  maldistributions.  The  honeycomb 
not  only  isolates  the  plate  from  contact  with  eddies  in  the  underbody 
region,  it  significantly  reduces  the  radiation  view  factor  between  the 
back  of  the  plate  and  the  underbody.  Essentially  all  that  the  back  of 
the  plate  can  'see'  through  the  honeycomb  is  the  pre-plate  directly  below 
it. 
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Thermocouples  for  measuring  the  plate  temperatures  are  set  into 
the  porous  plate  with  their  junctions  located  .068  Inches  from  the  ball 
crests  which  form  the  rough  surface.  The  surface  and  air  temperature 
thermocouples  are  0.010  Inch  diameter  iron-constantan  thermocouples. 

The  electric  heater  wires  are  stranded,  26  AUG  copper  wire.  Ir- 
radiated polyvinyl  chloride  insulation  for  the  wire  was  chosen  for  the 
temperature  tolerance  of  this  material  and  a series  of  glue  bond  tests. 
The  heater  wires  are  glued  into  grooves  spaced  0.43  inches  apart  on  the 
back  of  the  porous  plates  with  a thin  bead  of  Armstrong  A-31  epoxy.  Each 
plate  is  supported  along  its  long  edges  by  a 1/32  inch  thick  linen- 
reinforced  phenolic  stand-off.  These  strips  serve  to  thermally  isolate 
the  plates  from  the  castings  and  minimize  heat  losses.  The  phenolic 
strips  are  glued  to  the  sides  of  adjoining  plates  and  inserted  into 
slots  milled  into  the  web  dividing  two  adjacent  compartments  in  the 
casting.  Care  was  taken  to  provide  an  air  tight  seal  entirely  around 
each  compartment  to  prevent  leakage  between  adjacent  compartments. 

Assembly  of  the  plates  into  the  casting  is  a hand  operation  since 
adjacent  plates  must  fit  together  with  a minimum  disturbance  to  the 
surface  roughness  pattern.  Small  irregularities  between  plates  were 
corrected  by  hand-finishing  and  matching  each  group  of  six  plates. 

The  first  step  in  the  assembly  was  to  glue  the  hand  fitted  group  of  six 
plates  together  to  form  a continuous  sheet,  18  x 24  inches.  Between 
each  pair  of  plates,  the  phenolic  stand-offs  are  also  glued  in  place. 

This  gluing  assembly  is  made  by  clamping  the  plates  together  on  a pre- 
cision surface  table,  under  a high  clamping  load.  The  next  step  was 
installation  of  the  heater  wires  in  each  of  the  plates.  The  casting 
and  plates  were  next  assembled  together,  again  on  a surface  table, 
clamped  together.  Figure  2.5  shows  a photograph  of  the  clamping  arrange- 
ment used  for  gluing  the  casting  and  plate  assembly  together.  The  final 
step  was  the  installation  of  the  honeycomb,  thermocouples,  pre-plates 
and  the  bottom  plates  on  the  castings.  On  the  top  surface,  balsa  filler 
strips  were  used  to  fill  the  gap  between  the  plate  edges  and  the  top 
surfaces  of  the  casting  side  rails.  Figure  2.6  is  a photograph  of  the 
back  of  the  plate  assembly  showing  the  honeycomb,  thermocouples  and  pre- 
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plates.  A photograph  of  the  top  surface  of  the  casting  assembly  in 
place  on  the  support  structure  is  shown  in  Fig.  2.7. 

Five  surface  temperature  thermocouples  are  Installed  in  each  plate. 
One  is  in  thr<  geometric  center  of  the  plate.  The  remaining  four  are 
arranged  1 1/3  inch  upstream  and  downstream  and  three  inches  to  the 
right  and  left  of  the  central  thermocouple.  This  array  can  sense  tamper- 
ature  gradients  both  in  the  flow  direction  and  in  the  transverse  direc- 
tion. The  area  enclosed  by  this  thermocouple  pattern  is  considered  to 
be  the  measuring  area,  the  rest  of  the  plate  acts  as  a ‘guard'  to  mini- 
mize thermal  and  hydrodynamic  edge  effects. 

Casting  temperatures  are  measured  by  thermocouples  installed  in 
every  other  web  of  each  casting  at  the  test  section  center  line.  Temper- 
ature control  of  the  casting  ia  provided  by  cooling  wat  • tubes  in  the 
casting  webs.  To  maintain  casting  temperature  as  close  as  possible  to 
the  transpiration  air  temperature,  the  casting  cooling  water  is  first 
circulated  through  a copper  tube  heat  exchanger  Inside  the  transpiration 
header  box.  The  purpose  of  holding  the  casting  temperature  at  the  tran- 
spiration air  temperature  is  to  minimize  heat  transfer  between  the  tran- 
spiration air  and  the  casting.  Any  heat  exchange  after  the  temperature 
of  the  gas  stream  has  been  measured  in  the  casting  underbody  would  in- 
troduce an  energy  balance  error.  Casting  temperature  control  is  also 
useful  to  evaluate  conduction  losses  from  the  plates  as  will  be  discussed 
in  Section  E of  this  chapter. 

C.  Porous  Plates 

When  construction  was  began  on  the  Roughness  Rig,  considerably  effort 
went  into  the  problem  of  plate  fabrication.  The  porous  plates  used  in 
the  smooth  plate  HMT  rig  were  made  by  sintering  carefully  screened  bronze 
particles  whose  average  diameter  was  .005  Inches  in  a stainless  steel 
mold,  under  Stanford  supervision.  This  was  necessary  because  commercially 
sintered  plates  exhibited  permeability  variations  of  20  to  301.  The  de- 
cision to  use  0.0S0  spherical  particles  on  the  Roughness  Rig  virtually 
eliminated  sintering  es  a possibility. 


C.l  Plate  Fabrication 


When  construction  of  the  Roughness  Rig  began,  several  west  coaat 
vendors  were  contacted  and  asked  to  submit  samples  of  sintered  porous 
plates  for  porosity  testing.  Only  one  vendor,  the  Commercial  Filter 
Corporation  of  Los  Angeles,  submitted  a sample.  This  was  constructed 
using  bits  of  a .030  inch  bronze  wire  sintered  into  an  8 inch  square, 

1/2  inch  thick  plate.  Tests  revealed  that  the  porcsity  of  this  sample 
was  highly  non-uniform  and  it  was  not  acceptable.  At  this  same  time, 
the  decision  was  made  to  fabricate  the  plates  by  assembling  particles, 
rather  than  by  a bulk  process,  and  O.F.H.C.  copper  balls  were  chosen  as 
idsal  elements.  The  problem  then  became  one  of  joining  the  balls  to- 
gether. Brazing  seamed  most  likely  to  succeed,  providing  the  braze 
material  could  be  uniformly  applied.  Discussions  with  brazing  specialists 
revealed  that  a nickel-phosphorus  brazing  alloy  could  be  deposited  on  the 
copper  balls.  This  established  the  technique  which  was  used  to  fabricate 
the  porous  pistes  used  for  the  Roughness  Rig.  Tha  plates  are  constructed 
from  O.F.H.C.  copper  balls,  .030  lnchss  in  diameter.  The  braze  material 
was  provided  by  plating  each  ball  with  .0005  inches  of  elactroleas  nickel. 
The  balls  were  then  arranged  by  hand  into  rows  and  layers  Inside  copper 

molds  in  their  most  dense  array  ami  fired  in  an  inert  atmosphere  furnace 

to  Just  above  the  melting  temperature  of  the  nickel  plating.  This  re- 
sulted in  a brazing  together  of  the  ball  pack  into  a uniformly  poroua 

plate.  The  platiug  material,  as  it  melted,  formed  small  fillets  between 

adjacent  balls  at  each  contact  point.  A close-up  photograph  of  the  sur- 
face of  the  porous  plate  is  shown  in  Fig.  2.8.  The  final  dimensions  of 
each  plate  were  18.0  x 4.0  x 0.5  inches.  This  method  of  piste  fabrica- 
tion, although  tedious,  provided  s well  defined  surface  roughness  pattern 
and  uniformly  porous  piste  for  the  transpiration  experiments.  A more  de- 
tailed description  of  the  plate  fabrication  technique  is  given  in  Ap- 
pendix D. 

C.2  Plate  Heaters 

Eight  heater  grooves  were  formed  in  the  back  surface  of  each  plate 
by  omitting  eight  of  the  rows  of  balls.  The  allowable  heater  wire  spacing 
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is  governed  by  the  requirement  that  the  temperature  of  the  top  surface 
of  the  plate  remain  uniform.  The  uniformity  of  temperature  across  the 
top  plate  surface  is  determined  by  the  wire  spacing,  the  thermal  con- 
ductivity, the  plate  thickness  and  the  transpiration  flow  rate.  An 
analysis  of  this  problem  was  carried  out  by  Moffat  {1}  when  the  original 
HMT  Rig  was  designed.  Using  parameters  from  his  analysis,  the  tempera- 
ture non-uniformity  in  the  Roughness  Rig  plates  have  been  estimated  to 
be  on  the  order  of  0.01*F.  A key  parameter  in  this  analysis  is  the 
thermal  conductivity  of  the  porous  plate.  In  order  to  estimate  the 
thermal  conductivity  of  the  Roughness  Rig  plates,  the  ball  pack  was 
thermally  modeled  by  a series  of  resistances  representing  the  ball  layers 
and  nickel  connecting  spots.  Several  brazed  samples  were  examined  under 
a microscope  to  estimate  the  size  of  the  braze  connections  between  ad- 
jacent balls.  These  varied  somewhat  but  a conservative  size  for  the 
connecting  'spot'  diameter  was  .008  iDchea.  Using  this  as  an  average 
spot  size,  the  conductivity  of  the  ball  pack  becomes  30  BTU/hr  ft  *F. 

This  is  the  value  used  for  the  estimate  of  surface  temperature  uniformity. 
It  can  be  argued  that  this  is  a somewhat  conservative  estimate  since 
the  arrangement  of  the  balls  in  the  pack  are  such  that  a given  ball  is 
in  contact  with  6 surrounding  balls  on  the  same  layer,  but  with  only 
3 balls  each  in  the  balls  layers  above  and  below.  Since  the  conductivity 
of  the  ball  pack  is  controlled  by  the  size  and  number  of  contact  spots, 
thu  ball  pack  used  for  the  Roughness  Rig  has  a higher  apparent  thermal 
conductive  across  a given  ball  layer  than  it  does  between  ball  layers. 

C. 3 Thermocouple  Installation 

Plate  thermocouples  were  imbedded  to  a depth  that  located  their 
junctions  at  the  center  of  the  ball  layer  below  the  surface  ball  layer. 
Several  braze  samples  were  drilled  for  thermocouple  holes  and  then 
sectioned  to  examine  the  condition  of  the  ball  pack  around  the  ho'e. 

The  brazing  process  apparently  removes  all  the  temper  from  thn  OHFC 
copper  balls  and  they  become  very  soft.  It  proved  difficult  to  drill 
the  ball  pack  and  the  technicians  had  to  carefully  prepare  the  drill 
bit  used  for  this  operation  to  achieve  clean  holes.  The  soft  copper 
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actually  wiped  itself  at  the  hole  surface.  The  holes  appeared  to  be 
through  solid  material  when  examined  in  cross  section.  Upon  closer 
examination,  it  appeared  ae  if  the  Interstitial  regions  between  ball 
layers  were  covered  with  a very  thin  layer  of  copper  at  the  hole  surface 
so  the  hole  appeared  to  have  a solid  wall.  Further,  by  drilling  the 
holes  on  a milling  machine,  it  was  possible  to  locate  the  bottom  of  the 
hole  at  a ball  center  in  the  ball  layer  below  the  surface.  In  this 
location,  the  thermocouple  junction  was  >068  inches  below  the  top  of 
the  surface  layer  of  balls.  Bacause  of  the  high  conductivity  of  the 
ball  pack,  the  temperatures  recorded  by  these  thermocouples  were  taken 
to  be  the  plate  surface  temperatures.  The  thermocouples  were  installed 
by  partially  tilling  the  holes  with  epoxy  using  a hypodermic  needle, 
then  insert log  the  thermocouple  into  the  hole  until  It  was  felt  to  bot- 
tom out.  This  resulted  in  s small  collar  of  epoxy  squeezing  up  around 
the  thermocouple  at  the  plate  back  surface.  Every  Installation  was 
examined  from  the  front  surface  and  there  were  no  cases  where  bleeding 
through  of  the  epoxy  sround  the  thermocouple  holes  could  be  detected. 

C. 4 Plate  Permeability 

Uniformity  of  plate  permeability  is  an  important  requirement  for 
any  transpired  boundary  layer  experiment.  Non-uniformities  can  result 
in  variations  in  both  surface  temperature  as  well  as  local  transpiration 
rate.  Permeability  measurements  on  the  plates  used  in  the  smooth  plate 
HMT  rig  were  made  before  their  installation.  An  instrument  was  developed 
which  sealed  off  both  faces  of  the  piste  except  for  s 3/4  inch  diameter 
spot,  and  then  forced  a meterod  flow  of  air  through  that  spot.  The 
pressure  drop  across  the  plate  provide  the  measure  of  the  resistance 
to  flow  through  the  test  area.  With  this  device,  porosity  maps  were 
constructed  by  making  measurements  on  one  inch  centers  over  the  surface 
of  each  plate. 

Porosity  varistion  measureaants  on  the  Roughness  Rig  pistes  were 
made  after  the  plates  were  in  place  in  the  sub-assembly  castings,  with 
thermocouples  and  heater  wires  lnstallsd,  using  the  actual  rig  transpira- 
tion air  supply  system.  After  assembly  of  each  casting  was  complete,  It 
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was  Installed  on  the  rig  support  structure,  with  the  transpiration  air 
supply  connected  and  test  section  sidewalls  but  not  the  test  section  top 
xn  place.  A constant  current  hot-wire  anemometer  was  supported  on  a 
traversing  mechanism  mounted  on  the  sidewalls.  It  was  adjusted  to  a 
height  of  about  1/8  inch  above  the  plate  surface.  Output  from  the 
anemometer  was  run  directly  into  an  X-Y  plotter.  A resistance  bridge 
built  into  the  traversing  mechanism  drove  the  plotter  such  that  the  pen 
movement  across  the  graph  corresponding  approximately  to  probe  movement 
over  the  plate.  This  set-up  was  used  to  map  the  surface  velocity  of 
the  air  delivered  through  each  plate. 

Initial  results  from  these  tests  yielded  velocity  traces  which 
appeared  very  lumpy  at  high  blowing  rates.  The  non-uniformity  in  velocity 
was  several  percent  and  the  spanwise  period  of  the  velocity  disturbances 
appeared  to  be  several  ball  diameters.  A literature  survey  into  the 
subject  of  flow  through  porous  media  uncovered  several  interesting  papers 
which  seem  to  explain  the  phenomena  that  had  been  observed.  Two  of 
particular  interest  are  by  Bradshaw  [44]  and  Morgan  [45],  both  of  whom 
studied  the  flow  through  screens.  Both  of  these  studies  showed  the 
existence  of  a critical  velocity  in  the  flow  through  the  porous  material. 
Above  this  velocity,  the  flow  emerges  from  the  porous  media  as  a pattern 
of  jets.  These  jets  coalesce  into  random  groups  because  they  can  only 
entrain  fluid  from  each  other.  The  velocity  field  investigated  by 
Morgan  was  behind  a two  dimensional  grid.  It  showed  the  same  onset  of 
a lumpy  structure  with  a period  that  was  two  or  three  times  the  grid 
spacing.  To  confirm  that  the  "lumpy"  surface  velocity  measurements 
were  the  result  of  the  coalescence  of  the  jets  leaving  the  plate  surface, 
the  tests  were  repeated  over  a wide  range  of  surface  velocities.  Re- 
sults from  these  traces  are  shown  in  Fig.  2.9.  Here  the  surface  velocity 
traverses  have  been  superimposed  in  a single  graph.  The  top  graph  shows 
the  results  obtained  by  starting  at  a high  surface  velocity  and  reducing 
it.  The  top  curve  is  for  a surface  velocity  of  0.39  ft/sec.  It  shows 
the  lumpy  appearance  which  was  observed  over  all  the  plates.  The  mag- 
nitude of  the  non-uniformities  is  reduced  as  the  surface  velocity  is 
reduced  but  they  are  still  clearly  present  at  0.21  ft/sec.  The  velocity 
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is  much  smoother  at  0.18  ft/sec  and  at  lower  surface  velocities.  This 
indicates  there  is  a critical  velocity  at  or  near  0.18  ft/sec.  The 
lower  graph  in  the  figure  shows  results  obtained  when  the  surface 
velocities  are  increased  beginning  from  a low  value.  Here  smooth  traces 
are  obtained  up  to  0.24  ft/sec.  Above  this,  the  lumpy  traces  reappear 
and  persist  to  higher  velocities.  Apparently  the  onset  of  the  flow 
instabilities  have  a hysteresis-like  response  to  velocity  changes.  The 
jets  continue  to  coalesce  to  much  lower  surface  velocities  as  the  surface 
velocity  is  reduced  than  the  surface  velocity  at  which  coalescence  be- 
gins when  the  surface  velocity  Is  being  increased. 

As  a final  step  to  confirm  the  mechanism,  these  same  tests  were 
repeated  over  the  porous  plates  in  the  existing  HMT  rig.  Again  the 
flow  snowed  the  same  lumpy  structure,  but  the  velocity  traces  appeared 
much  finer  grained  than  had  been  observed  over  the  rough  surface.  This 
was  expected  siuce  the  jet  pattern  over  the  smooth  surface  porous  plates 
should  be  much  finer  to  start  with.  It  was  also  noted  that  the  onset 
of  flow  Instabilities  occur  at  surface  velocities  above  the  normal 
operating  range  of  the  rig. 

The  next  step  was  to  investigate  the  possible  effects  of  this  Jet- 
ting action  on  the  heat  transfer  experiments.  It  was  feared  that  the 
jetting  action  might  artificially  enhance  boundary  layer  heat  transfer. 

To  evaluate  this  possible  effect,  a series  of  eight  Stanton  number  runs 
with  blowing  were  made  at  a test  section  velocity  of  90  ft/sec.  At  this 
velocity,  it  was  possible  to  run  several  blowing  fractions  above  and 
below  the  critical  surface  velocity  observed  in  the  permeability  mapping. 
If  there  was  an  effect  on  heat  transfer  due  to  the  onset  of  a surface 
jetting,  it  should  clearly  be  seen  in  the  heat  transfer  data.  The  data 
was  taken  and  reduced  without  corrections  for  radiation  or  conduction 
losses.  Although  it  showed  some  scatter,  it  clearly  indicated  that 
there  was  no  abrupt  change  in  the  surface  heat  transfer.  Apparently, 
the  presence  of  the  boundary  flow  on  the  porous  surface  stabilized  the 
transpiration  flow  and  there  was  no  onset  of  instability  or  any  similar 
effect  that  can  be  seen  in  the  heat  transfer  data. 

The  porosity  maps  for  the  plates  were  completed  but  the  flow 
striatlons  could  be  eliminated  only  by  reducing  the  transpiration 


velocity  to  low  values.  At  these  low  levels , the  width  of  the  recorded 
line  was  equal  to  about  5Z  of  the  velocity  signal.  This  prevented  the 
careful  resolution  of  non-uniformities  that  had  been  hoped  for.  With 
that  reservation,  the  plate  porosity  appeared  entirely  uniform:  the 

recorded  traces  contained  no  discernible  evidence  of  non-uniformities. 

D.  Rig  Instrumentation 

D.l  Temperature  Instrumentation 

All  temperature  measurements  on  the  Roughness  Rig  are  made  using 
iron-constantan  (ISA  type  J)  thermocouples.  Many  precautions  have  been 
taken  to  reduce  the  spurious  OIF's  frequently  encountered  in  thermocouple 
circuits.  The  thermocouples  are  all  brought  together  at  a common  test 
console  zone  box  where  they  are  connected  to  rotary  thermocouple  switches 
for  read-out.  A diagram  of  the  thermocouple  circuit  is  shown  in  Fig. 

2.10.  The  thermocouples  were  made  sufficiently  long  so  that  the  thermo- 
couple wire  Itself  could  be  used  for  the  test  console  lead  wires.  To 
avoid  introducing  sharp  temperature  gradients  in  any  of  the  thermocouples, 
all  plate,  transpiration  air,  and  casting  thermocouples  were  thermally 
guarded  with  Polyflo  tubing.  The  test  console  zone  box  is  constructed 
from  1/2  inch  plywood,  lined  inside  with  1/16  inch  aluminum  plate  and 
insulated  outside  with  aluminum  foil  backed,  rock  wool  insulation.  All 
entry  ports  are  gasketed  and  all  busing  of  wire  connections  within  the 
zone  box  is  done  with  thermocouple  wire.  These  precautions  were  taken 
to  reduce  temperature  stratification  within  the  zone  box  and  to  minimize 
its  effect.  Thermocouples  are  mounted  at  opposite  ends  of  the  zone  box 
to  provide  a direct  measurement  of  the  zone  box  temperature  gradient. 

In  its  present  configuration,  temperature  gradients  greater  than  1/10°F 
have  never  been  observed.  The  entire  thermocouple  circuit  uses  a single 
ice-bath  reference  junction  which  is  also  thermally  guarded  in  Polyflo 
tubing  to  stretch  its  thermal  gradients.  Thermocouple  output  is  measured 
with  a Hewlett-Packard  Integrating  Digital  Voltmeter,  Model  2401C. 

The  most  difficult  temperature  measurement  problem  in  this  system 
is  the  plate  surface  temperature.  These  thermocouples  are  glued  into 
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close-f itcing  holes  In  the  bottom  of  the  plate.  The  wire  extending 
from  the  bottom  surface  Is  exposed  to  the  transpiration  air  which  is 
at  a different  temperature  than  the  plate.  This  tends  to  introduce 
an  error  in  the  measurement  of  the  plate  temperature.  To  estimate  this 
error  on  the  smooth  plate  HMT  Rig,  Moffat  [1]  developed  a thermal  model 
of  the  thermocouple  treating  it  like  a cylindrical  fin,  partly  exposed  eo 
the  plate  and  partly  exposed  to  the  transpiration  air.  The  problem  is 
complicated  by  the  fact  that  the  thermocouple  may  not  be  perfectly 
bonded  in  the  hole.  This  problem  is  modeled  by  introducing  a 'bond 
factor'  between  the  plate  and  thermocouple.  Because  of  the  small  thermo- 
couple wire  diameter  and  relatively  deep  immersion  into  the  plate, 
bonding  factors  can  be  taken  as  low  as  0.5  and  still  have  conduction 
errors  less  than  0.1°F  for  a 2Q#F  temperature  difference  across  the 
plate.  Subsequent  qualification  tests  on  the  plate  indicated  that  the 
thermocouple  error  due  to  imperfect  bonding  and  to  other  sources  are 
acceptably  small.  These  tests  are  discussed  in  Section  E of  this 
chapter. 

None  of  the  temperature  measuring  locations  on  the  Roughness  Rig 
are  so  restrictive  in  space  that  conduction  error  becomes  a problem. 

Ample  immersion  depth  is  available  everywhere. 

As  noted  earlier,  each  plate  contains  five  thermocouples  arranged 
in  a cross  pattern.  In  the  tests  performed  here,  these  thermocouples 
were  wired  in  parallel  in  the  instrument  console  zone  box  so  that  an 
average  of  their  collective  readings  were  taken.  Care  was  taken  to 
size  all  the  thermocouple  lengths  the  same  for  each  casting  so  if  they 
were  ganged  together,  they  would  give  a true  average  reading. 

D.2  Pressure 

Pressures  were  measured  on  a variety  of  manometers  and  transducers. 
A 3 inch  inclined  Merrlam  Manometer  was  used  to  measure  tunnel  static 
pressures  and  to  set  the  top  for  zero  pressure  gradient  conditions. 

This  same  manometer  was  used  to  record  mainstream  total  pressure  when 
test  section  velocities  were  below  100  ft/sec.  At  higher  mainstream 
velocities,  total  pressures  were  read  out  with  a transducer.  Two 
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St* than  unbonded  strain  gauge  differential  pressure  transducers  were 
used  for  velocity  profile  measurements;  a PM5  (pressure  range  0 to  0.5 
psi)  and  a PM97  (pressure  range  0 to  0.05  psi).  Both  units  were  equipped 
with  zeroing  bridges  and  individually  calibrated  in  the  Thermosciences 
Measurements  Center  against  a precision  30"  Merrlam  Micromanometer.  The 
H-P  2401  1DVM  equipped  with  an  external  quartz  crystal  oscillator  clock 
was  used  to  read  the  pressure  transducers. 

D.3  Flow  Rate 

Transpiration  flow  rates  for  ?ach  of  the  porous  plates  are  measured 
using  hot-wire  type  flow  meters  which  are  discussed  in  Appendix  B.  The 
flow  signal  from  these  meters  is  from  a differential  thermocouple  which 
is  read  out  with  the  H-P  2401  IDVM  using  the  external  clock  tc  extend 
signal  integration  time.  Before  each  flow  meter  reading  is  taken,  the 
flow  meter  heater  circuit  current  is  checked  to  insure  it  is  set  exactly 
at  the  calibration  value.  This  is  done  by  reading  the  voltage  drop 
across  a precision  Weston  1 amp  shunt,  again  with  the  IDVM. 

In  addition  to  the  flow  meter  signal,  transpiration  air  temperature 
for  each  supply  line  is  measured  so  the  appropriate  temperature  cor- 
rections can  be  made. 

D.4  Electric  Power  Measurement 

One  advantage  of  the  D.C.  power  supply  system  used  in  the  Roughness 
Rig  is  that  power  measurements  are  relatively  simple.  The  heater  voltage 
is  measured  directly  with  leads  which  are  attached  to  the  heater  wires 
just  as  they  leave  the  castings.  The  heater  ground  connection  is  made 
through  Individual  shunts  which  are  mounted  beneath  the  bus  bar  box. 
Measurement  leads  from  the  shunt  and  heater  are  taken  in  shielded  pairs 
to  the  same  Instrument  console  where  the  thermocouples  are  read,  but 
are  connected  in  a separate  * zoned'  selector  switch  station  for  read 
out.  The  heater  and  shunt  voltage  are  read  separately,  using  the  IDVM, 
and  plate  power  calculations  are  made  in  the  data  reduction  program 
using  the  individual  shunt  resistances  obtained  from  an  in-place  calibra- 
tion of  each  shunt.  The  plate  heater  power  supply  voltages  have  been 
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carefully  scoped  to  Insure  that  a constant  voltage  level  Is  applied 
to  each  plate  over  the  entire  power  range.  This  insures  that  the 
selective  sampling  of  the  1DVM  will  yield  a truly  representative  average 
measurement  for  the  voltage  across  the  heater  and  shunt. 

D. 5 Mainstream  Conditions 

For  all  runs,  mainstream  temperature  and  total-to-static  pressure 
were  measured.  In  addition,  mean  velocity  profiles  were  taken  at  sev- 
eral positions  along  the  test  section  length.  Mainstream  temperature 
was  measured  through  the  circuitry  already  described  with  a probe  con- 
structed using  .004  inch  iron-constantan  thermocouple  wire.  This  probe 
is  a fixed  position  version  of  the  traversing  probe  described  by  Kearney 
[6],  The  mainstream  total  pressures  were  measured  with  a Kiel-type 
probe  located  in  the  center  of  the  potential  flow  region  and  the  static 
pressures  were  taken  from  the  adjacent  wall  tap.  All  static  wall  taps 
were  .040  diameter  at  the  wall  plane  with  0.125  inch  diameter  tubing 
connections  outside  the  test  section. 

Mean  velocity  profiles  were  taken  using  a small  diameter  boundary 
layer  probe  mounted  in  a micrometer  driven  traversing  probe  holder. 

The  probe  holder  is  supported  from  an  instrument  sled  which  is  aligned 
with  the  instrument  ports  in  the  test  section  top  using  locating  pins 
and  is  attached  with  hold  down  bolts  to  the  side  walls.  A photograph 
of  the  probe  itself  is  shown  in  Fig.  2.11.  The  probe  stem  is  made 
from  1/8  inch  diameter  brass  tubing.  Into  the  brass  stem  is  soldered 
a .030  inch  stainless  steel  hypodermic  needle  which  has  been  bent  into 
a C-shape.  The  probe  mouth  has  been  flattened  down  to  a thickness  of 
.022  inches. 

In  use,  the  probe  was  lowered  until  it  was  in  visible  contact  with 
the  rough  surface.  This  could  be  confirmed  with  a resistance  reading 
from  a V0M  attached  to  the  probe  and  to  a thermocouple  in  the  plate 
the  probe  was  centered  over.  As  the  probe  was  traversed  away  from  the 
surface,  the  circuit  was  broken  when  the  probe  lifted  from  the  surface. 
This  method  was  ucmd  to  locate  the  first  reading  with  respect  to  the 
crest  of  the  surface  balls.  The  traversing  probe  was  checked  at  least 


once  during  each  run  by  comparing  It  with  a Kiel  probe  to  Insure  It 
was  not  clogged  with  dirt  particles  which  would  result  In  erroneous 
readings.  Care  was  also  taken  to  align  the  probe  before  every  set  of 
readings  to  avoid  error  introduced  by  probe  yaw. 

E.  Rig  Qualification  Tests 

The  Roughness  Rig  was  tested  in  detail  for  reliability  before  being 
approved  for  use.  There  were  three  types  of  qualification  tests  per-* 
formed:  tests  of  the  mainstream  condition,  instrument  system  qualifi- 

cation tests  and  energy  balance  tests  to  determine  validity  of  the  data 
reduction  program. 

E.l  Mainstream  Conditions 

Extensive  uniformity  and  stability  tests  were  conducted  on  the  test 
section  inlet  flow  by  M.  Crawford,  another  member  of  the  HMT  group.  The 
purpose  of  these  tests  was  to  measure  mainstream  uniformity  and  local 
free  stream  turbulence  intensity.  The  importance  of  a uniform  free 
stream  inlet  condition  is  to  insure  a two-dimensional  flow  in  the  test 
section.  Boundary  layer  skin  friction  and  enthalpy  thickness  are  inferred 
from  momentum  thickness  and  Stanton  Number  measurements  and  the  integral 
momentum  and  energy  equations.  A two-dimensional  boundary  layer  is 
necessary  if  these  techniques  are  to  be  used  in  reducing  the  experimental 
data. 

To  teat  mainstream  uniformity,  velocity  profiles  across  the  test 
section  inlet  were  taken  using  a differential  Kiel  probe  technique. 

The  differential  traversing  employs  two  Identical  probes  connected  to 
a differential  measuring  instrument.  This  technique  allows  separation 
of  space-wise  and  time-wise  variation  in  flow,  to  a first  order  at  least. 
If  the  velocity  field  can  be  represented  as  a product  of  two  functions; 
one  varying  slightly  in  space,  the  other  varying  in  time,  then  differen- 
tial traversing  (using  probes  of  the  same  time  constant)  will  display 
only  space-wise  variations  on  the  flow  field.  Two  test  section  velocities 
were  examined;  the  lowest  velocity,  32  ft/sec  and  maximum  test  section 
velocity,  242  ft/sec.  Measurements  were  taken  at  135  positions  in  the 
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test  section  Inlet  plane  at  each  of  the  two  velocities.  The  study 
showed  velocity  defects  no  greater  +0,162  at  the  low  test  section 
velocity  and  +0.102  at  the  higher  test  section  velocity. 

Free  stream  turbulence  intensity  measurement  was  also  made  with  a 
DISA  5SD05  constant  temperature  anemometer  read  through  a DISA  55D15 
linear izer.  The  velocity  fluctuations  about  the  mean  velocity  were 
measured  on  a Thermosystems  RMS  voltmeter,  Model  1060.  Measurements 
were  made  at  the  same  two  test  section  velocities  and  at  a total  of  120 
positions  in  the  test  section  inlet  plane.  These  results  showed  a 
very  uniform  turbulence  intensity  over  the  entire  inlet  plane  with  an 
intensity  level  of  about  0.4*  at  both  test  section  velocities.  This 
was  higher  than  expected  but  still  below  levels  measured  in  the  existing 
HMT  Rig.  The  design  objective  for  the  screen  pack  was  a free  stream 
turbulence  intensity  level  of  0.12  to  0.22.  This  higher  than  expected 
level  does  not  limit  the  rig's  performance. 

An  additional  check  on  the  two-dimensionality  of  the  test  section 
boundary  layer  was  made  by  measuring  its  uniformity  in  the  lateral 
direction  at  the  end  of  the  test  section.  Lateral  holes  are  provided 
for  this  purpose  in  the  test  section  top  over  plate  ^3.  Velocity  pro- 
files were  taken  at  the  test  section  center  and  at  two  positions  on 
either  side  of  the  center  for  a tunnel  velocity  at  90  ft/sec.  Momentum 
thicknesses  tn>m  these  profiles  are  plotted  in  Fig.  2.12,  the  variation 
can  be  seen  to  be  no  more  than  +32. 

E.2  Instrument  System  Qualification  Tests 

To  qualify  the  plate  temperature  measuring  system  a scheme  was 
used  that  had  been  developed  for  use  on  the  smooth  plate  HMT  Rig. 

Rather  than  calibrate  the  Individual  thermocouples,  they  were  calibrated 
in  place  for  each  plate.  This  was  accomplished  using  a plenum  box  which 
could  be  placed  over  each  test  plate  in  such  a way  as  to  collect  all 
the  air  passing  through  the  center  6-inch  span.  The  plenum  box  was 
surrounded  by  four  similar  guard  chambers  and  was  equipped  with  a 
guard  heater  system  so  that  the  walls  of  the  box  could  be  maintained 
at  the  same  temperature  as  the  air.  The  entire  unit  was  lined  with 


aluminized  Mylar  and  Insulated  with  1/2  inch  felt.  A mixing  section 
was  Installed  in  the  plenum  box  so  that  the  air  entering  from  the  test 
plate  would  be  mixed  thoroughly  before  its  temperature  was  measured  with 
two  independent  thermocouples.  The  guard  chamber  surrounding  the  plenum 
box  was  equipped  with  adjustable  area  vent  holes  so  all  chambers  had 
the  same  static  pressure  to  prevent  cross  flows.  The  bottom  surfaces 
of  the  box,  where  it  contacts  the  plates  are  covered  with  closed  cell 
foam  rubber  to  further  prevent  cross  flow  between  chambers.  In  opera- 
tion, the  guard  heater  system  power  is  adjusted  until  thermocouples  in 
the  box  walls  indicate  that  the  box  is  at  tha  same  temperature  as  the 
mean  mixed  air  temperature,  as  measured  by  the  mixing  chamber  thermo- 
couples. This  temperature  is  then  compared  to  the  average  of  the  five 
thermocouples  in  the  plate  beneath  the  plenum  box.  These  tests  were 
conducted  over  a range  of  blowing  rates  from  5 to  20  CFM  and  over  a 
range  of  air  temperature  differences  from  zero  (no  plate  power)  to  20*F. 
Over  this  wide  range  of  temperature  and  flow  conditions  none  of  the 
plate  average  thermocouple  readings  differed  from  the  mixing  chamber 
thermocouple  readings  by  more  than  0.2°F  (6  microvolts).  From  these 
tests  it  was  concluded  that  the  plate  temperature  measurements  as  re- 
corded by  the  five  plate  thermocouples  was  satisfactory  and  errors  due 
to  imperfect  bonding  of  the  thermocouples  and  other  sources  are  accept- 
ably small. 

While  these  tests  were  in  progress,  thermocouple  traverses  were 
made  under  the  plates,  just  below  the  honeycomb  using  a special  probe 
installed  through  fittings  in  the  side  of  the  castings  in  the  center 
of  each  plate.  These  traverses  shoved  a nearly  uniform  temperature 
distribution  beneath  the  plate.  Typically  these  profiles  showed  higher 
readings  at  the  plate  center  by  2 or  3 microvolts  than  at  3 inches  to 
either  side  of  center. 

A final  series  of  tests  were  performed  using  the  transpiration  air 
system  with  tne  top  removed,  but  without  plate  power.  These  teste  were 
run  with  both  hot  and  cold  transpiration  air,  by  either  by-passing  or 
using  maximum  cooling  in  the  transpiration  air  heat  exchanger.  In 
addition,  either  maximum  or  no  casting  cooling  water  was  used.  The 
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purpose  of  these  tests  wss  to  provide  s maximum  amount  of  thermal  dis- 
tortion, without  plate  power,  to  insure  that  the  temperature  measuring 
system  could  properly  record  these  unusual  renditions.  Operating  in 
this  mode,  the  plate  thermocouples  can  be  used  to  check  the  single, 
centrally  located  air  thermocouple  located  just  below  the  plate.  With- 
out plate  power,  the  five  plate  thermocouples  read  the  average  tempera- 
ture of  the  transpiration  air.  In  cases  where  there  is  a temperature 
mismatch  between  the  casting  and  transpiration  air,  the  air  thermocouple 
is  subject  to  error.  It  was  discovered  that  the  casting  temperature  had 
a much  more  dominant  effect  on  the  transpiration  air  temperature  than 
had  been  thought.  When  the  casting  inlet  air  was  colder  than  the  casting, 
the  five  plate  thermocouples  Indicated  a higher  reading  than  the  single 
air  thermocouple  behind  the  plate.  The  opposite  was  true  when  air  enter- 
ing the  casting  was  hotter  than  the  casting.  Apparently  the  air  next 
to  the  casting  is  heated  or  cooled  by  the  casting  and  enters  the  plates 
with  a temperature  gradient.  This  gradient  was  not  observed  in  the 
traverses  taken  beneath  the  plates  since  these  were  taken  in  the  longi- 
tudinal direction.  The  measurement  error  due  to  the  inlet  temperature 
gradient  is  related  to  the  mismatch  between  the  casting  temperature  and 
the  temperature  of  the  transpiration  air  as  it  enters  the  casting. 

The  temperature  of  the  transpiration  air  is  measured  at  the  header 
box.  Just  upstream  of  the  flow  meters  to  provide  the  temperature  cor- 
rections for  the  flow  meter  readings  but  not  at  the  inlet  to  the  casting. 
To  estimate  the  air  temperature  at  the  casting  inlet,  the  transpiration 
system  was  modeled  as  a heat  exchanger  between  the  air  and  the  laboratory. 
Thermocouples  were  Installed  at  the  casting  inlet  of  four  pistes  to  check 
this  model  and  good  agreement  wea  obtained.  Using  the  casting  inlet 
temperatures  estimated  with  this  model,  a correction  term  was  formulated 
for  the  air  thermocouples.  The  error  was  sssumed  to  vsry  liuesrly  with 
the  inlet  air-to-csstlng  temperature  difference.  This  correction,  while 
relatively  small  aided  in  closure  of  system  energy  balance. 

E. 3 Energy  Balance  Tests 

Energy  balance  tests  were  conducted  to  establish  the  validity  of 
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Che  experimental  system.  The  system  Includes  not  only  the  Roughness  Rig 
and  Its  measuring  system  but  the  procedure  used  to  reduce  data  from  the 
rig.  Experimental  Stanton  numbers  are  obtained  by  taking  plate  power, 
subtracting  off  the  lossea  and  the  energy  carried  away  by  the  transpira- 
tion air,  and  dividing  by  the  plate-to-free  stream  temperature  difference. 
In  equation  form  this  Is  given  by: 


St  - 


(plate  power)  - m"c  (Tw_Tt)  “ (losses) 
G c(Tw-Tw) 


(2-1) 


A key  part  of  the  data  reduction  program  is  the  modeling  used  to 
evaluate  energy  losses.  Depending  on  blowing  fractions  and  plate  temper- 
atures, these  losses  represent  from  1Z  to  10Z  of  the  total  power  supplied. 
The  loss  mechanisms  modeled  are:  radiation  from  the  top  and  bottom  sur- 

faces of  the  plates,  conduction  between  the  plates  and  the  casting  through 
the  phenolic  support  webs  and  conduction  through  the  stagnant  air  beneath 
the  plates  when  there  la  no  transpiration.  The  purpose  of  the  energy 
balance  tests  was  to  evaluate  the  modeling  used  for  energy  losses. 

Two  kinds  of  energy  balance  tests  were  conducted  on  the  Roughness 
Rig.  Each  of  these  tests  provided  a direct  measure  of  the  ability  of  the 
rig's  measuring  system  to  produce  zero  values  of  measured  heat  transfer 
when  teat  conditions  were  such  that  there  was  zero  heat  transfer. 

In  the  first  energy  balance  mode,  the  wind  tunnel  was  operated  with- 
out transpiration  and  without  main  stream  cooling.  In  this  mode,  the 
free  stream  temperature  equilibrated  near  100*F.  Plate  power  was  adjusted 
to  match  the  plate  temperature  to  the  free-stream  temperature.  In  this 
operating  condition,  the  plate  power  exactly  equaled  plate  losses.  These 
losses  were  due  to  radiation  from  the  back  of  the  plate  and  conduction 
to  the  casting  and  to  the  air  beneath  the  plate.  These  tests  were  run 
with  no  casting  cooling,  to  minimize  conduction  losses,  and  then  again 
with  maximum  casting  cooling  to  magnify  conduction  losses.  Since  the 
plate  and  free-stream  temperatures  were  the  same,  there  was  no  heat 
transfer  between  the  plate  and  the  free-stream.  The  difference  between 
measured  plate  power  and  calculated  losses  for  the  observed  conditions 
represent  a potential  Stanton  number  error.  Data  taken  from  this  energy 
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balance  run  ware  reduced  with  the  data  reduction  program  and  the  net 
Stanton  number  error  waa  evaluated  using  the  energy  balance  difference 
and  a typical  plate-to-free  stream  temperature  difference: 


St 


(plate  power)  - (calculated  losses) 

fled  -T) 

w 06  typ 


(2-2) 


The  typical  temperature  difference  used  waa  25*F  which  is  represen- 
tative of  toe  Roughness  Rig  test  conditions.  Figure  2.13  shows  a plot 
of  Stanton  number  error  versus  distance  along  the  teat  section.  There 
was  some  difficulty  in  closing  the  energy  balance  in  the  first  few  plates, 
attributed  to  an  inlet  thermal  boundary  layer.  With  the  tunnel  i .rating 
hotter  than  ambient,  a thin  thermal  boundary  layer  of  cool  air  next  to 
the  duct  walls  was  unavoidable.  Since  the  plate  and  free-stream  temper- 
atures were  equal,  the  first  few  plates  removed  the  inlet  thermal  boundary 
layer  and  acted  as  guard  heaters  for  the  rest  of  the  test  section. 

In  the  second  energy  balance  mode,  the  tunnel  was  operated  with 
transpiration  flow  only  and  with  the  test  section  top  removed  so  that  the 
transpiration  flow  moved  directly  upwards.  In  this  operating  condition, 
there  was,  again,  no  surface  heat  transfer.  All  the  plate  power  was 
taken  up  either  in  losses  or  by  the  transpiration  flow.  These  blowing 
energy  balances  were  conducted  over  a wide  range  of  transpiration  rates 
and  air-to-plate  temperature  differences,  again  to  demonstrate  the  ade- 
quacy of  the  data  reduction  program  to  predict  the  losses  and  achieve 
energy  balance  closure  at  all  operating  conditions.  For  this  case  the 
potential  Stanton  number  error  was  given  by: 


ASt  - 


(plate  power)  - m"c  (Tw*-Tt)  - (calculated  losses) 


Gtyp  c(VT->typ 


(2-3) 


Blowing  energy  balance  runs  were  made  at  2,  8,  15  and  22  CFM  which 

covers  the  range  over  which  the  Roughness  Rig  experiments  were  conducted. 

As  before,  the  typical  temperature  difference,  (Tw-T—)  • was  taken 

at  25”F.  The  typical  mass  flow,  G.  was  chosen  so  that  the  combination 

typ 

of  it  with  the  transpiration  rate  would  give  a blowing  fraction  of  .008, 
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the  «?»<«■■»  tested.  This  was  done  for  all  but  the  lowest  transpiration 

rate.  The  value  of  G_  uaed  for  both  2 CFM  and  8 CFM  corresponded  to 

typ 

the  m<n <jBum  test  section  velocity.  A plot  of  the  Stanton  mother  error 
from  these  tests  versus  distance  along  the  test  section  is  shown  in  Fig. 
2.13.  As  can  be  noted,  the  Stanton  number  error  was,  in  almost  every 
case,  within  +.0001  Stanton  nuaber  units.  In  those  few  cases  where  the 
error  exceeds  this  value,  subsequent  energy  balance  teats  did  close  to 
within  +.0001  level.  This  has  led  us  to  believe  that  the  Stanton  nuaber 
data  are  reliable  to  within  +.0001  Stanton  nuaber  units  over  the  range 
of  conditions  tested. 

E.4  Verification  of  the  Data  from  Plate  24 

Despite  the  excellent  closure  that  was  obtained  in  the  energy  bal- 
ance tests,  an  anomaly  was  noticed  in  the  heat  transfer  data  from  plate 
24.  At  conditions  of  high  blowing  and  high  test  section  velocity,  plate 
24  exhibited  a muen  lower  Stanton  nuaber  than  its  neighbor,  piste  23. 

This  behavior  waa  entirely  repeatable  and  in  fact  was  visible  in  earlier 
data  taken  when  the  rig  was  first  being  shaken  down,  but  was  not  as  ap- 
parent then.  After  carefully  confirming  the  flow  meter  calibration  and 
the  teaperaturc  measuring  system  it  was  concluded  that  the  reduced  Stanton 
nuaber  was  real  and  not  the  result  of  trouble  in  the  measuring  system. 

A review  of  the  possible  sources  of  this  problem  led  to  the  conclusion 
that  it  was  related  to  the  diffuser,  directly  downstream  of  plate  24. 

It  seeattd  likely  that  at  high  blowing  rates,  when  there  was  injection  of 
low  momentum  fluid  into  the  boundary  layer,  the  exit  diffuser  experienced 
stall.  This  stall  probably  began  right  at  the  diffuser  inlet  but  may 
have  propagatad  upstream  to  plate  24  itself.  The  presence  of  the  'stall 
bubble'  could  disturb  the  flow  in  the  region  near  plate  24  enough  to  re- 
duce the  plate  heat  transfer.  The  exit  diffuser  was  designed  to  operate 
without  stall,  however  it  is  known  that  the  inlet  velocity  profile  has 
an  important  effect  on  diffuser  performance.  A thickened  boundary  layer 
due  to  blowing  could  result  in  stall  that  would  otherwise  uot  have  been 
there.  To  confirm  this  possibility  a temporary  extension  of  the  test 
section  surface  was  installed  in  the  first  stage  of  the  diffuser.  This 
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extension  produced  e step  In  the  flow  path  at  the  inlet  to  the  second 
stage  of  the  diffuser  and  reduced  the  overall  performance  of  the  diffuser, 
but  did  remove  the  influence  of  the  diffuser  from  plate  24.  The  tunnel 
was  then  run  at  240  ft/sec  with  a blowing  fraction  of  .002,  without  the 
step.  After  reaching  equilibrium,  data  were  recorded  and  then  the  dif- 
fuser insert  was  installed  and  the  run  repeated.  The  insert  reduced  dif- 
fuser performance  so  that  the  test  section  velocity  only  reached  226  ft/ 
sec.  Since  the  flow  meter  settings  were  unchanged,  a slightly  higher 
blowing  fraction  resulted.  This  reduced  the  Stanton  number  slightly  be- 
low the  previous  data,  but  the  plate  24  Stanton  number  was  significantly 
increased.  It  was  still  slightly  low  (.0001  Stanton  number  units)  but 
much  improved  over  ths  Stanton  number  without  the  Insert.  Figure  2.14 
shows  the  Stanton  number  data  from  these  two  teats. 

Because  of  this  difficulty,  the  plate  24  data  has  not  been  shown 
for  data  runs  at  high  velocity  and  blowing  fraction  but  is  listed  with 
the  tabulation  of  experimental  results  in  Appendix  G,  for  completeness. 

F.  Data  Reduction  Program 

In  reducing  data  taken  from  the  Roughness  Rig  experiments,  three 
separate  data  reduction  programs  were  used,  two  for  the  hydrodynamic  data 
and  one  for  the  heat  transfer  data. 

The  first  data  reduction  program  was  used  for  the  boundary  layer 
mean  velocity  profiles  and  converted  the  pitot  tube  pressure  readings  to 
velocity.  It  then  executed  a simple  trapezoidal  rule  numerical  integra- 
tion to  get  momentum  and  displacement  thicknesses  and  finally  shape  factor. 
The  second  data  reduction  program  used  a least  square  curve  fit  to  the 
momentum  thickness  data  for  the  entire  plate  length  to  determine  the 
skin  friction,  using  the  integral  momentum  equation  for  the  two-dimensional 
boundary  layer.  It  was  assumed  that  over  the  relatively  small  momentum 
thickness  range  of  a single  run,  the  momentum  thickness  varied  as  a simple 
power  law  with  distance  along  the  test  section.  The  equation  which  was 
curve  fit  was: 

e - a(x-x  )b  (2-4) 

o 
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where  xQ  is  the  boundary  layer  virtual  origin  and  a , b ar.d  x are 
determined  by  the  curve  fit.  The  procedure  used  was  that  of  'weighted 
residuals'  as  described  in  Chapter  18  of  Scarborough  [46] . Since  no  trip 
was  used  in  the  Roughness  Rig  experiments  both  laminar  and  turbulent  re- 
gions existed  inside  the  test  section.  Only  momentum  thickness  measure- 
ments taken  in  the  fully  turbulent  boundary  layer  were  used  in  the  fit. 

The  Stanton  number  data  was  then  examined  to  ensure  that  each  of  the  vel- 
ocity profiles  used  had  been  in  the  fully  turbulent  boundary  layer.  To 
qualify  this  procedure  for  predicting  skin  friction,  the  momentum  thick- 
ness measurements  taken  by  Andersen  [8}  and  by  Simpson  [2]  were  used  as 
a trial  data  set.  The  skin  friction  data  reported  by  Simpson  were  based 
on  his  'best  estimate'  using  several  evaluation  methods.  Skin  friction 
predictions  by  Andersen  were  made  using  turbulent  shear  measurements  taken 
in  the  boundary  layer,  then  extrapolated  to  the  wall.  Figure  2.15  com- 
pares skin  friction  predictions  made  by  the  present  curve  fitting  method 
to  the  values  originally  reported.  For  both  data  sets  the  curve  fitting 
method  yields  a reasonably  good  check  of  the  skin  friction  values  that 
Simpson  and  Andersen  found  by  independent  means.  Since  both  of  these  data 
reduction  programs  are  straightforward  applications  of  well  known  princi- 
pals, listings  of  the  programs  are  not  given. 

The  data  reduction  program  for  the  heat  transfer  data  is  more  com- 
plicated since  it  must  convert  data  from  the  Roughness  Rig  into  Stanton 
numbers.  A flow  diagram  for  the  data  reduction  program  is  shown  in  fig. 
2.16  and  a Fortran  listing  is  given  in  Appendix  F.  The  program  has 
several  features.  It  containu  the  individual  calibrations  for  each  flow 
meter  and  perform?  .i  logarithmic  interpolation  to  obtain  the  flow  from 
each  flow  meter  re*,  tng,  based  on  the  calibration  data  points.  Thermo- 
couple readings  are  converted  into  temperatures  by  interpolating  in  tables 
taken  from  National  Bureau  of  Standards  Circular  #561.  Plate  voltage  and 
shunt  voltages  are  converted  into  power  using  calibrated  resistance  values 
for  each  shunt  and  correcting  for  excess  heat  release  due  to  the  turn- 
around wrappings  of  the  heater  vires  at  the  plate  ends.  It  also  estimates 
conduction  and  radiation  losses  from  each  plate  based  on  plate,  casting 
and  other  rig  temperatures.  The  program  also  converts  manometer  readings 
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Co  pressure  levels  and  the  free-stream  thermocouple  reading  into  the  free- 
stream  total  temperature  using  a 0.36  recovery  factor  for  the  thermocouple 
probe.  The  program  has  tw«-  output  options:  the  summary  output  shown  In 

Appendix  E and  a more  detailed  output  with  additional  diagnostic  Informa- 
tion. 


G.  Experimental  Uncertainties 

To  investigate  the  possible  effects  of  measurement  uncertainty  on  the 
experimental  Stanton  number,  the  following  procedure  was  used.  The  data 
reduction  program  was  used  to  calculate  the  derivative  of  each  experimental 
variable  with  respect  to  the  Stanton  number.  This  was  done  by  sequentially 

varying  each  variable  by  the  amount  of  its  uncertainty  as  input  to  the 

data  reduction  program  and  calculating  the  change  in  Stanton  number  which 
resulted.  Assuming  that  these  uncertainties  are  all  independent  they  can 
be  combined  by  the  procedure  of  Kline  and  McClintock.  [47]. 

The  results  of  this  uncertainty  analysis  is  shown  in  Fig.  2.17  for 

all  test  section  velocities,  over  the  range  of  blowing  in  these  experi- 

ments. As  can  be  seen,  the  lowest  test  section  velocity,  at  high  blowing 
is  the  most  affected  by  the  propagation  of  uncertainties  in  these  experi- 
ments. The  following  b^slc  uncertainty  intervals  were  assumed  for  the 
input  values: 


Plate  Voltage 
Shunt  Voltage 
All  Temperatures 
Stagnation  Pressure 
Flow  Meter  Readings 


1 millivolt 
10  microvolts 
10  microvolts 
.002  inches  of  water 
25  microvolts 


It  can  be  concluded  from  this  study  that  except  at  high  blowing  the 
uncertainty  in  Stanton  numbers  due  to  random  errors  is  less  than  the  +.0001 
Stanton  number  units  cited  earlier  based  on  the  energy  balance  tests. 

Anotht  area  of  experimental  uncertainty  associated  with  these  ex- 
periments is  associated  with  the  velocity  profile  measurements.  A problem 
which  is  common  to  all  rough  surface  boundary  layers  is  the  location  of 
the  origin  of  the  velocity  profile.  Profile  measurements  begin  with  the 
probe  resting  on  top  of  the  rough  surface.  The  location  of  this  data 
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point  with  respect  to  the  surface  is  known,  but  examination  of  the  pro- 
file data  indicates  that  the  origin  of  the  profile  is  somewhere  below 
the  ball  crests.  Several  authors  have  discussed  systematic  procedures 
for  determining  'origin  error'.  What  has  been  done  in  the  present  work, 
is  to  examine  each  profile  on  an  individual  basis  and  extrapolate  s*  »o 
what  appears  to  be  a reasonable  choice  for  its  velocity  origin.  Generally 
this  turned  out  to  be  .010  inches  below  the  ball  crests.  The  bounds  on 
this  number  would  be  from  zero  (the  profile  origin  at  the  ball  crests) 
to  a maximum  of  .025  inches,  the  profile  origin  at  half  of  the  ball 
diameter.  With  these  bounds,  it  seems  unlikely  that  the  experimental  un- 
certainty of  the  profile  origin  could  exceed  £.010  inches.  To  investi- 
gate the  sensitivity  of  the  integral  boundary  layer  parameters  to  an  un- 
certainty in  the  profile  origin  of  +.010  inches,  the  velocity  profile 
data  reduction  program  was  used.  A portion  of  the  profile  data  was  re- 
duced again  with  a .010  inch  shift  in  the  profile  origin.  It  was  not 
surprising  to  find  that  the  displacement  thickness  was  very  sensitive  to 
the  shift  while  the  momentum  thickness  was  not  sensitive.  There  was  no 
detectable  effect  of  test  section  velocity:  this  is  not  surprising  since 

the  boundary  layers  were  all  very  similar  at  all  test  section  velocities. 
A more  important  parameter  was  the  thickness  of  the  boundary  layer.  Near 
the  front  of  the  test  section  where  the  probe  could  not  penetrate  as 
deeply  into  the  thin  boundary  layer,  an  origin  error  of  .010  inches  pro- 

ducted  errors  on  the  order  of  10%  in  the  boundary  layer  displacement 

thickness  and  1%  in  the  momentum  thickness.  At  the  end  of  the  test  sec- 
tion, where  the  boundary  layer  was  thicker,  the  uncertainty  in  profile 
origin  produced  errors  of  7%  to  8%  in  displacement  thickness  and  again, 

1%  in  momentum  *■  hickness.  From  these  results,  we  have  concluded  that  the 
experimental  uncertainty  in  the  profile  virtual  origin  results  in  un- 
certainties on  the  order  of  10%  in  the  boundary  layer  displacement  thick- 
ness and  1%  in  the  momentum  thickness. 

Some  additional  work  has  been  done  to  translate  the  momentum  thick- 
ness uncertainty  into  uncertainty  in  the  skin  friction  prediction.  This 

can  be  done  by  varying  the  inputs  to  the  program  which  uses  a fit  to  the 

-''men turn  thickness  data  to  predict  skin  friction.  The  results  from 
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study  have  been  Inconclusive.  The  effect  of  the  momentum  thickness  un- 
certainty depends  on  how  many  other  points  used  by  the  curve  fit  remain 
unchanged  and  even  on  which  point  is  perturbed.  Another  very  important 
parameter  is  the  blowing  fraction,  since  it  is  subtracted  from  the  mo- 
mentum thickness  before  the  skin  friction  fit  is  made.  A general  result 
which  was  obtained  is  that  a 1Z  variation  in  momentum  thickness,  for  the 
non-blown  boundary  layer  could  produce  a varation  of  8Z  to  10Z  in  skin 
friction  in  a boundary  layer  with  5 or  more  velocity  profiles  as  basis 
for  its  least  square  fit.  In  many  cases,  the  variation  was  less  and  in 
cases  with  high  blowing  the  variation  was  greater.  Because  of  the  many 
difficulties  Inherent  in  taking  the  derivative  of  experimental  data 
points,  it  Is  felt  that  no  further  claims  about  the  accuracy  of  the  skin 
friction  data  generated  by  fitting  the  momentum  thickness  measurements 
was  Justified.  Experimental  work  now  in  progress  on  the  Roughness  Rig 
should  produce  much  better  skin  friction  data,  with  much  less  experimental 
uncertainty. 
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Fig.  2.2  Photograph  of  Roughness  Rig 
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Fig-  2.3  Photograph  of  aluminum  casting 


Fig.  2.5  Photograph  of  clanping  arrangement  Fig.  2.6  Photograph  of  the  back  of  the  casting 

used  for  gluing  the  plate  and  casting  assembly 

assembly 
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Fig.  2.9  Velocity  distribution  above  the  Roughness  Rig  surface 
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Fig.  2.10  Roughness  Rig  thermocouple  schematic 


Fig.  2.11  Photograph  of  the  Pitot  probe 


from  iiuuifl  oviitt-r  llttu  (liirlitiv) 


Fig.  2.12  Momentum  thickness  distribution  across  the  tunnel  at  plate  23 
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Fig.  2.14  Stanton  number  versus  x-Reynolds  number  for  u» 
F - 0.002  (standard)  and  for  u,*,  ■ 226  fps,  F 
(with  an  insert  in  the  diffuser) 


o Simplon' • date  i:j  utj  - 42.6  ft/mec* 

— Predicted  from  fit  to  Simpson's 
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Fig.  2.15  Comparison  of  stein  friction  predic 
tions  to  data  by  Andersen  [8]  and 
Simpson  [2] 
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EXPERIMENTAL  RESULTS 

Measurements  taken  in  these  experiments  can  be  divided  into  two 
types  — Stanton  number  data  and  mean-velocity  profile  data  which  were 
used  to  infer  boundary  layer  momentum  thickness  and  skin  friction. 

The  range  of  conditions  for  the  Roughness  Rig  tests  can  be  summa- 
rized as  in  Table  3-1.  The  data  taken  can  be  organized  into  the  follow- 
ing blocks:  fully  developed  turbulent  boundary  layer  skin  friction  and 
Stanton  number  without  blowing,  the  effects  of  blowing  on  rough  surface 
Stanton  number  and  friction  factor,  and  rough  surface  transition  with 
and  without  blowing.  The  presentation  of  the  data  will  be  organized  in 
the  same  manner. 

A,  Rough  Surface  Boundary  Layer  Data  without  Blowing 

Rough  surface  boundary  layer  friction  factors  and  Stanton  numbers 
are  plotted  against  x- Reynolds  number  in  Fig.  3.1  for  each  of  the  five 
test-section  velocities  used.  It  is  apparent  that  there  is  a roughness 
effect  and  that,  in  these  coordinates,  the  friction  factor  and  the  Stan- 
ton number  increase  with  test  section  velocity.  This  is  an  important 
difference  from  smooth  plate  data,  which  do  not  show  a velocity  effect 
in  these  coordinates.  Smooth  plate  data  for  air  in  this  range  of  x- 
Reynolds  number  would  correlate,  for  all  test-section  velocities,  with 
equations  of  the  form: 

Smooth  cf/2  - .0295  Re”0,2  , (3-1) 

Smooth  St  - .0295  Pr“0,4  Re-0,2  - .034  Re“0,2  . (3-2) 

x x 

For  comparison,  these  two  corrections  are  shown  in  Fig.  3.1.  Only  the 
32  fps  data  seem  to  be  approaching  the  smooth  plate  correlations. 

With  transition  occurring  at  different  locations  for  each  velocity, 
it  might  be  argued  that  the  failure  to  correlate  is  due  to  different 
locations  of  the  boundary  layer  'virtual  origin'.  To  avoid  this  problem, 
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boundary  layer  data  are  frequently  presented  In  local  coordinates. 

Fig.  3.2  shows  friction  factor  and  Stanton  number  as  functions  of  local 
momentum  thickness  Reynolds  number  and  enthalpy  thickness  Reynolds  num- 
ber, respectively.  In  these  coordinates,  too,  the  effect  of  the  rough 
surface  has  been  tc  separate  the  data  taken  at  different  test  section 
velocities.  Also  shown  are  typical  correlations  for  a smooth  plate  in 
local  coordinates: 

Smooth  cf/2  - .0128  Re^*25  (3-3) 

Smooth  St  - .0128  Pr-0*5  Ra“0,25  - .0153  Re"*25  (3-4) 

It  is  apparent  that  boundary  layer  thickness  Reynolds  number  is  not 
a sufficient  descriptor  for  the  rough  surface  data.  At  any  boundary 
layer  thickness  Reynolds  number,  higher  velocities  yield  higher  values  of 
friction  factor  and  Stanton  number.  While  the  data  at  32  fps  lie  close 
to  the  smooth  correlations,  at  the  highest  test  section  velocity  the 
friction  factor  is  nearly  twice  the  expected  smooth  plate  value.  At  that 
velocity  the  Stanton  number  is  up  by  about  one  and  one-half  times.  This 
behavior  is  similar  to  that  reported  by  Nunner  [23]  for  rough  pipes.  Uis 
results,  given  by  Eqn.  (1-1),  show  the  lncreaee  in  heat  transfer  due  to 
roughness  approximately  equal  to  the  square  root  of  the  increese  in  skin 
friction. 

The  present  data  are  clearly  affected  by  the  eurface  roughness.  As 
already  noted  in  Chapter  I,  the  roughness  Reynolds  number,  ReT  - u^k^/y, 
is  often  used  to  categorize  the  behavior  of  rough  surfaces.  Values  less 
than  5 to  10  define  'smooth*  behavior,  while  between  10  and  70  the  sur- 
face is  described  as  'transitionally  rough',  and  above  70  it  is  referred 
to  as  'fully  rough'.  The  roughness  Reynolds  number  range  for  these  tests 
is  shown  in  Table  3-1.  The  32  fps  data  lie  between  24  and  29,  in  the 
'transitionally  rough'  range,  while  all  other  velocities  have  'fully 
rough'  roughness  Reynolds  numbers.  For  these  calculations,  ka>  the 
equivalent  sand-grain  roughness  was  taken  as  0.625  times  the  ball  diam- 
eter, yielding  an  equivalent  sand  grain  roughness  of  .031  inches.  This 
multiplier  Is  the  value  recommended  by  Schlichting  [51]  for  densely 
packed  spheres. 
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The  effects  of  roughness  are  also  evident  In  the  rough  surface 
velocity  profiles.  Fig.  3.3  shows  velocity  profiles  from  these  experi- 
ments at  different  test  section  velocities.  In  Fig.  3.4,  velocity 
profiles  are  shown  for  one  velocity  and  several  different  stations  along 
the  test  section,  along  with  a typical  velocity  profile  for  a smooth 
surface.  The  rough  surface  velocity  profiles  show  the  characteristic 
offset  or  depression  when  plotted  in  wall  coordinates.  The  amount  of 
the  depression,  Au+,  can  be  used  to  determine  the  equivalent  sand 
grain  roughness  of  the  surface,  as  described  by  Schlichting.  Applying 
this  method  to  the  non-blown  velocity  profiles  from  these  experiments, 
the  ratio  of  the  equivalent  sand  grain  roughness  to  the  ball  diameter  is 
found  to  be  very  nearly  the  same  for  all  the  velocity  profiles  lying 
between  0.60  and  0.63.  This  confirms  the  value  recommended  by  Schlicht- 
ing for  a similar  surface,  but  tested  in  duct  flow  instead  of  boundary 
layer  type  flow. 

One  difficulty  in  interpreting  rough  surface  velocity  profiles 
arises  in  determining  the  profile  origin.  Typically,  the  apparent  sur- 
face of  the  rough  plate  lies  below  the  tops  of  the  roughness  elements. 
Different  investigators  have  handled  this  problem  in  different  ways. 

Tsuji  and  Iida  [19]  measured  velocity  profiles  from  the  crests  of  the 
roughness  elements.  Others,  such  as  Moore  [12],  Perry  [14,15]  and  Liu 
[16],  place  the  profile  origin  below  the  crests  of  the  roughness  ele- 
ments. Perry  haw  argued  that  the  distance  below  the  rough  element 
crests  should  be  adjusted  until  the  velocity  profile  plotted  in  semi-log 
coordinates  exhibits  the  familiar  'log'  region.  This  approach  estab- 
lishes both  the  all  shear  and  the  'logarithmic  asymptote',  which  is  the 
distance  below  the  roughness  element  crests  to  the  profile  origin.  All 
of  the  methods  which  have  been  suggested  require  some  judgment  by  the 
Individual  examining  the  data.  In  the  present  work,  still  another  method 
was  used:  origins  for  the  velocity  profiles  taken  in  these  experiments 
were  determined  by  extrapolating  the  data  to  zero  In  cartesian  coordi- 
nates, by  eye.  Typically,  the  origins  for  the  turbulent  velocity  pro- 
files fell  about  .010  inches  below  the  crests  of  the  balls.  Laminar  and 
transitional  profiles  were  depressed  less.  Using  the  origins  determined 
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In  this  manner,  the  profiles  were  integrated  to  obtain  momentum  thick- 
nesses which  were  then  used  to  evaluate  akin  friction.  This  approach 
provies  a convenient  basis  for  determining  origins  for  all  the  profiles 
and  yield  results  similar  to  what  would  have  been  obtained  using  Perry's 
method  for  the  fully  turbulent  profiles. 

One  interesting  result  obtained  from  the  velocity  profile  data  is 
that  the  apparent  roughness  of  the  plate  seems  independent  of  both  dis- 
tance along  the  plate  (boundary  layer  thickness)  and  of  free-stream 
velocity  (viscosity  effects).  Thus  the  roughness  effects  seem  to  scale 
on  roughness  particle  size,  not  on  'relative  particle  size'.  This  is 
different  from  rough  pipe  data,  where  roughness  effects  scale  on  (k^/d), 
the  ratio  of  the  roughness  size  to  the  pipe  diameter. 

Early  rough  surface  boundary  layer  experiments  demonstrated  that  the 

skin  friction  becomes  a function  of  x/k  alone.  It  follows  from  this 

s 

that  the  momentum  thickness  also  must  be  a function  of  x/k  alone  for 

tt 

the  zero -pressure  gradient  flow  if  the  surface  roughness  is  not  changing. 

Fig.  3.5  shows  measured  values  of  0 versus  (x-x  ) for  each  of  the 

o 

test  sectiou  velocities,  being  the  appropriate  virtual  origin.  The 

growth  curves  of  all  the  boundary  layers  are  substantially  the  same. 

Also  si  town  is  the  expected  behavior  of  a smooth  plate  at  the  same  test- 
section  velocities.  Although  the  32  fps  data  are  somewhat  low  and  the 
190  and  242  fps  data  show  a slight  shift  at  the  downstream  end,  the  co- 
herence is  good  compared  to  the  smooth  plate  curves.  The  two-dimensional 
Integral  boundary  layer  equation  for  flat  plate  unblown  flows  is: 


From  Fig.  3.5  it  can  be  seen  that  d9/dx  is  substantially  indepen- 
dent of  velocity,  above  32  fps,  being  the  same  function  of  0 for  all 
velocities.  From  this  observation  c^/2  is  mainly  a function  of  0 
alone.  This  might  have  been  anticipated  from  early  data  which  shows! 
cf/2  was  only  a function  of  x/k8  and  independent  of  velocity. 

Skin  friction,  deduced  by  differentiating  the  individual  curves  of 
0 versus  x,  is  shown  in  Fig.  3.6  plotted  against  the  ratio  of  momentum 
thickness  to  ball  radius.  The  32  fps  data  are  clearly  low,  and  there 
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remains  a small  velocity  dependence  In  the  other  data.  At  the  tvo  high- 
est velocities,  the  skin  friction  lies  102  to  152  above  the  other  data 
at  the  downstream  end  of  the  plate.  It's  difficult  to  assess  the  sig- 
nificance of  this,  based  on  the  present  data.  An  error  in  measurement 
of  0 of  only  32  to  52  could  produce  the  difference?  observed  in  the 
friction  factors.  Further,  smooth  surface  data,  if  plotted  in  these  same 
coordinates  would  place  the  low-velocity  skin  friction  above  high- 
velocity  values.  This  suggests  that  the  velocity  dependence  observed 
here  is  probably  not  an  approach  to  the  'fully  rough  state'.  Finally, 
the  Stanton  number  data  show  no  such  tendence,  and  they  are  less  sensi- 
tive to  errors  which  may  be  present  in  the  skin  friction  data. 

Figure  3.0  also  shows  the  Stanton  number  data  for  all  five  velocities 
used,  plotted  versus  A/r,  where  A is  the  enthalpy  thickness  of  the 
boundary  layer.  The  data  are  coherent  for  all  velocities  in  the  turbu- 
lent region,  including  the  data  for  32  fps.  After  an  initial  overshoot 
of  a few  percent,  the  32  fps  Stanton  number  data  fall  on  top  of  the  data 
for  all  higher  velocities.  The  turbulent  data  show  no  velocity  depen- 
dence; the  Stanton  number  is  a function  only  of  A/r,  for  all  veloci- 
ties. It  is  Interesting  to  compare  the  values  of  Stanton  number  and 
friction  factor  at  the  same  boundary  layer  thicknesses  (i.e.,  A/r  » 

0/r).  Friction  factors  for  90  and  140  fps  tend  to  lie  about  52  above 
the  Stanton  number  values  ail  along  the  plate,  rather  than  152  low,  as 
in  the  smooth  plate  case.  The  agreement  between  the  32  fps  data  and 
those  for  higher  velocities  was  surprising,  in  view  of  the  friction  fac- 
tor results.  Furthermore,  the  roughness  Reynolds  number  for  that  veloc- 
ity lies  between  22  and  30,  hardly  'fully  rough'. 

To  further  investigate  the  Stanton  number  behavior  at  low  roughness 
Reynrxds  numbers,  a special  test  was  conducted  at  a free-stream  velocity 
of  18.8  fps.  The  boundary  layer  was  augmented  by  blowing  at  F ■ .004 
over  the  first  two  feet  of  the  test  section.  The  boundary  layer  was  then 
allowed  to  relax  into  its  natural  state.  Stanton  number  data  for  the  un- 
blown plates  are  shown  in  Fig.  3.7,  along  with  the  mean  of  the  90  fps 
data.  It  is  notable  that  the  low  velocity  Stanton  number  data  form  a 
natural  extension  of  the  unblown  90  fps  data.  The  roughness  Reynolds 
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number  at  the  end  of  the  test  section  for  the  18.8  fps  run  was  only  14, 
as  determined  from  hot  wire  measurements  of  the  shear  stress  near  the 
surface.  This  is  far  below  the  'fully  rough'  condition  and  nearly 
'smooth'  by  the  usual  criteria.  It  seems  clear  that  'fully  rough'  beha- 
vior of  the  heat  transfer  data  persists  to  much  lower  roughness  Reynolds 
numbers,  for  the  present  surface  geometry,  than  had  been  expected.  The 
classification  of  'rough'  or  'smooth'  based  on  roughness  Reynolds  number 
does  not  seem  to  be  a reliable  indicator  for  the  heat  transfer  perfor- 
mance of  rough  surface  tested  here. 

In  the  turbulent  range,  the  present  data  fit  the  following  correla- 
tion: 

...-0.25 

Rough  St  - .0043 . (3-6) 

The  usual  smooth-plate  correlation  in  terms  of  enthalpy  thickness  Rey- 
nolds number  can  be  cast  in  a similar  form: 


Smooth 


St 


(3-7) 


In  this  form  the  Stanton  number  for  both  rough  and  smooth  behavior 
is  substantially  the  same  at  32  fps.  At  18.8  fps,  the  smooth  plate  cor- 
relation predicts  a value  nearly  502  higher  than  the  rough.  This  is 
because  of  the  inverse  dependence  of  the  smooth  correlation  on  free- 
stream  velocity.  Clearly,  the  smooth  plate  correlation  in  this  form 
does  not  provide  a suitable  limit  for  the  rough  surface  behavior. 

In  summary,  the  unblown  data  in  Figs.  3.5,  3.6  and  3.7  seem  nearly 
independent  of  free-stream  velocity.  Stanton  number  appears  entirely 
independent  of  velocity,  being  only  a function  of  enthalpy  thickness. 
Skin  friction  may  be  Independent  and  at  most  has  a small  dependence. 

Skin  friction,  when  Q/r  is  fixed,  Is  only  about  52  higher  than  Stanton 
number  at  the  same  value  of  A/r.  The  effective  value  of  the  Reynolds 
analogy  factor  is  about  0.95,  Instead  of  1.15,  which  is  typical  of  a 
smooth  plate.  The  boundary  layer  seems  to  be  completely  turbulent,  with 
no  discernible  molecular  effect.  The  effective  value  of  the  turbulent 
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Prandtl  number  seems  near  unity.  Roughness  effects  on  heat  transfer 
persists  to  very  low  values  of  the  roughness  Reynolds  number,  with  no 
evidence  of  the  ’transitional  roughness  range'  in  the  Stanton  number 
data.  Momentum  thickness  and  friction  factor  data  do,  however,  seem  to 
show  a transitional  effect. 

B.  The  Effects  of  Blowing 

Blowing  diminishes  the  Stanton  number  and  skin  friction  for  a rough 
plate,  just  as  it  does  for  a smooth  plate.  Roughness  Rig  skin  friction 
and  Stanton  number  data  at  each  of  the  five  teat  section  velocities  are 
shown  in  Fig.  3.8  through  3.17.  The'>e  data  exhibit  very  similar  trends 
to  the  smooth  plate  skin  friction  a ad  heat  transfer  data  presented  by 
Moffat  [1]  and  Simpson  [2]. 

An  important  feature  of  the  unblown  turbulent  boundary  layer  Stanton- 
number  data  is  its  lack  of  dependence  on  free-stream  velocity.  This  in- 
dependence is  preserved  with  blowing.  Fig.  3.18  shows  Stanton  number 
data  for  all  the  blown  boundary  layers  at  all  velocities  plotted  versus 
A/r  parametric  in  F.  Although  the  scatter  is  increased,  the  fully 
rough  blown  turbulent  boundary  layer  exhibits  the  same  independence  of 
velocity  as  does  the  unblown  boundary  layer. 

In  view  of  this  agreement,  it  is  of  interest  to  examine  the  roughness 
Reynolds  number  range  of  the  data,  since  this  is  so  often  used  as  the 
measure  of  the  rough  surface  behavior.  At  32  fps  and  F - 0.002,  the 
roughness  Reynolds  number  is  about  17  at  the  point  farthest  downstream. 

For  all  higher  velocities  or  lover  blowing  fractions,  ReT  is  greater. 
Higher  blowing  reduces  the  wall  shear  still  further,  and,  for  32  fps, 
takes  the  roughness  Reynolds  number  down  to  10  (F  - 0.004  at  the  20th 
plate).  No  reliable  friction  data  are  available  for  F * 0.008  at  32 
fps i but  Stanton  number  data  for  F - 0.008  are  shown  in  Fig.  3.18.  The 
32  fpi?  Stanton  number  data  at  both  F * 0.004  and  0.008  are  higher 
than  the  remaining  data,  and  it  is  felt  that  they  may  represent  a differ- 
ent state  of  the  boundary  layer  than  the  higher  velocities.  The  scatter 
in  these  data  sometimes  (though  not  often)  exceeds  ± 0.0001  Stanton- 
number  units.  No  cause  has  been  assigned. 
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Values  of  Stanton  number  with  blowing  are  predictable  from  the  un- 
blown values  bya  Couette  relationship  used  for  smooth  plate  flows,  with 
proper  interpretation: 


[M±SU]l'25a  + B)-25 


(3-8) 


The  fully  rough  prediction  is  to  be  made  at  constant  A , whereas  the 
smooth  plate  prediction  is  made  at  constant  Re^ . Predictions  based  on 
this  equation  are  shown  as  solid  lines  in  Figs.  3.19  through  3.23.  The 
agreement  is  excellent  for  all  data  more  than  about  30  boundary  layer 
thicknesses  past  the  peak  of  the  transition  hump.  Data  points  more  than 
30  boundary  layer  thicknesses  past  transition  are  accented  in  Figs.  3.19 
through  3.23.  The  application  of  the  Couette  flow  estimator  baaed  on 
enthalpy  thickness  instead  of  enthalpy  thickness  Reynolds  number  is  very 
reasonable  in  view  of  the  success  of  this  coordinate  in  organizing  the 
Stanton  number  data. 

The  skin  friction  data  with  blowing  show  a high  degree  of  scatter 
and  like  the  unblown  data  it  is  difficult  to  draw  any  definitive  conclu- 
sions from  them.  It  is  Interesting,  however,  to  examine  the  momentum 
thicknesses  of  the  blown  boundary  layer.  These  measurements  are  less 
prone  to  uncertainties  and  provide  a better  basis  for  interpreting  the 
rough  surface  behavior.  Fig.  3.24  shows  all  of  the  momentum  thicknesses, 
plotted  against  distance  from  the  virtual  origin,  parametric  in  F.  It 
is  clear  that  in  these  coordinates,  the  boundary  layers  for  all  veloci- 
ties are  substantially  the  same.  This  leads  one  to  suspect  that  in  terms 
of  local  thickness  coordinates  the  skin  friction  should  be  the  same. 

There  is  so  much  scatter  in  the  blown  skin  friction  data  that  this  cannot 
be  cunfirxttd.  It  does  seem  likely,  Siowever,  that  the  skin  friction  has 
only  a small  dependence  on  free-stream  velocity,  since  that  was  observed 
in  the  unblown  data.  More  accurate  skin  friction  measurements  will  be 
taken  in  the  continuing  experimental  program  on  the  rough  surface  bound- 
ary layer,  and  this  question  will  undoubtedly  be  answered  in  the  future. 

The  curve  fitting  technique  used  determine  skin  friction  did  not 
work  well  at  blowing  fractions  greater  ,han  .004.  The  skin  friction  is 
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so  small  at  these  high  bloving  rates  that  even  small  uncertainties  in 
the  momentum  thickness  measurement  result  In  differences  In  skin  friction 
that  are  percentage-vise  large.  Because  of  this,  skin  friction  predic- 
tions for  the  .008  bloving  fraction  runs  have  not  been  included  in  the 
figures.  The  skin  friction  at  90  fpo  and  F ■ .002  vas  accorded  some 
special  treatment.  The  curve  fit  technique  predicted  a virtual  origin 
for  the  boundary  layer  that  was  inconsistent  with  the  other  data  at  that 
velocity.  In  this  one  particular  data  set,  the  virtual  origin  vas  selec- 
ted by  extrapolating  the  boundary  layer  momentum  thicknesses  to  zero  in 
cartesian  coordinates.  The  skin  friction  vas  then  evaluated  with  the 
same  curve  fit  routine  but  with  fixed  virtual  origin.  In  all  other  data 
sets,  the  curve  fitting  routine  was  allowed  to  select  both  the  virtual 
origin  snd  the  curve  from  which  the  skin  friction  data  were  determined. 

The  effects  of  blowing  on  heat  transfer  and  skin  friction  to  a 
'fully  rough*  turbulent  boundary  layer  can  be  summarized  as  follows. 
Stanton  number  is  a function  of  A and  F only,  and  the  blown  value 
can  be  predicted  from  the  customary  predictive  equation,  but  evaluated  at 
constant  A instead  of  constant  Re^.  The  effects  of  blowing  on  skin 
friction  are  less  well  defined,  but  boundary  layer  thickness  measure- 
ments indicate  that  the  growth  of  the  momentum  boundary  layer  is  sub- 
stantially independent  of  velocity  and  a function  only  of  the  blowing 
fraction.  Baaed  on  this  and  the  behavior  of  the  unblown  skin  friction 
data,  it  seems  likely  that  the  rough  surface  blown  skin  friction  has  at 
most  small  dependence  of  free-sU'cam  velocity.  In  view  of  the  data  scat- 
ter, more  definitive  conclusions  are  not  possible. 

C.  Roughness  and  Blowing  Effects  on  Transition 

All  boundary  layers  in  these  experiments  were  allowed  natural  tran- 
sitions wi‘Sout  a boundary  layer  trip.  There  are  two  advantages  of 
'tripping'  the  boundary  layer  and  having  the  virtual  origin  at  the  test- 
section  inlet:  the  momentum  baundary  layer  origin  is  matched  to  the  ther- 
mal boundary  layer,  and  it  makes  maximum  use  of  the  test  section  length. 
However,  when  the  well-defined  laminar  boundary  layer  and  transition 
region  were  found  at  the  32  fps  test  section  velocity,  it  was  decided 
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that  the  opportunity  to  observe  the  effects  of  roughness  and  blowing  on 
transition  more  than  offset  any  advantages  offered  by  tripping  the  bound- 
ary layer.  There  is  very  little  experimental  information  in  the  litera- 
ture about  transition  of  the  blown  boundary  layer  on  a rough  wall. 

The  present  experimental  program  was  not  specifically  designed  to  study 
transition,  but  some  interesting  observations  could  he  made  from  the 

data  concerning  rough  surface  boundary  layer  transition. 

The  effect  of  the  laminar  boundary  layer  on  Stanton  number  can 
clearly  be  seen  in  the  non-blowing  data  in  Fig.  3.1.  Increasing  the 
free-stream  velocity  moves  the  transition  upstream.  It  is  interesting 
to  note,  however,  from  Fig.  3.2  that  the  transition  is  well  correlated 
by  momentum  thickness  Reynolds  number.  Transition  occurs  at  a momentum 
thickness  (and  enthalpy  thickness)  Reynolds  number  of  about  400. 

The  effects  of  blowing  on  transition  can  be  seen  using  the  data  for 
32  fps  in  Fig.  3.10.  Shown  on  this  figure  are  momentum  thickness  Rey- 
nolds numbers  at  the  onset  of  transition.  The  data  clearly  show  that 
blowing  moves  the  transition  upstream  but  does  not  much  affect  the 
momentum  thickness  Reynolds  number  (about  400)  at  which  the  transition 
occurs.  It  is  important  to  note  that  the  thermal  protection  offered  by 
blowing  does  not  offset  the  effect  of  an  early  transition.  Blowing 
through  a region  which  would  otherwise  have  remained  laminar  can  result 
in  a much  higher  heat  load  than  no  blowing  at  all. 

The  present  data  indicate  that,  like  pipe  flow  data,  below  a certain 
criticai  Reynolds  number  the  rough  surface  boundary  layer  remains  laminar. 
For  these  experiments,  transition  occurs  at  or  near  the  same  Reynolds 
number  that  would  be  expected  for  a smooth  surface  transition.  There 
are  signs  of  transition  in  all  of  Sranton  number  data,  even  at  the  higher 
test  section  velocities. 

A final  remark  that  can  be  made  about  the  laminar-to-turbulent  tran- 
sition observed  in  these  tests  is  that  once  initiated,  it  occurs  in  a 
relatively  short  distance.  Smooth  plate  boundary  layer  transition  is 
often  spread  over  a larger  Reynolds  number  range.  Fig.  3.25  compares 
smooth  plate  transition  data  by  Reynolds  [53]  to  the  present  experiment. 
The  smooth  plate  Stanton  number  after  transition  approaches  the  asymp- 
totic turbulent  value  slower  than  the  rough  surface  data  from  these 
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experiments.  A possible  explanation  of  this  is  that  the  surface 
roughness  elements  may  aid  the  transition  process  once  initiated  and 
help  bring  it  to  completion  over  a shorter  section  of  the  plate. 

D.  Comparison  with  Other  Rough  Surface  Experiments 

Direct  comparison  of  skin  friction  and  heat  transfer  data  to  similar 
rough  surface  experiments  is  useful  because  it  demonstrates  that  the  pres- 
sent  results  are  similar  to  what  has  been  found  in  other  rough  surface 
experiments.  Several  investigators  have  compared  their  rough  surface 
skin  friction  to  the  'fully  rough'  skin  friction  correlation  suggested  by 
Frandtl  and  Schlichting  [10]. 

c,  - (2.87  + 1.58  log  x/k  )“2,5  . (3-9) 

I 8 

Moore's  [12]  data  compare  reasonably  well  with  this  equation,  but  he 
found  different  size  - 'ughness  elements  were  not  organized  in  these  coor- 
dinates. The  data  by  Liu  [16]  were  reasonably  well  organized,  but  the 
predicted  skin  friction  was  about  202  higher  than  the  data.  Wu  [18] 
measured  average  akin  friction  over  a 10  * 21  inch  test  plate  which  agrees 
well  with  the  Prandtl-Schlichtiag  equation  for  the  average  skin  friction. 
Skin  friction  measurements  from  these  experiments  are  shown  in  Fig.  3.26. 
Like  Liu’s,  the  present  data  lie  about  202  lower  than  the  Prandtl- 
Schlichting  curve. 

In  a more  recent  and  extensive  study,  Lakshoan  and  Jayatilleke  [52] 
examine  nearly  every  available  set  of  rough  surface  skin  friction  data 
from  both  pipes  and  plates.  They  correlate  the  data  in  E and  Re^ 
coordinates,  where  E Is  defined  by  the  equation: 

u+  - ~ 2.n(Ey+)  . (3-10) 

Data  from  these  experiments  are  shown  in  Fig.  3.27  in  F.  versus  Re, 
coordinates.  Also  shown  is  the  line  for  sand-grain  roughness.  Although 
the  data  lie  above  the  sand-grain  roughness  line  in  the  transition  rough- 
ness Reynolds  number  range,  it  is  reasonably  well  organized  in  these 
coordinates  and  agree  well  at  ’fully  rough’  roughness  Reynolds  numbers. 
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A direct  comparison  of  the  heat  transfer  results  with  other  rough 
surface  experiments  is  more  difficult  because  previous  rough  surface 
heat  transfer  experiments  have  been  performed  in  pipe  geometries.  The 
data  correlations  of  these  tests  typically  include  pipe  dimensions  or 
have  their  rougliness  scaled  by  pipe  diameters.  The  heat  transfer  data 
by  Dipprcy  and  Sabersky  [24]  were  correlated  by  the  following  expression: 


g(ReT,Pr)  - A(Re^) 


(3-11) 


For  luily  rough  flow,  the  q and  A functions  are  simplified  to  give: 


+ 8.4d  - 5.19  Re°*2  Pr°*44  - gfr<ReT,Pr)  . (3-12) 

Dipprcy  and  Sabersky  were  able  to  correlate  all  of  their  fully  rough  data 
in  and  Kct  coordinates.  Owen  and  Thompson  [25]  recommend  the 

following  expression  for  heat  transfer  from  a flat  plate: 


K/2 
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where  8*  is  the  sublayer  Stanton  number.  Based  on  their  experiments  and 
data  by  others,  the  following  correlation  was  recommended  for  the  sub- 
layer St an Lou  number: 


B* 


Pr 
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(3-14) 


with  a lying  between  0.45  and  0.7  but  represented  best  by  0.52.  The 
Owen-Thumpson  expression  can  be  cast  into  a form  similar  to  the  Dipprey- 
Sabersky  correlation: 

1 

VV1 

These  two  models  for  rough  surface  heat  transfer  are  plotted  In 
Fig.  3.28  along  with  data  from  the  present  experiments.  Both  correla- 
tions agree  reasonably  well  for  roughness  Reynolds  numbers  greater  than 
70.  The  data  from  the  present  experiments  fall  below  these  predictions. 


- 1 | + 8.48  - 0.52  Re°*45  Pr0,8  + 8.48 


(3-15) 
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Both  of  these  correlations  employ  a skin  friction-to-Stanton  number 
ratio,  minus  one.  The  correlation  in  this  form  emphasizes  experimental 
uncertainties  of  the  present  data  by  looking  at  differences  between 
nearly  equal  experimentally  determined  values. 

Another  useful  data  comparison  can  be  made  using  the  correlation 
suggested  by  Nunner  [23].  This  correlation  compares  the  increase  in 
rough  surface  skin  friction  to  the  increase  in  heat  transfer.  Data  from 
these  experiments  is  shown  in  Figure  3.29.  Also  shown  is  Nunner's  cor- 
relation and  the  more  recent  correlation  recommended  by  Norris  [30].  The 
data  from  these  experiments  is  best  organized  by  Nunner's  original  corre- 
lation. 

Based  on  these  comparisons,  it  appears  that  the  rough  surface  skin 
friction  and  heat  transfer  data  from  the  present  experiments  are  compar- 
able with  other  similar  experiments. 


66 


STREAM  VELOCITY 
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Fig.  3.2  Rough  surface  skin  friction  and  Stanton  number  versus  momentum  and  enthalpy 
thickness  Reynolds  number 


Fig.  3.3  Velocity  profiles  at  different  test  section  velocities 
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surface  skin  friction  versus  (momentum  thickness)  / (ball  radius)  and 
surface  Stanton  number  versus  (enthalp:’  thickness) / (ball  radius) 
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Fig,  3.8  Rough  surface  ekln  friction  versus 
x- Reynaldo  number  at  u ■ 32  fps 
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Fig.  3.S  Rough  surface  Stanton  number  versus  x~  Reynolds  number  at 

3'  * 


r*T*] 


u« 

- 139 

ft/sec 

o 

r $ 

x X 

o 

P 

o 

p 

o 

X 

o 

<•> 

X 

B 1,0'. 'INC  FRACTION’ 

&> 

Fig.  3.12  Rough  surface  skin  friction  versus 
x-Reynolds  number  at  u - 139  fps 
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Fig.  3.13  Rough  surface  Stanton  number  versus  x- Reynolds  number  at 
u * 139  fps 


Fig.  3.18  Stanton  number  versus  (enthalpy  thickness)/ (ball  radius)  for  all  boundary  layers 


13S  ft/sec 


Fig.  3.21  Stanton  number  versus  (enthalpy  thickness)/ (ball  radius)  at  u - 139  fps 
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Fig.  3.24  Momentum  thickness  versus  distance  from  virtual  origin  for  all  boundary  layers 
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(Note:  lines  shown  as  visual  aid  only) 


Fig.  3.25  Comparison  of  rough  surface  transition  Stanton  number  to 
smooth  plate  transition  data  by  Reynolds  [53] 
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Fig.  3.26  Comparison  of  skin  friction  data  with  P rand tl~ Schlichting 
correlation 
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Fig.  3.27  Comparison  of  skin  frictiou  data  in  E versus  Re^  coordinates 
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Fig.  3.28  Comparison  of  heat  transfer  data  tq  correlations  by  Dipprey 
and  Sabersky  [24]  and  Own  and  Thompson  [25] 
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CHAPTER  IV 


PREDICTION  OF  THE  ROUGH  SURFACE  DATA 

Prediction*  of  the  rough  surface  boundary  layers  havs  been  made 
using  a finite-difference  prediction  scheme  based  on  the  computer  program 
originated  by  Patankar  and  Spalding  [37 J.  The  turbulence  model  used  in 
this  program  lias  undergone  substantial  development  at  Stanford  in  order 
to  equip  it  to  predict  smooth  surface,  transpired  turbulent  boundary 
layers  with  pressure  gradients.  A brief  description  of  the  program  is 
given  by  Kays  [38],  who  also  gives  examples  of  its  use. 

A.  Hydrodynamic  Predictions  of  the  Rough  Surface  Boundary  Layer 

Closure  of  the  boundary  layer  equations  is  provided  by  a mixing- 
length  scheme  using  a modified  van  Driest  model  in  the  Inner  regiou. 

This  model  lias  been  extensively  and  successfully  used  by  several  similar 
prediction  programs.  In  wall  coordinates.  It  is  given  by: 

1+  - k/(i  - e-y+/A+)  . (4-1) 

Efl 

Near  a smooth  wall,  the  mixing~length  is  'damped'  to  allow  viscous 
effects  to  become  impatient,  since  molecular  viscosity  dominates  the 
flow  in  that  region.  Farther  away  from  the  wall,  the  mixing-length  be- 
comes a linear  function  of  distance  from  the  wall,  following  Prandtl's 
assumption  for  mixing-length  behavior  in  the  logarithmic  region. 

The  van  Driest  parameter  A+  may  be  Interpreted  as  a measure  of  the 
sublayer  thickness.  Based  on  this  interpretation,  several  authors  have 
represented  as  a function  of  pressure  gradient  and  blowing.  These 

parameters  are  known  to  affect  the  thickness  of  the  wall  region  of  the 
boundary  layer.  Andersen  [8]  gives  a brief  discussion  of  the  role  of 
A*  and  reviews  some  interpretations  which  have  been  given  to  it. 

In  the  outer  region  of  the  boundary  layer,  the  mixing  length  is  held 
at  a constant  fraction  of  the  992  boundary  layer  thickness. 
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Values  of  X - 0.08  to  0.09  have  been  demonstrated  to  be  satis- 
factory for  smooth  wall  boundary  layers  by  several  authors  [35,38,50] . 

In  these  predictions,  X - 0.08  was  used. 

It  was  suggested  by  van  Driest  [39]  that  roughness  effects  could  be 
simulated  by  changing  the  damping  expression,  since  the  presence  of 
roughness  elements  In  the  flow  seemed  to  reduce  the  sublayer  thickness. 

This  can  be  simulated  in  the  van  Driest  mixing-length  model  by  reducing 

+ 4* 

A . In  the  present  program,  the  relationship  between  A and  roughness 

was  investigated  by  systematically  predicting  pipe  flows  over  a series  of 
Reynolds  numbers.  It  was  felt  that  the  relative  simplicity  of  the  pipe 
geometry  would  provide  results  which  could  then  be  generalized  to  the 
boundary  layer  cose.  Pipe  flow  predictions  using  different  dancing  con- 
stants produced  results  which  verified  that  roughness  effects  could 
indeed  be  simulated  by  making  A*  a function  of  the  roughness  Reynolds 
number.  Fig.  4.1  shows  the  functional  forms  of  A*  which  were  deter- 
mined for  sand-grain  roughness,  commercial  roughness,  and  the  rough  sur- 
face used  in  the  present  experiments.  The  A+  correlation  for  sand- 
grain  roughness  was  obtained  using  Nlkuradse'a  expression  [9]  relating 
friction  factor  to  the  roughness  Reynolds  number.  The  commercial  rough- 
ncsu  A correlation  was  based  on  relative  roughness  read  directly  from 
a Moody  chart  at  friction  factors  calculated  for  various  pipe  Reynolds 
numbers.  The  A*  correlation  for  the  rough  surface  used  in  these  exper- 
iments was  based  on  the  friction  factor  measurements  made  at  the  32 
ft/sec  tunnel  velocity.  This  was  the  only  tunnel  velocity  at  which  the 
roughness  Reynolds  number,  based  on  an  equivalent  sand-grain  roughness, 
was  below  55.  It  is  not  surprising  that  the  A+  behavior  of  the  uniform- 
ball  roughness  used  in  these  experiments  is  different  than  either  commer- 
cial or  sand-grain  roughness.  Roughness  elements  in  a commercially  rough 
surface  are  irregular  and  may  Include  some  very  large  elements  which  c?»u 
become  'active*  and  begin  to  produce  roughness  effects  at  relatively  low 
values  of  the  nominal  roughness  Reynolds  number.  The  roughness  elements 
in  the  sand-grain  roughened  surface  are  more  uniform,  and  the  onset  of 
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’transitional  roughness'  is  delayed  to  higher  roughness  Reynolds  numbers. 
The  extremely  uniform  nature  of  the  roughness  used  in  these  experiments 
still  further  delayed  the  onset  of  transitional  roughness.  This  is  the 
behavior  shown  in  the  A correlations  for  the  different  surfaces.  As 
would  be  expected,  the  different  A+  correlations  come  together  as  the 
roughness  Reynolds  number  approaches  the  'fully  rough'  value  and  at  a 
roughuess  Reynolds  number  of  55,  A+  has  reached  zero  for  all  surfaces. 
This  is  consistent  with  the  idea  that  at  the  onset  of  'fully  rough'  be- 
havior, the  roughness  elements  protrude  outside  the  sublayer.  The  sub- 
layer thickness  has  effectively  gone  to  zero  when  A*  has  gone  to  zero. 
Without  a sublayer,  form  drag  replaces  viscous  shear  and  the  flow  resis- 
tance becomes  independent  of  viscosity  — an  important  characteristic  of 
'fully  rough'  flow. 

In  the  prediction  program,  the  turbulent  shear  transmitted  to  the 
wall  is  calculated  from  an  eddy  diffusivity  which  is  based  on  the  mixing- 
length  model.  One  way  to  represent  the  'fully  rough*  behavior  at  the 
wall  is  to  make  the  mixing-length  and,  therefore,  the  eddy  diffusivity 
non-zero  at  the  wall.  There  is  evidence  that  the  rough  wall  affects 
mainly  the  How  very  near  the  wall.  Clauser  [48]  has  shown  that  velocity 
profiles  taken  over  both  smooth  and  rough  surfaces  are  similar  when  plot- 
ted in  defect  coordinates.  His  conclusion  was  that  the  outer  flow  (which 
is  some  90Z  of  the  boundary  layer)  is  not  affected  by  the  surface  rough- 
ness. In  these  coordinates  both  smooth  and  rough  wall  velocity  profiles 
are  similar.  The  prediction  scheme  of  the  present  work  is  based  on  these 
observations. 

The  mixing- length  was  made  non-zero  at  the  wall,  becoming  asymptotic 
to  the  smooth  wall  behavior  away  from  the  wall  by  expressing  the  mixing- 
length  in  the  following  manner: 

- v/(icyV  + (A**)2  . (4-3) 

m o 

By  again  predicting  pipe  flows,  a correlation  of  the  value  of  the  mixing 
length  at  the  wall,  Ai.*,  versus  roughness  Reynolds  number  was  obtained. 
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(4-4) 


Figure  4.2  shows  the  variation  of  A £ with  roughness  Reynolds  number 

along  with  the  data  predictions  used  to  obtain  Eqn.  (4-4).  It  is  inter- 
esting to  note  that  for  large  values  of  roughness  Reynolds  number  (Re^  > 
100)  the  functional  relationship  reduces  to  a simple  linear  variation 
between  mixing-length  offset,  Ai,Q,  and  equivalent  sand-grain  roughness. 
The  effect  of  roughness  on  mixing-length  is  related  only  to  the  size  of 
the  roughness  elements,  an  important,  characteristic  of  'fully  rough' 
behavior.  With  this  correlation,  it  is  possible  to  predict  the  behavior 
of  all  of  the  surface  roughnesses  which  require  different  A+  correla- 
tions. This  was  expected,  since  in  fully  developed  roughness  the  beha- 
vior of  many  different  surfaces  can  be  expressed  in  terms  of  an  equiva- 
lent sand-grain  roughness. 

The  mixing-length  variation  as  a function  of  roughness  Reynolds 
number  is  shown  in  Fig.  4.3.  Fig.  4.4  shows  predictions  of  Nikuradse'c 
[9]  rough  pipe  data.  These  predictions  were  made  by  modeling  the  Nlku- 
radse  experiment  with  the  prediction  program  using  the  mixing-length 
model  Just  described.  The  sand-grain  roughness  A+  correlation  was  used 
for  roughness  Reynolds  numbers  below  55  and  the  correlation  for  mixing- 
length  value  at  the  wall,  A&+,  above  55.  The  rough  pipe  behavior  was 
well  predicted  over  a wide  range  of  roughnesses  and  pipe  Reynolds  num- 
bers. 


Predictions  were  made  of  the  rough  surface  boundary  layer  from  the 
present  experiments  using  these  same  correlations.  For  roughness  Rey- 
nolds numbers  below  55,  it  was  necessary  to  use  an  A+  correlation 

selected  for  this  particular  type  of  roughness.  For  roughness  Reynolds 
numbers  greater  than  55,  the  same  mixing-length  offset  correlation  al- 
ready described  was  used.  Figs.  4.5  and  4.6  show  a comparison  of  mea- 
sured and  predicted  velocity  profiles  at  velocities  of  32  fps  and  139  fps. 

In  general,  measured  and  predicted  profiles  in  wall  coordinates  compare 
favorably  as  long  as  the  measured  and  predicted  skin  frictions  are  com- 
parable. This  is  generally  the  case  over  the  last  two- thirds  of  the  test 
section.  Fig.  4.7  compares  the  measured  and  predicted  skin  friction  for 
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three  test  section  velocities.  In  r early  every  case  the  predicted  values 
are  greater  than  measured  near  the  front  of  the  test  section.  This  same 
trend  has  also  been  observed  in  the  smooth  wall  predictions  and  probabJy 
results  because  initially  the  boundary  layer  is  still  recovering  from 
transition  and  attempting  to  reach  an  equilibrium  state. 

In  making  these  predictions,  it  was  found  that  blowing  decreased 
the  predicted  skin  friction.  In  order  to  keep  this  reduced  skin  fric- 
tion from  reducing  the  effect  of  the  roughness,  a first-order  blowing 
correction  has  been  added  to  the  correlation  for  the  mixing-length  at 
the  wall,  A£*.  With  the  blowing  correction,  it  becomes: 
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One  would  expect  that  blowing  should  play  the  same  role  in  the 
mixing-length  formulation  no  matter  if  the  wall  is  rough  or  smooth.  The 
expression  for  roughness  Reynolds  number  with  its  blowing  correction  can 
be  rewritten  as  the  sum  of  two  roughness  related  Reynolds  numbers: 


• «•:>)  -(¥*)•■  (£) 


(4-6) 


It  is  reasonable  that  if  there  is  blowing  at  the  wall  this  should  affect 
the  mixing-length  correlation  and  that  the  effect  be  similar  to  the  cor- 
rection that  other  investigators  [3,35,49,50]  have  applied  to  A+  to 
ensure  correct  prediction  of  smooth  surface  skin  friction  with  blowing. 
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B . Heat  Transfer  Predictions  for  the  Rough  Surface  Boundary  Layer 

Heat  transfer  predictions  were  made  by  solving  the  boundary  layer 
energy  equation  using  a turbulent  Prandtl  number  model  for  closure.  A 
constant  turbulent  Prandtl  number  is  often  used  for  smooth  surface  bound- 
ary layer  predictions.  There  is  not  complete  agreement  among  experimen- 
ters about  the  behavior  of  the  turbulent  Prandtl  number  with  distance 
from  the  wall,  but  most  agree  that  it  increases  near  the  wall.  One 
turbulenl  Prandtl  number  expression  which  has  been  used  successfully 
for  air  [38]  and  which  reflects  this  trend  is  given  by: 


Pr  . - (1.43  - 0.17  y+i/A) 

t, smooth  1 


for  Pr  , > 0.86.  If  Pr  ..  < 0.86,  then 

t,  smooth  — t, smooth 


(4-7) 


smooth 


0.86  . 


A plot  of  this  model  is  shown  in  Fig.  4.8,  and  it  can  be  noted  that 
near  the  wall  the  Prandtl  number  increases  to  1.43,  which  is  approximately 
the  reciprocal  of  the  molecular  Prandtl  number  for  air.  Also  shown  is 
the  constant  turbulent  Prandtl  number  model  often  us-2d  for  air  as  well 
as  otlser,  higher  Prandtl  number  fluids. 

Figure  4.9  illustrates  the  results  obtained  using  the  smooth  si  rface 
turbulent  Prandtl  number  model  for  rough  pipe  heat  transfer  to  air.  The 
predictions  are  plotted  in  the  coordinates  used  by  Nunner  [23]  to  organ- 
ize his  rough  pipe  heat  transfer  data;  rough- to-smoouh  Stanton  number 
versus  rough-to-smooth  skin  friction.  While  the  predictions  are  reason- 
ably well  organized  in  these  coordinates,  probably  due  to  the  correct 
modeling  of  the  skin  friction,  the  general  trend  of  the  data  is  high.  Use 
of  this  smooth  wall  turbulent  Prandtl  number  over-predicts  the  heat  trans- 
fer. This  was  not  unexpected.  The  smooth  surface  turbulent  Prandtl  num- 
ber model  forces  the  momentum  and  thermal  eddy  dif fusivities , and  therefore 
the  transport  process,  themselves,  to  be  similar.  In  solving  the  momen- 
tum problem,  the  mixing-length  has  been  increased  near  the  wall  to 
account  for  form  drag  and  increased  mixing  at  and  near  the  surface  due 
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to  the  presence  of  roughness  elements.  In  the  heat-transfer  problem, 
the  Increased  mixing  will  increase  the  thermal  transport  across  the 
boundary  layer  near  the  wall,  but  at  the  wall  there  is  no  heat  transfer 
mechanism  analogous  to  the  pressure  force  transmitted  to  the  wall  as 
form  drag.  At  the  wall,  heat  can  be  transferred  only  by  conduction.  To 
account  for  this  difference  using  only  the  parameters  available  to  the 
prediction  program,  the  turbulent  Prandtl  number  must  be  increased  both 
at  and  near  the  wall.  To  accomplish  this,  a thermal  mixing-length  has 
been  defined.  The  rough  surface  turbulent  Prandtl  number  is  then  ex- 
pressed as: 

prt, rough  ‘ prt,Smooth(r)  <4-S> 

The  thermal  mixing-length  Is  defined  with  the  same  modified  van 
Driest  model  used  for  the  momentum  mixing-length  and  uses  the  same  rough- 
ness correlations.  However,  when  evaluating  the  thermal  mixing- Length,  a 
reduced  roughness  Reynolds  number  is  used.  Since  the  roughness  effects 
are  built  into  the  models  only  near  the  wall,  the  thermal  mixing-length 
is  less  than  the  momentum  mixing-length  near  the  wall,  but  both  approach 
smooth  surface  behavior  away  from  the  wall.  This  produces  an  Increased 
turbulent  Prandtl  number  near  the  wall,  but  approaching  the  smooth 

surface  correlation  away  from  the  wall.  It  was  decided  to  first  apply 

# 

this  ujudel  to  rough  pipe  heat  transfer  data  and  then  generalize  these 
results  to  the  boundary  layer  case.  The  Dipprey  and  Sabersky  [24]  pipe 
data  were  chosen  because  of  the  wide  range  of  roughness  and  fluid  Prandtl 
numbers.  It  soon  became  evident  when  predicting  these  data  that  the 
thermal  mixing-length  correlation,  to  be  successful,  must  be  a function 
of  Prandtl  number  as  well  as  roughness  Reynolds  number.  At  high 
roughness  Reynolds  numbers,  the  increased  mixing- length  used  for  the 
hydrodynamic  predictions  almost  entirely  masked  any  molecular  Prandtl 
number  effects  on  heat  transfer.  To  correctly  predict  the  heat  transfer 
results,  the  amount  by  which  the  roughness  Reynolds  number  is  reduced 
when  evaluating  the  thermal  mixing- length  must  be  a function  of  molecular 
Prandtl  number.  Predictions  of  the  Dipprey  and  Sabersky  data  are  shown 
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in  Figs.  4.10  and  4.11.  It's  apparent  that  the  Prandtl  number  effects 
persist  even  in  the  very  rough  data.  In  making  these  predictions,  the 
roughness  Reynolds  numbers  were  reduced  by  the  factor  C^,  which  was 
given  by: 


- 0.5  for  Re^  < 40  , (4-9) 

CH  - 0.5  - . 0454 (1.9+Pr) (log  Re^-1.602)  . 

In  using  this  correlation,  is  never  allowed  to  become  less 

than  0.1.  A plot  of  is  shown  in  Fig.  4.12.  It  should  be  pointed 

out  that  simply  provides  a convenient  means  by  which  to  reduce  the 

effects  of  roughness  on  the  thermal  mixing-length  at  and  near  the  wall 
and  therefore  produce  a higher  turbulent  Prandtl  number  at  and  near  the 
wall. 

To  prHict  the  heat  transfer  data  from  the  present  tests,  it  is 
sufficient  to  simply  reduce  the  roughness  Reynolds  number  by  half.  As 
can  be  seen  from  Fig.  4.12,  over  the  limited  roughness  Reynolds  number 
range  of  these  experiments,  this  is  essentially  the  same  which 

would  be  produced  by  the  correlation  used  for  prediction  of  the  Dipprey 
and  Sabersky  data.  Fig.  4.13  compares  the  Stanton  number  predictions 
for  the  present  experiments  with  the  data.  As  con  be  seen,  the  agree- 
ment of  the  predictions  with  the  data  is  very  good. 

Finally,  it  should  be  pointed  out  that  the  models  used  here  are 
based  primarily  on  the  observation  that  the  roughness  effect  is  confined 
to  a region  very  near  the  wall,  with  the  outer  region  of  the  boundary 
l&yer  being  relatively  unaffected.  In  the  present  heat  transfer  model, 
if  the  reduced  roughness  correction  were  not  used,  the  increased  turbu- 
lent transport  at  the  wall  would  overshadow  any  molecular  effects.  This 
has  been  shown  experimentally  not  to  be  the  case  — the  Prandtl  number 
of  the  fluid  plays  an  Important  role  in  rough  surface  heat  transfer.  As 
a result,  a Prandtl  number  effect  must  be  included  in  the  correlations 
used  for  the  Stanton  number  predictions.  A possible  reason  why  molecular 
effects  persist  in  rough  surface  heat  transfer  lies  in  the  fact  that  con- 
duction is  the  only  mechanism  available  for  the  exchange  at  the  wall. 
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As  the  height  of  the  roughness  elements  increaser  it  is  likely  that  the 
inner  regions  can  be  modeled  as  a nearly  stationary  fluid  film.  Heat 
must  be  transferred  through  this  film  by  conduction,  hence  the  molecular 
Prandtl  number  becomes  important.  This  fluid  film  is  not  part  of  the 
present  finite-difference  prediction  program.  An  improved  modeling 
might  include  a thermal  conduction  layer  in  the  energy  equation  to  better 
represent  the  molecular  effect  that  is  present  in  the  rough  surface  heat 
transfer  data. 


97 


REYNOLDS  NUMBER 


o 

o 

o 

o 

o 

o 

o 

o 

o 

0) 

to 

T-i  4J 

o 

o 

o 

o 

o 

Pd  tw 

o 

o 

m 

o 

•k 

rH 

o 

•* 

. 

«k 

o 

o 

to  CM 

Csl 

o 

ITl 

CM 

C“4 

•k 

to 

in 

r-i 

CNl 

IC| 

m 

r-i 

0) 

B 

0 

II 

II 

0 

« 

£ “ 

'O 

XI 

•o 

•u 

0) 

0) 

QJ 

01 

<U 

0) 

o 

Pd 

Ct! 

a. 

pd 

Oi 

pd 

O o 0 o +^J  x 


tn 

U) 

0) 

(0 

0) 

a) 

3 

3 

•H 

pC 

00 

3 

00 

3 

o 

S3 

(2 

(Hi 

rH 

h* 

0) 

a 

3 

<D 

<0 

>rl 

•H 

Pd 

0) 

«5 

O 

a 

M 

M 

0 

Ja 

00 

0) 

H 

00 

■a 

B 

3 

3 

| 

O 

3 

O 

OS 

0Q 

r*. 

i 

Fig.  4.1  Rough  surface  damping  constant,  A , versus  roughness  Reynolds  number 


Fig,  4.3  Rough  surface  mixing-length,  £ , versus 


BLOWING  FRACTION 


7 


yy 


*-R) 


a/ 


5 6 7 


Fig.  4.9  Rough  pipe  heat  transfer  predictions  with  smooth  surface 
turbulent  Prsndtl  number  model 
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CHAPTER  V 


SUMMARY  AMD  CONCLUSIONS 

The  principal  results  and  conclusions  of  the  present  rough  surface 
experiments  can  be  summarized  as  follows: 

1.  The  present  study  has  used  a deterministic  roughness,  0.050  inch 
diameter  spherical  elements  arranged  into  a most-dense  array.  The  ther- 
mal and  hydrodynamic  performance  of  the  boundary  layer  on  this  surface 
was  quite  likely  affected  by  its  highly  regular  roughness.  In  particular, 
the  marked  absence  of  any  'transitional'  behavior  may  be  due,  in  part,  to 
this  effect. 

2.  The  range  of  roughness  Reynolds  numbers  u k /v,  for  the  no- 
blowing  tests  was  from  24  to  200,  covering  the  transitional  roughness 
range  and  extending  well  into  the  fully  rough  range  by  conventional  defi- 
nitions. The  fully  rough  state  seems  to  have  been  attained:  it  certainly 
was  for  the  heat  transfer  data  and  probably  was  for  the  friction  factor 
data.  The  measured  Stanton  numbers  are  not  dependent  on  free-stream 
velocity,  but  are  functions  only  of  the  enthalpy  thickness  of  the  bound- 
ary layer  and  the  bloving  fraction,  p v /p  u . The  unblown  Stanton 
number  data  for  all  five  velocities  tested  lie  within  a few  percent  of 
each  other  when  compared  at  the  same  enthalpy  thickness.  With  blowing, 
the  uncertainty  in  the  measurements  Increases  but  the  tendency  to  be 
independent  of  velocity  remains,  even  at  very  high  blowing. 

These  comments  seem  also  to  apply  to  the  friction  factor  behavior, 
but  with  less  certainty,  <*lnce  the  data  are  more  sensitive  to  errors  and 
show  more  scatter. 

3.  The  behavior  of  the  boundary  layer  at  low  roughness  Reynolds 
numbers  was  Investigated  by  an  additional  test  with  the  boundary  layer 
thickened  by  blowing  through  the  first  two  feet  of  the  eight  foot  test 
section.  This  action  reduced  the  roughness  Reynolds  number  to  14  at  the 
downstream  end  of  the  plate.  Even  at  this  low  roughness  Reynolds  num- 
bers the  heat  transfer  data  continued  to  exhibit  'fully  rough'  behavior: 
Stanton  number  was  the  same  function  of  enthalpy  thickness  as  observed 
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at  higher  roughness  Reynolds  numbers.  It  remained  independent  of  stream 
velocity. 

The  expected  approa:h  to  smooth  plate  behavior  was  not  observed. 

If  smooth  plate  heat  transfer  represents  the  asymptotic  behavior  of  the 
rough  surface  as  the  roughness  Reynolds  number  is  reduced,  it  was  not 
apparent  from  these  experiments. 

4.  Blowing  through  a rough  surface  diminishes  the  Stanton  number. 
The  effect  of  blowing  is  predictable  using  the  same  relationship  found 
valid  in  smooth  plate  studies: 


where  St  and  StQ  are  evaluated  at  the  same  enthalpy  thickness,  A 
(for  a smooth  plate,  the  Stanton  numbers  are  to  be  evaluated  at  the  same 
enthalpy  thickness  Reynolds  number) . 

5.  Transition  on  the  rough  surface  began  at  approximately  the  same 
momentum  thickness  Reynolds  number  (350-450)  for  all  of  the  conditions 
tested.  For  the  surface  tested  here,  transition  occurs  at  about  the 
same  momentum  thickness  Reynolds  number  that  would  be  expected  for  tran- 
sition on  a smooth  plate.  The  transition  moves  upstream  with  increasing 
free-stream  velocity  and  with  blowing. 

6.  Blowing  in  a region  which  would  otherwise  have  remained  laminar 
may  greatly  increase  the  local  heat  load.  The  turbulent  heat  transfer 
coefficient,  even  with  moderately  strong  blowing,  is  much  larger  than 
the  laminar,  unblown  value. 

7.  The  effects  of  roughness  on  skin  friction  can  be  accounted  for 
in  a finite-difference  boundary  layer  prediction  scheme,  using  a mixing- 
length  model  for  mean  field  closure,  by  modifying  the  mixing-length 
variation  very  near  the  wall.  Successful  prediction  of  rough  surface 
heat  transfer  requires  modification  of  the  turbulent  Prandtl  number  dis- 
tribution very  near  the  wall. 

In  conclusion,  it  is  safe  to  say  that  the  work  reported  here  is  just 
the  beginning  of  the  rough  surface  investigations  at  Stanford.  Already 
under  way  is  an  extensive  program  to  make  detailed  hydrodynamic 
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measurements  in  the  rough  surface  boundary  layer.  When  complete,  these 
should  provide  a much  better  basis  for  formulating  realistic  prediction 
models. 
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APPENDIX  A 


MAIN  AIR  SUPPLY  SYSTEM  COMPONENTS 


The  following  Is  a brief  description  of  the  construction  details  of 
the  major  components  in  the  main  air  supply  system. 

A.l  Inlet  Header  and  Screen  Box 

The  heat  exchanger  inlet  header  and  screen  box  are  shown  with  their 
overall  dimension  in  Fig.  A.l.  One  side  of  the  box  is  removable  to  allow 
removal  of  the  screens  and  screen  spacers.  The  screen  material  is  HO 
mesh,  .0065  stainless  steel  wire  stretched  on  1 x 4 inch  clear  pine  frames. 
Four  such  frames  are  used  with  2x4  inch  clear  pine  spacer  frames  be- 
tween each  screen.  Both  sides  of  each  screen  and  spacer  frame  are  lined 
with  felt  gasket  material  to  avoid  leakage  around  the  screen  pack.  Align- 
ment pins  are  used  to  insure  that  the  screen  and  spacer  frames  are  properly 
aligned  in  the  box.  A photograph  of  the  header,  heat  exchanger  and  screen 
box  assembly  before  the  nozzle  was  installed  is  shown  in  Fig.  A. 2. 

A. 2 Wind  Tunnel  Nozzle 

The  Roughness  Rig  nozzle  shape  is  shown  in  Fig.  A. 3.  The  nozzle  was 
fabricated  in  two  halves  from  fiberglass  and  then  hand  fitted  together. 

The  nozzle  wall  shapes  were  first  layed  out  on  1/16  inch  aluminum  tem- 
plates from  which  the  mold  for  the  fiberglass  build-up  was  made  using 
high  density,  closed  cell  polyurethane  foam.  The  approximate  nozzle 
shape  was  carved  from  a block  of  polyurethane  foam  with  a hot  chromel 
wire  using  the  aluminum  templates  as  guides.  The  final  shape  was  ob- 
tained by  hand  sanding  the  foam  block  and  then  reinforcing  its  surface 
with  a thin  layer  of  polyester  resin  which  was  sanded  to  conform  to  the 
final  nozzle  surface  shape  and  surface  smoothness. 

Before  the  fiberglass  layup  was  made,  the  mold  was  waxed  and  sprayed 
with  a release  agent.  The  first  coat  of  resin  was  a white  gel  coat, 
sprayed  on  the  mold.  This  was  followed  by  two  layers  of  fiberglass  mat, 
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a layer  of  fiberglass  cloth  and  two  additional  layers  of  fiberglass  mat. 

An  additional  layer  of  mat  was  added  to  the  nozzle  flange  for  reinforce- 
ment. 

Reinforcing  ribs  made  from  one-eighth  inch  thick  masonite  were  at- 
tached to  provide  added  strength  and  the  nozzle  halves  were  bolted  to- 
gether along  the  horizontal  midplane.  The  final  shape  of  the  nozzle  sur- 
face conformed  to  the  design  shape  within  .020  inches  total  run  out. 

Figure  4.4  is  a photograph  of  the  two  nozzle  halves.  The  male  nozzle 
mold  Is  shown  in  the  background. 

A. 3 Multistage  Diffuser 

The  Roughness  Rig  diffuser  is  a multistage  vaned  diffuser  unlike 
any  previously  built  at  Stanford  in  that  several  two  dimensional  diffusers 
are  arranged  in  series.  Figure  A. 5 shows  a sketch  of  the  diffuser  with 
its  overall  dimensions.  A photograph  of  the  diffuser  before  it  was  in- 
stalled in  the  rig  is  shown  in  Fig.  A. 6.  The  first  stage  of  the  diffuser 
has  a movable  top  which  can  be  adjusted  to  align  itself  with  the  test 
section  top  and  form  a smooth  transition  at  the  test  section  exit.  Fol- 
lowing the  adjustable  inlet  section  is  a two-dimensional  expansion  which 
employs  five  vertical  vanes  to  expand  the  diffuser  width  from  20  to  24 
inches.  Next  are  two  separate  two-dimensional  expansions  employing  five 
and  then  eight  vanes  to  expand  the  diffuser  height  to  24  inches.  This 
empties  into  a plenum  box  connected  to  the  blower  inlet  with  a flexible 
connection.  As  can  be  noted  from  the  photograph,  the  diffuser  variable 
area  inlet  section  and  first  vaned  expansion  has  been  constructed  from 
plexiglass  with  aluminum  vanes.  At  the  bottom  of  the  plenum  box  is  a 
flange  for  connection  to  the  small  charging  blower  used  to  control  the 
tunnel  static  pre isure  level.  The  estimated  efficiency  of  this  diffuser, 
had  each  of  the  stages  performed  in  the  assembly  as  it  would  have  in- 
dividually, would  have  been  in  excess  of  502.  The  actual  diffuser  ef- 
ficiency is  closer  to  402  and  even  less  at  lower  test  section  velocities. 
This  performance  is  however  acceptable  and  the  wind  tunnel  does  achieve 
its  test  section  design  velocities.  Part  of  the  difficulty  with  the 
diffuser  performance  is  that  the  inlet  flow  into  the  diffuser  is  almost 
all  boundary  layer  fluid,  particularly  when  there  is  blowing  in  the  test 
section.  This  contributes  to  its  reduced  performance. 
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Constant  current,  hot  wire  type  flow  meters  are  used  to  measure  flow 
rate  in  each  of  the  individual  plate  transpiration  air  supply  lines. 

These  flow  meters  consist  of  two  side-by-side  probes  installed  in  the 
center  of  a 3 foot  long  tube  used  as  a riser  between  the  transpiration 
header  box  and  the  control  valves.  Each  unit  is  individually  calibrated. 
The  probes  are  small  tubes  which  span  the  delivery  pipe  and  contain  the 
junctions  of  a differential  iron-constantan  thermocouple.  One  junction 
is  in  a thin  wall,  0.025  inch  diameter  glass  tube  which  is  a passive 
probe.  The  other  junction  is  in  the  center  of  a heated  rod.  The  current 
ir.  the  heater  is  held  constant  during  flow  meter  operation  and  the  dif- 
ferential thermocouple  output  is  the  signal  used  to  measure  flow  rate  in 
the  tube.  A schematic  of  the  circuit  is  shown  in  Fig.  B.l.  The  heated 
probe  consists  of  a #2-56  nylon  screw,  center-drilled  with  a #50  drill. 
Each  end  of  the  screw  is  threaded  into  a short  piece  of  #8  brass  screw 
stock,  which  serves  as  an  end  terminal  for  the  heater  power  connections. 
Nichrome  heater  wire  is  wound  onto  the  nylon  screw  using  the  threads  to 
space  the  heater  wire  and  is  soldered  at  each  end  to  the  brass  terminals. 
The  heater  wire  is  then  potted  in  place  with  Permabond  cement.  A photo- 
graph showing  the  soldering  of  the  heater  wire  to  the  brass  end  terminals 
on  a flow  meter  probe  is  shown  in  Fig.  B.2.  Typical  resistance  of  a probe 
heater  is  5 ohms,  although  this  value  varies  a few  percent  from  heater 
to  heater. 

Each  flow  meter  unit,  consisting  of  the  riser  tube  with  differential 
thermocouple  and  heater,  was  individually  calibrated  on  the  flow  bench  In 
the  Thermosciences  Measurement  Center.  The  meters  were  calibrated  in 
both  upward  and  downward  flow,  for  blowing  and  suction,  using  a bank  of 
Meriam  Laminar  Flow  Meter*?  as  secondary  standards.  These  laminar  flow 
meters  had  been  previously  calibrated  against  standard  ASME  flow  orifices. 
In  the  controlled  environment  of  the  Measurement  center,  it  Us  felt  that 
the  flow  measurements  were  accurate  to  within  ± 1%,  the  accuracy  normally 
claimed  for  the  ASME  orifice. 
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At  each  flow  calibration  point,  the  heater  current  was  set  and  the 
differential  thermocouple  signal  allowed  to  stabilize.  A Hewlett- 
Packard  2401  IDVM  using  an  external  clock  to  allow  Integration  of  the 
flow  meter  signal  over  a 10-second  Interval  was  used  to  read  the  differ- 
ential thermocouples.  A typical  flow  meter  calibration  curve  is  shown 
In  Fig.  B.3. 

The  flow  meters  Installed  In  the  rig  can  operate  In  two  modes.  In 
the  first  mode,  the  power  and  differential  thermocouple  connections  are 
made  through  a series  of  double  pole  switches  such  that  flow  meter  opera- 
tion Is  restricted  to  a single  channel  at  a time.  Activation  of  each 
successive  flow  meter  heater  circuit  also  activates  its  read-out  circuit. 
The  flow  meters  are  operated  in  this  mode  while  setting  up  the  rig- 
running conditions.  In  the  second  mode  of  operation,  power  is  supplied 
to  the  meters  in  blocks  of  6 at  a time.  This  mode  is  used  during  data- 
taking.  A short  computer  pror'-ram  is  used,  called  FLOMET,  which  calculates 
the  flow  meter  signal  requir  d for  the  desired  transpiration  flow  for 
each  individual  meter.  The  flow  control  valves  must  be  adjusted,  one  at 
a time,  such  that  each  flc  / meter  produces  its  desired  signal.  The  same 
flow  meter  calibration  U jles  in  FLOMET  are  also  contained  in  the  data- 
reduction  program,  and  i ne  flow  meter  readings  are  re-recorded  when  the 
actual  heat  transfer  d t ta  are  taken,  so  accurate  flow  measurement  is 
assured. 

The  flow  calibration  tables  contain  flow  and  flow  meter  response  in 
terms  of  SCFM  and  EMFO  — that  is,  cubic  feet  per  minute  of  air  at  stan- 
dard conditions  and  the  emf  reading  the  flow  meter  would  have  produced  if 
air  at  standard  conditions  had  been  flowing  through  the  transpiration 
system.  To  convert  from  actual  to  standard  conditions,  the  following 
arguments  must  be  used.  The  flow  meter  differential  thermocouple  signal 
is  proportional  to  the  temperature  difference  between  the  heated  and  the 
passive  probes  of  the  meter. 


emf  a AT 


(B-l) 
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The  temperature  difference  can  be  related  to  the  heat  transfer  coef- 
ficient between  the  flow  and  the  heated  probe,  and  finally  to  the  Nusselt 
number. 


emf  “ AT 


(B-2) 


A similar  expression  can  be  written  for  a 
standard  conditions,  denoted  by  the  subscript 


flow  meter  operating  at 


emf 

o 


ec 


q d 
o o 

A k 
o o 


(B-3) 


Taking  the  ratio  of  the  two  and  noting  that  the  area,  diameter  and 
heat  flux  will  remain  the  same  between  the  two  conditions. 


M 1 

o Nu  k 

■ — ■ — . h — . ■ — ■ » — i 

emf  Ku  k 

o o 


(B-4) 


From  the  flow  meter  calibration  data,  the  Bignal  is  seen  to  vary 
with  the  0.45  power  of  the  flow  Reynolds  number.  Using  the  normal  C.3 
dependence  on  Frandtl  number,  we  can  write 


Nu 


Re0'45  Pr0’3 


„ 0.45  ....  .0.3 
■ Re  (k/yc) 


(B-5) 


Also, 


Re  “ SCFM 


(B-6) 


With  this,  holding  SCFM  constant  yields 
emf 


emf 


/SCFM  \°*45  / M°-45  (Kj^  c_\°-3  /k_\  /M°‘15  /c_\°-3  /U°-7 

" \scFmJ  \y  ) \k  yo  c J \kj-  \u  ) \cj  \kj 


(B-7) 

The  following  temperature  and  humidity  dependence  is  assumed  for  the 
air  properties: 


k/k  - (T/T  ) 
o o 


0.74 
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\i/Uo  - (T/To) 


0.725 


c/c  - not  temperature  dependent  in  the  range  of 
0 the  test  data 


u - absolute  humidity  of  the  air 


(k/k  ) - not  humidity  dependent  in  the  range  of  the 
test  data 


<u/yQ)  - (l  - 0.7  u) 


(c/cq)  « (1  + 0.9  oj)  . 


Substituting  these  into  the  expression  for  emf  gives 


/l_\0,41  (1  + 0.9  u)°" 
\To)  (1  - 0.7  w)°‘- 


(B-8) 


Using  binomial  expansions  for  the  humidity  terms,  this  can  be  simplified 


emfQ  - emf (T/T  ) * A (1  + 0.38  w) 


(B-9) 


All  flow  meter  calibration  data  have  been  converted  into  standard 
conditions  with  this  expression,  for  use  in  the  data  reduction  program. 
Flow  meter  readings  taken  at  test  conditions  are  then  converted  to  stan- 
dard conditions  before  entering  the  calibration  tables. 

Experience  with  this  metering  system  has  been  highly  satisfactory. 

At  relatively  high  flows,  10  cfm  and  over,  adjustment  of  the  valves  to 
give  the  desired  flow  meter  signal  is  relatively  easy.  At  lower  flow 
rates,  flow  valve  adjustment  is  a bit  more  tedious  because  of  the  slower 
response  of  the  meter  and  the  coarse  adjustment  that  the  ball  valves 
afford  in  their  nearly  closed  positions.  The  meters  themselves  seem 
stable  and  do  not  show  any  appreciable  drift  during  a data  run.  On  at 
least  two  occasions  during  the  rig  shakedown,  flow  meters  were  removed 
from  the  rig  and  recalibrated  tn  be  sure  that  the  heaters  and  differential 
thermocouples  were  not  'aging'  and  changing  with  use.  In  only  two  cases 


were  changes  in  the  meter  calibrations  noted , and  each  of  those  occurred 
because  of  overheating  of  the  heater  probe.  The  two  probes  had  been 
inadvertently  activated  when  there  was  no  flow  in  the  transpiration  sys- 
tem. In  all  other  cases,  there  had  been  virtually  no  change  in  the  flow 
meter  calibrations. 
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APPENDIX  C 


PLATE  POWER  CONTROL  CIRCUIT 


Plate  power  control  is  achieved  by  controlling  heater  voltage  using 
a solid-state  control  circuit  for  each  plate,  as  shown  In  Fig.  C.l.  The 
circuit  was  designed  by  T.  R.  Mogel,  of  the  Civil  Engineering  Department 
at  Stanford,  who  also  did  the  printed  circuit  board  layout  and  acted  as 
general  consultant  in  shaking  down  the  Roughness  Rig  power  control  sys- 
tem. The  following  description  is  paraphrased  from  a discussion  of  the 
control  circuit  written  by  T.  R.  Mogel  for  the  seventh  quarterly  progress 
report  on  the  research  contract  which  funded  the  Roughness  Rig  project. 

The  power  control  circuit  consists  of  three  main  parts:  an  opera- 
tional amplifier,  a high-power  output  amplifier,  and  a feedback  network. 
The  heater  voltage  is  set  by  a 10  turn  pot  with  +0.2%  linearity,  which 
sends  a proportional  amount  of  the  -10  volt  reference  to  the  heater  con- 
troller. The  voltage  appearing  across  the  heater  is  determined  by  the 
feedback  network.  The  whole  circuit  can  be  thought  of  as  an  operational 
amplifier  whose  D.C.  transfer  characteristic  will  be: 

v„ut  - - Vl„  VR3  <R1  ■ R2*  R3  ■ V <c-l) 

The  output  voltage  will  be 

Vheater  * ’ \ef  0 VR3  <C-2) 

For  the  resistances  used  in  the  circuit, 


V 


heater 


22.1  a 


(C-3) 


where  a is  0 when  the  pot  resistance  is  at  the  counterclockwise  stop 
and  when  1 is  at  the  clockwise  stop. 

The  capacitor  is  used  to  stabilize  the  amplifier  by  rolling  off  the 
high  frequencies.  A diode  is  used  to  prevent  the  output  of  the  operational 
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amplifier  from  going  negative.  A second  diode  is  used  to  allow  the  high- 
est of  the  heater  voltages  from  all  the  channels  to  be  known. 

A feature  of  this  system  is  that  the  heater  voltage  stability  is  a 
function  of  the  stability  of  the  -10  volt  reference  and  not  the  heater 
supply.  The  voltage  drop  of  the  lead  wires  of  the  heaters  is  compensated 
for,  since  the  amplifier  'looks  at*  the  voltage  directly  across  the 
heater. 

A shunt  is  placed  in  series  with  each  heater  element  to  measure  the 
current.  The  heater  voltage  is  also  measured  so  that  heater  power  can 
be  determined  independently  from  the  control  circuit.  The  accuracy  of 
the  data  is  not,  therefore,  in  any  way  related  to  the  operation  of  the 
control  circuit. 

The  circuitry  in  the  dotted  box  is  on  a printed  circuit  card  (6  chan- 
nels per  board) . 

A protection  system  has  also  been  designed  which  provides  for  shut- 
down of  the  heater  power  supply  in  any  of  the  following  cases:  over- 
voltage on  the  bus  bar,  over-voltage  on  any  individual  plate  or  over- 
current on  the  bus.  The  system  guards  against  exceeding  a total  bus 
current  of  750  amps  or  bus  voltage  of  22  volts  using  fixed-level  detectors. 
A variable-level  detector  circuit  is  used  as  a check  against  individual 
plate  over-voltage.  One  potentiometer  setting  guards  all  of  the  plate 
circuits  against  the  same  maximum  plate  voltage.  If  any  of  the  individ- 
ual plate  voltages  exceeds  the  set  value,  the  D.C.  power  supply  is  shut 
down.  Circuit  diagrams  for  these  protection  circuits  are  shown  in  Fig. 

C.2.  They  have  been  wired  on  a single  printed  circuit  board  and  are 
included  with  the  control  circuits  in  a common  card  box  behind  the  circuit 
control  panel.  A photograph  of  the  back  of  the  plate  power  control  panel 
is  shown  in  Fig.  C.3.  One  of  the  four  printed  circuit  boards  is  shown  on 
an  extender  card  above  the  card  box.  The  bundle  of  wires  at  the  bottom 
of  the  box  Include  both  the  control  leads  to  the  power  transistors  mounted 
on  the  water-cooled  bus  bar  and  the  heater  voltage  sense  leads.  A photo- 
graph of  the  bus  bar  box  is  shown  in  Fig.  C.4.  The  heavy  cables  at  the 
rear  of  the  box  are  the  power  leads  from  th  • overhead  bus  bars. 
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HEATER  VOLTAGE 


Fig.  C.2  Plate  power  control  protection  circuit  schematic 


APPENDIX  D 


POROUS  PLATE  FABRICATION 

In  order  to  evaluate  the  practicality  of  brazing  the  poroua  plates , 
the  brazing  of  ' logpiies'  made  c * segments  of  copper  wire  was  attempted 
first.  This  was  done  by  plating  (with  electroless*  nickel)  several  1/2 
inch  long  pieces  of  0.050  inch  diameter  copper  wire  which  were  then 
placed  into  j stainless  steel  mold  and  heated  for  seven  minutes  at 
17o0°F  in  a hydrogen  atmosphere  furnace.  The  result  was  a success:  a 
brazed  'logpile'  of  copper  wires.  With  this  established,  a search  was 
begun  for  a supplier  of  copper  balls  for  the  porous  plates.  A 5 lb. 
sample  of  copper  balls  was  obtained  from  the  Pioneer  Steel  Ball  Company 
of  Uniouville,  Connecticut.  Finding  a plating  contractor  to  plate  the 
balls  proved  more  difficult  than  had  been  anticipated.  A part  of  the 
first  sample  of  copper  balls  was  sent  to  the  Naval  Weapons  Lab  Plating 
Shop  at  China  Lake.  The  remainder  were  distributed  among  several  local 
plating  shops  who  were  asked  to  bid  on  the  job  of  plating  100  lbs.  of 
the  balls.  As  each  sample  was  returned,  a small  plug  was  brazed  in  the 
Hansen  Labs  furnace  using  a quartz  tube  as  a mold  to  test  the  quality  of 
the  braze.  The  ball  sample  sent  to  China  Lake  was  used  to  assemble  a 
A inch  square  plate  for  the  first  large-scale  demonstration  of  the  braz- 
ing technique.  The  mold  for  this  sample  was  constructed  with  OFHC  side 
rails  and  a sheet  of  1/8  inch  quartz  for  the  bottom.  The  advantage  of 
the  transparent  bottom  was  that  the  ball  rrrangemsnt  could  be  examined 
to  ensure  the  ball  array  on  the  bottom  of  the  mold  had  not  been  disturbed 
us  successive  layers  were  added.  There  was  some  hope  initially  that  if 
the  balls  were  carefully  shaken  or  jarred  they  would  arrange  themselves 
in  their  most  iense  array;  but  after  considerable  experimentation  this 
hope  was  finally  abandoned  and  the  balls  were  arranged  by  hand,  one  layer 
at  a time,  until  the  desired  thickness  was  achieved.  The  A -inch  mold  was 
successfully  fired  in  the  Hansen  Lab  furnace. 

Work  on  four  full-size  mold  cavities  was  then  begun.  High  cost  and 
slow  delivery  for  the  quartz  bottoms  mitigated  against  their  use  for 
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these  molds  and  OFHC  copper  was  used  instead.  An  order  for  the  full  100 
lbs.  of  balls  was  placed,  and  Electro-Coatings,  Inc.  of  Emeryville, 
California,  was  selected  as  the  plating  contractor.  They  had  developed 
a small  stainless  steel  basket  to  tumble  the  balls  to  improve  the  plating 
uniformity.  Samples  were  then  obtained  with  various  plating  thicknesses 
and  used  to  determine  the  plating  thickness  which  would  provide  the 
strongest  bond  without  excessive  buildup  of  the  brace  material  at  the 
ball  connections.  A small  mold  which  produced  1 inch  square  samples  was 
used  for  these  test  braze  samples.  A plating  thickness  of  0.0005 
Inches  was  finally  selected.  The  plating  contractor  agreed  to  plate  the 
balls  in  9 lb.  batches,  which  was  just  sufficient  for  a single  4 x 18 
inch  plate  1/2  inch  thick.  It  was  found  that  difficulty  was  encountered 
If  a plate  was  assembled  with  balls  from  different  plating  batches. 

Small  differences  in  plating  thicknesses  resulted  in  small  differences 
in  the  ball  sizes  which  didn't  show  up  until  several  hundred  of  them  had 
been  arranged  in  rows  in  the  mold. 

The  firing  of  these  large  molds  was  not  possible  In  the  Hansen  Lab 
tube  shop,  but  space  and  furnace  time  were  provided  to  us  at  the  Stanford 
Linear  Accelerator  (SLAC)  tube  shop.  The  SLAC  tube  shop  has  a 24”  hydro- 
gen atmosphere  brazing  furnace,  one  of  the  few  of  this  size  on  the  West 
Coast.  The  tube  shop  supervisor.  Hr,.  Bob  Boscnburg,  and  his  people  spent 
considerable  time  educating  us  in  'brazing  technique'.  Keeping  the  molds 
and  balls  clean  and  use  of  the  proper  release  agent  to  keep  the  plates 
from  sticking  in  the  molds  were  only  part  of  this  help.  The  balls  were 
assembled  In  the  molds  by  hand,  filling  four  molds  each  time  before 
firing.  To  assist  in  this  phase  of  the  project,  two  temporary  employees 
were  hired.  The  ball  pack  used  consisted  of  11  layers  of  balls,  each 
layer  91  rows  wide  at  the  test  surface  with  357  balls  in  each  row.  The 
ball  layer  behind  the  surface  layer  contained  90  rows  and  all  other  layers 
used  89  rows.  This  construction  provided  a lip  at  the  top  surface  of  the 
plate  bo  that  when  two  plates  were  laid  side  by  side  they  would  form  a 
continuous  surface  but  there  would  be  space  between  them  beneath  the  top 
layer  for  the  phenolic  support  strips.  As  already  noted,  the  final  layer 
of  balls  was  arranged  so  that  a ball  row  wae  omitted  between  every  ten 


rows  to  form  grooves  for  the  heater  wires.  A photograph  of  a ball  assem- 
bly nearly  completed  in  a copper  mold  before  firing  is  shown  in  Fig.  D.l. 
Fig.  D.2  shows  a photograph  of  the  molds  being  filled.  As  shown  in  this 
phorograph,  most  of  the  ball  arrangement  was  done  with  small  brushes. 

With  some  practice  it  became  possible  to  literally  paint  the  ball  rows 
into  place.  Each  plate  contained  on  the  order  of  350,000  balls.  Assem- 
bly of  the  ball  pack  took  between  three  and  four  days  for  an  experienced 
worker. 

Another  problem  that  the  SLAC  tube  shop  personnel  worked  out  for  us 
was  the  correct  'heating  pattern'  for  the  brazing  of  the  full-iize  plates. 
The  heating  pattern  finally  adopted  was  to  heat  the  plates  to  84G°C  and 
allow  the  molds  and  plates  to  equalize  at  this  temperature.  This  appears 
to  be  just  under  the  melting  temperature  of  the  phopuickel  plating  on 
the  copper  balls.  After  equalizing  all  the  furnace  thermocouples  to 
within  a few  degrees,  the  plate  and  mold  temperature  was  increased  to 
910°C  and  held  for  ten  minutes.  During  this  finel  heating,  all  mold 
thermocouples  were  held  within  a ten  degree  band  to  avoid  uneven  heating. 
At  the  end  of  the  ten  minute  cycle  the  furnace  heating  element  was  re- 
moved and  the  steel  retort  was  air-quenched > Typical  cooling  times  for 
this  process  was  usually  in  excess  of  twelve  hours  before  the  plates 
could  be  removed  from  their  molds.  This  particular  heating  pattern  seemed 
to  optimize  the  brazing  process.  At  lower  temperature,  the  braze  fillets 
formed  between  balls  were  smaller  and  often  would  not  hold.  At  higher 
temperatures,  the  plating  seemed  to  evaporate  or  be  absorbed  into  the 
balls  and  there  was  less  available  for  the  braze. 

Figure  D.3  is  a photograph  of  the  arrangement  used  for  the  molds  on 
the  pedestal  of  the  24  inch  SLAC  furnace.  A photograph  of  a tube  shop 
technician  removing  the  heating  element  from  the  steel  retort  after  fir- 
ing Ik  shown  in  Fig.  D.4.  A photograph  of  the  finished  plate,  after  the 
copper  molds  were  removed  is  shown  in  Fig.  D.5.  The  black  material  on 
the  mold  is  the  release  agent  used  to  avoid  having  the  balls  braze  them- 
selves to  the  mold.  It  was  found  after  several  firings  the  copper  molds 
became  extremely  soft  and  lost  all  temper.  Inspection  was  needed  after 
each  firing,  before  reassembly,  to  ensure  that  the  side  walls  were 


straight  and  square.  Fig.  D.6  is  a close-up  of  the  plate  edge  after 
firing.  The  ball  layer  lip  is  clearly  shown.  As  described  earlier,  this 
lip  made  it  possible  to  construct  the  rough  surface  test  section  without 
interruptions  in  the  roughness  pattern  at  plate  joints. 

The  successful  fabrication  of  the  Roughness  Rig  plates  was  possible 
thanks  to  the  help  of  several  people.  Particularly,  the  author  is  thank- 
ful to  Nick  Andrews,  foreman  at  the  Hansen  Lab  Tube  Shop  who  first  helped 
me  successfully  braze  the  trial  samples.  To  Manny  Gill  of  Electro- 
Coatings,  Inc.,  who  worked  out  the  batch  process  which  allowed  the  suc- 
cessful plating  of  the  8 million  balls  that  were  used.  To  Bob  Bosenburg, 
Supervisor  of  the  SLAC  Tube  Shop,  and  his  people,  who  gave  us  space, 
furnace  time  and  the  benefit  of  their  many  years  of  experience  in  brazing 
the  full  size  plates.  Finally,  we're  thankful  to  the  two  temporary 
employees  who  helped  assemble  the  plates,  Jim  Burlison  and  George 
Zanetti,  whose  ball-stacking  ability  wan  surpassed  only  by  their  patience 
and  gpod  humor. 


ti>;.  D.2  Photograph  of  balls  being  arranged  in  the  molds 


Fig.  D.3  Photograph  of  the  molds  on  the  Fig-  D.4  Photograph  of  furnace  being  raised 

furnace  pedestal  before  firing  after  brazing  the  ball  assembly 


APPENDIX  E 


TABULATION  OF  EXPERIMENTAL  RESULTS 


This  appendix  contains  a summary  of  the  Stanton  number  data  taken 
in  these  tests,  along  with  a limited  number  of  velocity  profiles.  In- 
cluded are  Stanton  data  for  each  bloving  fraction  and  for  each  free- 
stream  velocity  tested.  For  each  of  the  test  conditions,  an  Initial 
velocity  profile  has  been  included.  This  is  typically  the  first  profile 
measured  once  the  boundary  layer  was  fully  turbulent.  The  purpose  of 
including  this  profile  is  to  provide  a 'starting  point'  for  prediction 
schemes,  for  those  wishing  to  use  these  data  to  check  a predictive 
method.  The  following  is  a summary  of  the  abbreviations  used  in  the 
data  listings. 
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Stanton  Number  Data  Listings 

VEL  Free-stream  velocity  for  the  test. 

T,  TO  Free-stream  static  and  total  temperature. 

RHO  Free-stream  density. 

P Free-stream  static  pressure. 

TDB,  TUB  Dry  and  wet  bulb  temperatures. 

PAMB  Ambient  pressure. 

VIR  ORIG  Turbulent  boundary  layer  virtual  origin,  determined 

from  the  curve  fit  to  momentum  thicknesses.  Measured 
from  the  test  section  inlet,  positive  values  are 
downstream. 


PL  NO 
X 

STN  NO 
MOM  REY  NO 
ENTH  REY  NO 
X REY  NO 
(X-XO)  REY  NO 
F 
B 


Plate  number. 

Distance  along  test  section,  from  inlet. 

Stanton  number. 

Momentum  thickness  Reynolds  number. 

Enthalpy  thickness  Reynolds  number. 
x-Reynolds  number,  measured  from  test  section  inlet. 
x-Reynolds  number  measured  from  the  virtual  origin. 
Blowing  fraction. 

F/St. 
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TPL,  TAIR  Plate  surface  and  transpiration  air  thermocouple 

readings . 

MONTH,  ENTTH  Momentum  and  enthalpy  thickness  for  each  plate. 

Momentum  thicknesses  vere  obtained  by  plotting  the  6 
or  7 measured  values  of  momentum  thicknesses  and 
Interpolating  to  find  values  for  each  plate.  En- 
thalpy thickness  for  each  plate  was  obtained  from  the 
Integral  energy  equation  using  the  measured  Stanton 
numbers. 


DIS  TH 


Displacement  thickness  from  measured  velocity  pro- 
files . 


MOM  TK 
SHAPE  FACTOR 
X-XO 

MOM  TH  (FIT) 


Momentum  thickness  from  measured  velocity  profiles. 
Displacement  thickness/momentum  thickness. 

Distance  from  virtual  origin. 

Momentum  thickness  from  least-squares  fit  to  momentum 
thickness  data. 


DIFF  Difference  between  the  measured  momentum  thickness 

and  that  from  the  curve  fit. 

CF2  (FIT)  Skin  friction  based  on  momentum  thickness  fit. 

CF2  SMOOTH  Smooth  surface  skin  friction,  taken  from 


.0125  Re 


-0.25 

6 


(E-l) 


RATIO 


Rough  surface  skin  friction/ smooth  surface  skin 
friction. 


ROUGH  RE  NO  Roughness  Reynolds  number. 


E.2 


Velocity  Profile  Listings 


F S VEL 
BLO  FR 
STAT  TEMP 
STAT  PRES 
PORT 
DIST 

Y - YTOP 


Free-stream  velocity. 

Nominal  blowing  fraction  for  the  run 
Free-stream  static  temperature. 

Free-stream  static  pressure. 

Plate  over  which  the  velocity  profile  was  taken. 
Distance  from  test  section  inlet. 

Distance  from  crest  of  rough  surface  balls  to  center 
of  probe,  with  no  corrections  applied. 


V/VINF  Local-to-free-stream  velocity  ratio. 
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4 14  U.0U2H6  6430  bU31  1033430  104b3b4  0.  U.  100.0  100.  7 U.UbO  O.I'bU 
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12  46  0.00222  16473  14433  43681.27  4630046  0.  0.  104.4  104.4  0.142  0.1211 
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24  44  0.0U2U2  50461  264 13104 7UU 721 12121 01  0.  0.  104.11  106.8  0.261  0.227 
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VELOCITY  PHOFlLtS  WITH  (SLOWING  AT  90  F T/SLC 


F S VCL 

VU  «olN  1/itC 

VO.«/t  1/StC 

09. /IF r/stc 

88.82F  T/SLC 

86.22F  T/SLC 

ULO  FK 

(1. 

O.OOl 
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0.004 

U * 000 

ST  AT  TLMF* 
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77 « OfjLGN 
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V/VINF 

Y-VIOP 

V/VlNf- 

Y-YTGP 

V/VINF 

Y-YTGF* 

V/VINF 

Y-YTUF 

V/VINF 

u. 

0, 

0. 

0. 

U. 

0. 

o,. 

0 . 

0. 

0. 

u.uu 

0.4  40 

0.011 

U.4S8 

0.011 

U * 429 

U.01 1 

0.3/b 

0*011 

0.374 

u.uu 

0.4S0 

0.013 

0.466 

U.013 

0.437 

0.013 

U • 30H 

0 • U 1 3 

0 . 392 

U . U 1 b 

0.462 

O.Qlb 

0.461 

0.01b 

0.4b2 

0 • U 1 b 

0.400 

0.01b 

0.40/ 

U . U 17 

0 . 4C>ti 

0.017 

0.494 

0.017 

0.464 

0.017 

0.411 

0 . U 1 7 

0.417 

U.llio 

0.482 

0.019 

0.606 

0*019 

0.4/4 

U . 0 1 9 

0*416 

0.019 

0.424 

U . U«!S 

U.SOl 

0.021 

0.S17 

0.021 

0.464 

U.U21 

0.432 

0.024 

U .449 

U.OJU 

0.623 

0 • 026 

0.S36 
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0 . 496 

0.026 

0.444 
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O.Ujb 
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U.  02b 

O.bOb 
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0.46/ 
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U.bl4 
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0.467 
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0.693 
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O.bOb 
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U • b 1 9 
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0.U49 
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U.UbD 
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U.Ubl 
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0.0 /b 
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0.6/7 

U.UOU 
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0.0  71 
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0.0/6 
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0 . 0 79 
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APPENDIX  F 


LISTING  OF  DATA  REDUCTION  PROGRAM 

The  following  pages  contain  a source  listing  of  the  data  reduction 
program  used  to  reduce  the  Stanton  number  data.  The  program  is  written 
in  Fortran  IV  and  employs  fixed-field  input.  Output  includes  the  re- 
duced data  and  a listing  of  the  raw  input  data  as  well.  Input  require- 
ments and  card  formats  are  documented  with  comment  cards  in  the  listing 
of  the  DATIN  subroutine  of  the  code. 
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,A1A  REDUCTION  PROGRAM  FOH  pUu&HNES<-.  RIO  StANIOM  NC  « ENEH6Y  UAL  PU>  S 
OlMEI'SlUN  CONI)  <24 1 .RK24I  .TCAST  (24)  »f TIN(24)*ul  1N<24) .LLLST ( ] 6) 
CuMON/UAI/TAtR.PSTAT.  A6.PUOX.TPL(24)  . ICS  1 12)  • SYR ( 24  I .CFM<24) 

1 . IaH(24)  ,TCFM(24)  ,EH(24)  .ES124)  ,KNtXT.M0REC»X2E.R0.U;iT.FM0MT(24) . 
XSTo(24) ,REX0<24) .ITAB. t AMU < 4 ) . 1 PUT , If R, IRHOH 
CuMM0N/UA0/mEX(24),STN(24>  .OfE(?4) 

COMMON/MU/TAM8TI4) . VA1H . T INF  , Uf  KO.  R.  iOINF  .PSl 
common/ amu/pamr . tub  . Tuo.uEni i it> ) > ki  . iEnbal. vah 
CummoN/AFLO/XXX.SCFM 

C LOuuUCTlyN  LOSS  CONSTANTS  > EUALUAIE  8/23/73 

DATA  C0ND/.145, . lbO> .1  JO. .U90» . 1 00 . . 2B0 » . IBS. . I90 . • 149 . . 1 SO > . 1 31 > 
1.16b..l77».172,.124i.l21..0U2*.|bO».l2o».12S**l2b».lOb..ll4..14?/ 

c shunt  resistances,  calibrated  4/73 

DATA  HI /3,044.2.96b  >3.036 >2.983.3. 1)3? *2. 928,2. 998.2.984. 2. 9BS> 
12.984>2.990.2.o84,2.S6u.2.9b9<2.Pfcbi,-'.9b7>  3.01S.2 .9d0 .2 .967. . 9A4. 
22.9b6>2.9Sb.3.o23.3.02u/ 

E CONSTANTS  FOR  CORRECTION  OF  AlR  INLFt  TEMP  FOH  CAST  I No  TEMP  MISMATCH 
DATA  AT  IN/  .30.  ,21  • .lb, . 14  • .30  . .26.  .Sr,.  ,20>  .29.  .24.  .38,  .lb  1,4  7,.  3b. 
1.27. .40. .25 1 .17> . 36. .2b< .18. ,32. .20. ,24/ 

DATA  UTIN/.097..1S0’  .Ob4  ,.194*  ,0S5».('.37>  .02S.  • USB# .04u < .060 . ,121. 
l.u/Ui  -076i  .040,  .029. . 0 /6 , . 038 • .uS7 . . ,l3u > . OU0 . . 02b. . J43 . . uSO . . U20/ 
IHROH=-l 

24  call  uatain 

C Hi.  is  AUS  HUMIDITY,  LBS  aATER/LB  air 

Hiis.01969l4-.OOOfc7952b*T28AU.91bb7E-nb*T»ib»*2-.000232»t  lUB-T  will 
HkeS-HAMB*8S0 . 74/ 12 . 

OEuCOHsl  .-.OO'/SMTOB-SU. ) 

IF  ( (Rl.iiT.31.ANO.  <Kl  .LT.b)  )PSTaT=0. 

PSTsPSTAT/12.*62. 3*0ENC0H 

IF ( (K1.E0.3) .0H.IK1.E0.6I )PU0X=0. 

PUX=PB0X/12, *bp *3*UENC0H 
LPpU.241«(1.*,9*HHI 

C Lol  is  A CU,/STadT  USED  IN  CALC  of  EL*  Cl  POaER 
Cul-J.41 5/ibOO.* (144 ./14b. I 
Pr4:0«71»»0,4 
UU  9 J=1.4 
N=i*(J-l)*l 
1 ib» ( 0-1 ) *1 

IF  ( (Kl.NE.AI  .AND.  (Kl.LT.6I)  CAl-L  EMF  T ( TAM*) I J)  > I AhdT  ( j ) < 613) 

It  I (Kl.tU.J)  .OH.  (Kl.EQ.b)) IAMBI  I J)=TmB 
CASl=(TCS(N*l)*TCS(N)l/2. 

CaS2=ITCS(N*1I*ICS(N42) >/2. 

Las3sICS(N)*(TcS(N!-CAS1)/2. 

Call  EMFT(CAS3,TCaST(I) .»b3) 

CAS4F(TCS(NI*CASl)/2. 

CaEL  EMFT(CAS4,TCAST(I*l)»*b3> 

LASS— ( CAS1 ♦ 1 CS (N*l ) )/2« 
call  LMF  T (CASS. TCAST (1*2) *663) 

CASbs<TCS(N*l)*cAS2)/2. 

C„LL  EMFTICASb.TCAST (1*3). bb 3) 

CaS7=<CAS2*TCS(N*2) )/2. 

CALL  EMFT(LAS7.TCAST(I*4).*63) 

CASd=.ICb(N*2)*(TCS(N*2)-CAS2  l/t. 

4 call  LMFT (CASH*  TCAST ( I *6) .663) 

c calc  of  free  Stream  temp  from  ic  RtA^lt.o 

CALL  EMFT(TAIH,TINF.*6J) 

Ki.0INF  = lPHtS‘PST)/(63.33»<Tll)F*4S9.6)  ) • 1 1 >-.2b*HH) 
lFIIKl.OT. OI.ANU.IKi.LT. 6))  GO  10  12 
RCF  = . Bb 

WAIH=SORT(VAH»SG*32.1T»62.4U/(o.»Rh0iNF) ) 
YtLH0=VAlH*VAlH/(77B.«b4,34*CP) 

L Cua'H  CT ION  OF  FREE  STREAM  TEMP  FOR  Vi  L HE*0 
T2tH0=T INF* l 1 .-HCF ) •VCLMU 
III,)  =TINF-HCF«vtLMU 
RhuOLOpRHOINF 

RnUlNFF(PRtS*PST)/(b3,33*(Tl)lF*4S9.6))*(i»*.2b*HH) 

PS1=(PHLS*PST)/144. 

VlSlNFP(li.O*.OlTtj*TINF)/<l.tb*RHOINF) 
l constant  for  retnoeds  no  calc 

UA1UsVA1r*SURT(RH00L0/HH01NF) 

CunS=VAIR/vIS1nF 

c constant  for  Stanton  no  calc 

EU2-0 . S*CP*RHO InF* VA IR 

12  IF  (( IOOT.EU.O)  ,ANU.  ( IRrOH.EE.  0)  ICAuE  HEAL) 

IM  ( IOUT.GT  .0)  ,ANq,  l IHROR.Lt.O)  )CAeL  MCAl)2 1 IOUT , x2e'(0«U2) 

So'MSO* 

Fe=1 .73. 

) i.07 1 . / 6 • 

*4UJ=A2ER0/U. 

Dw  b N- 1.24 

c Plate  poaER  calculation 

PmULN=C01*CS(N)*EH(N)/HUN) 

CrLL  EHfT(TAR<N><T.*S3l 
CALL  tMFTt TPLINI »TPLATt.»b3) 

IF  lxl.EU.31UU  lu  SOI 
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c 


l 


L 


l 


L 


c 


c 


l 


4 


L 


L 


CALC  OF  TUANS  FLO*  FHOd  FLO  MEiL'H  tMe.KUl  dLU«,Kls2  LUCK 
CALL  EMFTtTCFM(N).TT,*b3> 

Hho=lPFJEb*PHX)/(b3.J3*(TT  + 4b9.b))*(l.-.2o*HH) 

C«LL  FLOMIF'AMQ,  TT  .KM.  HllO.Kl  i N>  CFM  (At) , ACFMPLF'LO) 
cuKWfccriou  OF  aim  tlmp  foh  Cast  to  ai«  ump  milmatcn 
IF  (.J.Lt.9)  TTAMs.7b«TAMDT  ( 1 ) ♦ , 2b*  T AMD  t (21 
IF  ( (M.G1.9I  ,ANO.  (IJ.LT.lt,)  1 T TAM= I TAMHl  (2)  Ft  AMU  I lx)  wo. 

IF  lN.GL.16)TTAM=,7'>*TAMHT^3)■^.2‘^*^A^,  T(4) 

TlN=TT*(TTAM-TI)*U.-tXP(-1.17b/ACF«)l 
1 = 1 + ( A UN(iJ)-Bl  IN(IJ)  * ( ( IN-1CAST  (N)  ) ) 

FiAU  FHOM  FHoNT  OF  PLATtS 
bOi  IF  ( (ILNUAL.tO.l ) .AND.  (M.NE.3) ) GO  To 
HL0.'J=U.1p3 

1 KAijF  4s(llXtH0*4b4«b)/l(IU«)**4 

GO  TO  42 

41  HLoU=d  • 199 

TliAijF  4:  ( ( T Ju*4b9.b  I 2100,  ) **4 

42  0FK,jT=HC0N*.17I4*(  ( I TPlA  r£*4b9.f,  1 /idr . 

COjUOCTION  LObbfb  F POM  PLATES.  LflUOdt 

GLO-4l.'=COUOH.)»  ( TPLA  It- 1C  ALT  (II)  I / lOUO, 

UAS4=(  ( TCAbT  (I|)44‘j9.6)/1UU.  I *.4 
T A 1 1,4=  l ( T F4S9.6  )/  lo  U,  ) * * 4 
1P(k1.LT.3)uO  TO  21 
Kl;J  110  THAUb 

1 PLO=rPLATL*JbOO.*F'POFN/(4l),»lt, ) 

TPl  'AS  ( ( TPLt(*4',9.G)  /10U,  ) **4 
1ii(.4=<  ( <TPLu*T)/2.*459,b)/10P.)«*4 

H»U  ♦ COHO  FKO.,  HALK.fTO  THA.Ib.PL  A|  xVO  OF  Alh  AIJLI  FLAIL  UAU  T . V.P 
Up*CF.  = 12  ■ • .01b*  . b*  ( TPLd-T  I 

GoALK'OllAL^/iioO11*'  ■UU<J7*r,'U,4*,171,i*lHC4*-i4bb*ICAb4-.u3bU*lAIIJ4) 
t C0uV=0 . 

ACFA.SO. 

Sol " = U. 
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. )««4-T|JAdF4)/JbU(., 
( SOL  EXPlIJ, 


AXA-.il. 

PlFLOiO. 

1 Iso. 
bL0F=U. 

(ILL  (0)  =11. 

GO  TO  23 

21  F LAP-F’LP  LO*LP*  XoU0./(40.*l2.  I 
U (Kl.LT.2loO  TO  22 

KIsi  bLILKU,. 

T PLo=  rPLATt  *PH*  LO/  ( PLF  LO*CP  1 • ( l aP  t -F  EaF'J  I 
Ti'L04=(  ( TPL,I«4',9.0>  /10O,  ) **4 

HAp  FKOM  HACK  FOH  T.l  >Ck  • i|C  A T IIATl  ■ AlK  II  *..)• 

Go»cK  = . I 71'**(  .iHlb*TPLH4-.l4bb*TCAL4-.UiGU»TAI., 41/3  .00. 

k co-iV=CP*PLf  LO.lTPLb-TIUF  1 
Go  10  23 
K 1 si  ULO*  1 IjG 

22  I F'LjA  I PLA  TL-PHlL  It/  (PLF  LU*CP  )*(!•-£  xP(FlXPJ) 

Tf'Li  4 S | ( TPLJ*4..G.o)/1(’0.  | * * 4 

Ftp  FHpM  HACK  F pH  III  0* . nC  AT  ..JH  IF  P 

OuAp  ■*  - . 1714*1  .ni|G7»rF>L04«.13bkt*TAlH4-.l4b>*TCAj4)/3pU0. 
1 PO  :*  = CF‘*PLFL'I*  ( 1 t'LATk  - I I 

23  IF  ( iPACK.lt ,0.1  UpA  ,20. 

GLOGFUPIJ'jT * ,UAck*4C0UU 

G».k  at =ph:>l,  j-olos 

G...,  TSUHt  AT-k  CON/ 

IF ( lLNDAL.Lo.lluC  TO  l J 
IF  I'l.GT  . I lblOLe=iVj('J-i)  .HlOF 

bT  AH  ION  NO  CALL.  bAbLU  UN  MAG  TO  »A.  L TpFP  ['ll  > 
Sli.Mlsu  il  T/ICO.-*(TPLATL-T2LHO|  ) 

IF  (k  1 .Fu,  1)  G‘l  TO  lb 
HuoF  iALFM*M,.U/  IttU.  • .b*  ««ll<*Hii01i  I I 
HLl  (l.)=ULO»  /‘oTNlUl 
lb  IF  ('I.GT.l)oP  1(,  lu 

l .1  Ths.  . • (blUIN  | *hLOF  1 

bOHsCU.Nb^L'J)  Til/  12  . 

Go  Pi  3 

10  lLl'.=  tbT'.(:ji  FIlLOFFSTOLolbG, 

t.<1  r..=l  l(TTH*tL  1../2. 

bclx:SUM*LOI.L*tL  I M/24 . 

Fuo=FlG*fl 
3 FLOTsFi.g-FU.  u 

If  (FLGl  .lf.  IFLGT=Ij. 

HtMliISFLGl  »LONb 
IF (1HHOH.GT.OIgp  H b 

IF  ( (K  1 .III  ,31  .Oh  . (Ht»  (ill  ,LU.  0, 1 I GO  10  b9 
1 l*il=Il  AU*  1 
MOI  1 lAHISilNUu 
He  XO I IT  AU) =«EX ini 
»9  IhE  *0=F  M0M1  lUI»L0Nb/12.»U.b 
l * »-CONS*FL&*b .b 
1 XN-Ht  X ( 14 1 *0.3 
l.*LttJ=bUM*U,b 
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c 


c 


Ih  (HfcAUil  .u'J.O.  I OO  TO  7* 

M .j=.U^<*b/  (l‘K4*HX  * Hi)  ♦#U.D 

IHKl.UL  .3)Slb=IU0l'  / <L  AtMhLOt  /Slbl-l  • » 

SUUil=Mb/  4 i U'L/Wt*Vjy.u>/l  ULH(U4bV.O)  ) **.4 
Go  10  /b 
7*  bl  i-U. 

SlKlf))=n, 

7 b COulluUL 

H UOUI  .01 .1)1  MH  1 TL  (0*4  Hi*  U LlitjIMul  • 1KLMU*  1KLLN*  1 HI  **•  iKLXDULOX  * 
llU_L  ( N>  * JPLm  II  * I *XKOXJT  Hi)  *1  111  TM 

U UUUI .01 . 1 U'UfiUl  Ml  11*  II  LO*blNUil  * IHLXiU*  IHtLN*  l HL A A » IhLXM UU*  * 
IiIllINJ  . IKLmTI  *1*1  MOD  Hi)  *lUI  IM 
H I oi\MA  I l 1 3 * 14*1  ij  • b • t!  1 l»*  ,*  1 ?*  • I 7 . 4 ♦ X O « 3 * i'l  11 » 1 »2Mnil 

1M  10UI  .LU.OUKl  IL  U>*  MU*  bill  Hi)  *1  lil  TM*  IHlLN.X  MOM  l ( N ) * iHLMu*  ULo.  IH 
lLXA.I  Luf  * iKLXli.HtO*  • I3t  L Hi)  * f T * ft  AS  M IH  * 1MLA  II  * f * Hi.  ML  M*  LtO.JV  * QLOS  • U 
b i ot\MA I (13*1  ‘i.b»IL.3»lO*Xo.w*lO»14»lh»l  ti.4HiO»Hi.b»X0.3»4lu.l*3fH.b 
1#  !<*) 

00  10  o 

lb  l-»thC  = U.il  1/M|tf  AT  ♦ lUO. 

bl..LM  = bNLT/ < j,‘>***b»CP*«<Hollit  *V*1K>  *1UUUU. 

:)i.cM=0Ut  1/  ( c*' ♦ l I I ‘LA  1 L “ 1 I **U/b/OlJ.l 
II-  ( 1HH0H.0I  * II I oU  (O  11 

II-  ( l oui  *1.0.01  ttxl  It  C #«•  7)N*0HI  A 1 *1  LOflV’Oul  I *H  ML  * XXX  * bU  mi ULO1 • 11*11* 
1 LA  1 1 * I » TCAbl  (14)  *Ul  M*J  I iOHAI*  *OCoU‘*  T Hi*  SII  ILK 
7 X UKMAI  ( lb*  31-  111 Ui*x II.  3*  M t ♦ 3*41-  /.  1 * 3f  H.b.X  «*.  1 * Ko.«J  I 

II-  ( UHH.OI  • U ) KUl  1 L i O*  «? ) N#  liSH.fi*  LCoUV  * ULQb  * 01*1.  1 • HtHL  • Sltll-  h«  AtX  M.  TI'L 
1 A 1 1 * I »U 

IX  I iUUJ  *ol  • 1 IHUUUl  *?  *N*  I'KUL  *4*LwUI1V  * uLOO*  oUL  I • HLHC  » b 1 ULK  * ACX  M*  ILL  A I 


a • 1 *n 

w i Ui«*'.A  I ( 13*1  44 ♦ b * 31  b • b * l-tt • 3 • I 7.2*1  U.3.21  /.l  * 14) 
t>  L0..I1I4UI 

03  II-  ( l luiil  .L  4. CM  . AMU.  I IhKUii.Lt  .01  It  ALL  L)A  1 AO  l 

LIU  L*  ^ ) St,  I.  II  LMHOM  ANAL  t bib  lb  AbK  XOH#  ILK-OMIO  lLK=l»Ylb 
U (UK, l U.iMOU  10  1U4 
13  CJ..I1NUI 

H I l«l(oH.L,..l'io1  ll-bTfi(  1.  ) 

II  ( IKHbK.OT  .11  luLLb  I ( iHKOKUSTl  l-SHillD 
H l ( K i ,f  0.3)  .Auu.  I 1HK0H.U  . M)  >oO  10  l Du 
U 1 1*1.1  I .«,)  liHKOU.ol  *121)00  10  1JU 

OO  10  24 
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1 Uu  SLMK-U* 

(>0  101  1 = 1.  IKMoH 
jU  <*LKii  ^UML.n  4UI  Lb  H 1 1 KLLM  ( 1 1 
SU,  .1  H=b'4M  (SUMlHJ 
*n*  il  iu*  m 

H ( 1 Hk  jW  . Lv.  . K I 1 II  H.*  lll.HST  1 I * U'cLb  f ( 1 ) * 1=1  * lKMyl<  I ♦ bUMU. 
1*  I 1HKuM.Lu«1**)«N1  11  U * Iu3lbn  T*  luLlbl  U > * 1 = 1*  IHIMMI 
XUKMAI  U*1uA*  *UNILWIA1NI  1 A.iAUblb  bUMMAM  i I uK  LLATL  IJ*/ 
1uA*'I  LA1L  I-TAHIOI,  NU*4|«I  K =**l«f(l.s*/ 


1U 


loA*  • 
lOA*  * 

1 li  A * * 
luX*  * 
luX  * * 
1 u A * * 
luX*  ' 
li/X*  f 
1 LX  * * 


<U,S  III  01  11  ill  =•  *LJO.‘»*/ 

i • U 1 1 liil  ) =♦  *i  *j»i*,»*  A 

< • UdH  410*0  VA1HI  = • *IVO.b*  / 
<•01  XV  fLLAlL)  =**dU.b*/ 
4 • U 1 XiV  1 A 1 K I ibKMJ,*.*/ 
(.UU1  VOLI  1**1  f ='*LJU.b*/ 
(.01  MV  IbnUll)  =*.l*u.*>*/ 
l.ul  Vv  FLAbT  ) =•  •K/O.M/ 

S i At  it  (if*  flUN.Mt.H  LOKuWr  • •(  4,‘U.bl 


1U>  V Uu*».A  i ( / 1 1 u A * * VJfiLL  K I AlfiT  1 At. ALT  bib  bUM.IAHT  MU4  l*LAK  U**A 


1 


luX.'LLhlt  S 1 Ah  101*  4UMHIM  =**l<*U.b*/ 

IbX**  Ul.b  U4.0I  11  »D  = * *1.  *.0.b*/ 

lux**  I.Ul  Xw  TIul  I 2**Uj«'i*/ 

1 o A * * (.00,  Hii'O  V A lu  I = • »L*JO.b*  / 

1 U A * * (.01  Vv  1 l*L  A I L I -*»t:>U.b*/ 

1 O X * * (.01  MW  T A a 1 =**LLU.‘|*/ 

luX»*  (.001  VOL  1 i I'L  I =•♦!  ru.s*/ 

1 1 A * * (.01  MV  IbMUlill  =**l,'U.*»*/ 

lux*'  (.oi  Xv  Kami  — * * I fcU.s*/ 

lux*'  l.ul  MV  ll-LOfc/  =«*t^0.s#/ 

lux**  (*0."j  V X L OX’  lKl  = * •l«?U,*j*A 

lux.»  l.ul  XV  1A-MU  £ftl*.'U.W 

i u A * * (0.1  4n.'J  I'l'OAl  =«*X*»0.’»#/ 

lull*  blAUTOW  I.U.-UI  U<  LKKUM  = * *L*iU.b) 


1U4  IhuOKi-l 

(.MtL^  XOW  ,-OKL  L AULb  1.4  SLwUt.l«LL  * 1**1  X I - 1 00  UAtR  * KLAJ  bijMl 
IX  ( Klit  a I . t v..l  lyO  I u J4 
mK  1 ! L ( i* * 1 1 I 
11  XOi.MAlUHU 
b I wl* 


NyML 


LI«o 


C L)A  IAIN  bLliWOu  1 1 Sl  LiAIAIL  I Oi.  bl  UAT  A Ml  OOC  I lull  KHOoKAM 

'.OoKOUIlML  I At  AIN 

L TlUb  b'jMMOuTlf*L  HtMUb  114  UA I A M’OM  TIH.  MOUoMNLbS  KlO 

COKMOi4/liAl/1AlM*LbU(*  ^0*ljOUA*  ll*L  1^4  I . UbUi  1 *blK  4^4  I *LX«(^4  I 
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rrrrrrrr 


1.  Tan  (24)  ,TCRM<24)  .EH<24)  ,ES(24)  .HNEXT  .MOmU  , X2URU.NST  .FHuMT  I -.u  I 
2S)0(24),hex0(24).1TAB.TAMB(4),10UT,IeH.IHHqr  * * HUMTUl4)’ 

T0U'  T"tt'OENt  (la)  ,K1 . lEMBAL , VAH 

1 KHUN  ~ 1 KIaOH  ♦ 1 

IE ( IHHOH.tU.OlGO  TO  100 

' inU"v*ls  >0R  tWw0K  *««•»»* 

Tou:TUB*0.5 
00  To  54 

1 IK(IHHOH.6T.4)60  TO  2 
TUa=TUB-0,5 

YA1H:I AIHt.Ul 
Go  TO  54 

2 If (IHHOH.GT.JlGO  .0  3 
Ta1N=Ta1R-,01 
vAh=VAH*.002 

GO  TO  54 

a IE ( 1HHOH.GT.4IGO  TO  4 
VAHHV.AH-.002 
Tt  c 1 12 ) = .PL ( 12 ) t .01 
GO  TO  54 

4 IE  (lHHOH.Gt .5)00  TO  5 
U'L(12)=TPI.«12)-.01 
ThhU2)sTah(12)4.01 

Go  TO  54 

a IE (IRH0H.GT.6lG0  TO  6 
TmH(12!sTAH(12)-.01 
t><(  12)  =EH(  12  ) * ,001 
GO  TO  54 

6 IE ( 1HH0H.GT.7IG0  TO  7 
t.i(  1?  )=£!((  12  l-.OOl 
tSU2l=tSlU)-,0l 

Go  TO  54 

7 IE ( 1HH0H.G1 .fllGO  TO  8 
1. 5(12)  =£5(12)4. 01 
T(.S<bl=TC5lt.)4,ul 

GO  TO  54 

a 1K1HH0H.GT.4IG0  To  4 
105(0)=TCS(g)-.U1 
T(.(TH(12|RTCFM(i2)4.01 

Gg  10  54 

V li  (IHHOH.GT.IOIGO  TO  10 
10«(  12)=TIFM(  12)-, 01 
CE  m(  12)=CFm(12)4,Q25 
Go  Tl)  54 

10  1K1HH0H.GT.1UG0  TO  11 
CE  M < 12)=CFM( 12)-. 025 
T-ugl J)=lAM4(J)*,0l 

GO  TO  54 

11  CoNTlHUt 

TamK  JlsTAMul  J)-.01 

HoU«iHbO*«,l 
Gl  To  54 

loo  oui.t  i nut 

HLGlf.  MtlDl  ,G  INPUT  HtHE.  ALl  INPUT  t Ixfo  UtLO  1 Oh  AT  AG  IN.  lUTtO  utLi 
HtAU  INPUT.  ElHST  CARD  HAS  U LOL  UP  TjTLt  AMU  (.OL  .,0  IS  Al  luuICATOR 
1)  a 1 SI  STA'lTON  HUN  61 TH  OLOAlf.G 
Ie  41=2  STAuTon  HUN  *1  In  SULTIUN 
I)  41=3  STANTON  HUN  V1IH0UT  TRaNSPIHxT ION 
!►  41=4  ULOolNG  ENtHGY  BALANCE  HUN 
l)  4115  SUCTION  LNFHOT  BALANCE  RUN 
41=6  NO  .(Low  ENERGY  BALANCE  HUN 

HtA0(5,25)otNT,Kl 
25  FORMAT ( 18A4. 7x, 1 1 ) 
iE (4l.LT, 3)  GO  TO  Gl 

llAHsU 

UO  57  1=1,24 
SIHIIISO. 

siu( I )=U, 

57  HE  AO  1 11=0. 

IE  41  = 1 OH  2*  HtAU  TOH.  TwH  .PAMo  # TA1» , PSTa  t . VAH  .SG.PitOM , UH0  eORPAT 
ol  Ie (41«LT,J)HEAU(5.26)T06.Tgb,Pa4B*TAtH,PsTAT. VAH*  SG. P0UX 

Ie  M = J OH  h Head  TOB.T»H.PAua,TAIH.ESTA|.VAH.5o.  7F10  FORMAT 
lEUKl.ta.a).OH.(M.EO.b)  )NtAD(5.24)T0t(.I6B.PAM().  lAIH.I'STAT.VAH.SG 

IE  41=4  OH  i He.AU  T()B.T»(1>PAE'R,tAIH.pH0X.VAH,  giu  POKE  A I ^ 

Ie  ( IKI.&T.J)  .And.  (Al.LT. 6)1  HEA015.29 ) Tub.  I kU.PAPU*  | AIH.Puox.vaH 
T.0i*ET  OULU  TtHP  IN  DtGF  V 

iou=iiHY  outu  temp  in  uesf 
PAMB.-.ANii  PHLSSUHt  IN  tNCEEES  OF  MG 

,*1^=tNERGyTBA*ANCEMHl.'*S  "W  ^ ST*"t0n  *•«  ***»  *N  *V  F«* 

PS14T=TUNNfcU-ro-AMfl  PRESSURE  OlfFEHEf  CE  IN  INCHES  OE  H2o 
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O HtMu^  bimHotlTlNk  HOW  ST  DATA  HtOUtUOfJ  l'UOOKA.4 
SUuHOUUfJL  Ml  Al)2(  loUT  # A/t  HO.  • I 
CU,<f40N/HU/ I AM01  (4 ) # VA 1M # 1 INf  # l*#l  HC#HHolNf  *I»S1 
CU« .MON/AMH/PAMlM  TOli#  TAb#OLni  Uo>*KlilLUUAL»VAJ4 
If  (KIWaL.O)OO  10  A 
*Mi It io#  1U) 

io  iu..w*rt///*THt  ouit'ur  houunl  mouLbito  can  ul  ubi  l>  fuw  staihoii  mu 

INb  OH  LNLKoV  UAL  HUNS  ONL»  1 ) 

1001-U 
MLIUIIIJ  1 

1 If (Kl.Of.JIGO  10  J 

*N1  lt(..«Hlul  III  .VAIH#  1INI  »l/(  HU#hllOIIII  #l’bl>  tim.TWll'I’AMH.X/CHO 
If  I IOUI.GI  .lll’UNUf  biULHI  • VAIH#  I INf  • 1/LHOiNIIOIUf  il'bl  * TOO#  IWliif'AMU 
I#X#.LHO 

0 f0,<MAlUHl#IUA4//'  VCLS'if  b.  1 • ' f T/fiLC  ls'.fb.l#*  l)t(>f  TU='#fS.l 

I#'  DLGf  linos'  #1 1. .4#  • LUO/flJ  l's'#fG,2i'  f •»,  J • / • TIIHS'#fb.l#'  Ul  Of 
/ 1 oil-'  #f  b.  1 # • Ut  Of  I’AMIts' >f'*>.2#  • IN  NO  WIN  OHIO#  X04 * <f  <>•■! # * IN' 

J//<  I'L  X b TN  MOW  LNfH  X IX-XOI  f I)  H'L 

4 1 A lit  MOMfll  tNTIM'/'  NO  INLM  NO  HLYNU  HLYNU  KLYNO  lltYUO'tlO 

bX  > ' Ul  Of  UIGf  INCH  INCH') 

b fUKMAI  (111  » 1HA4//  ' VI. Ls ' #f  b.  1 # ' I T/!>LC  TS'ifb.l#'  DtGf  T0=  * * ♦ *» . 1 
1#'  III  Of  HnOS'.fu#4#'  LIlb/fM  I'='#f0.1'.'  I'bl'/'  rulis'ilb.li'  Ul  Of 
2 | .K-.  ..'-S.  1 • * ULOf  l‘AMli='.ff>.i!»'  IN  HO  VIII  OHIO#  *0='  «f  O.  it  ' IN' 

i//<  I’L  X brtl  MOM  LNTII  X (X-XO>  f H H’C 

4 1 A IK  WOW  1 1 1 I HUH'/'  No  INCH  NO  HL1N0  ItLYNO  IIL1N0  ULVNO'»lU 

bX.'HLOf  ULOf  INCH  INCH' I 

oo  ro  4 

1 COuTlNUL 
ILl.HALbl 

If (Kl#L0.4)KlbI 
If (Kl.LO.bJKl=H 
VAiKbVAH 

WHi  ftt o#  7I(IUNI>  IDI)'  I XIHI’AMIIi  VAIK 

7 I 0HWAT  1 1H1  # IHA4// ' fl)ll='.fb.I#'  ULOf  I WIJ- * *f  b.l#  * DLof  I’AMIIb'ilU 
l.L#'  IN  HO  VLLlblH  L KHl  s'  if  (.#  l # 1 IT/SLC'//'  I’L  HL  I'OriH  l COIIV 
/ LOUUL 0 tllff  I’tHCNf  bUI  (LOW  tl’L  TA1H  I’L'/'  NO  (II 

JU/bLC  UIU/bLC  IIIU/bt(,  UHJ/bLC ' # 1DX#  'LIIKOH  At!  M OLOf  ULOf 

4 HO' I 

If  ( I0U1 .01 . 1 1 HUNCH  'I.ULNT.  11)11#  rmilil’AMH#  VAIN 
VfOHMAKlM  *l«A4//»  Toils'  #fb#  I # ' Nt(.f  TWH-*  »f  b.  1'  • lit  (if  |'AWH='#lu 

I.id#'  III  HO  VtllblN  l,HH)s»»f(>#i( » f T/SLC ' // ' I’L  l’l.  I’OWK  t CONV 
'i  LobbLi  (lift  I'fHCNl  SIN  f 10.  TI'L  1AIH  I’L'/'  NO  111 

JU/btC  UTU/bLt  UlU/bU  li!U/bLC'#ll)».  'LHIIOM  ACfM  OLOf  Ul  Of 

4 HO'  I 
4 IILTUHN 
LHil 


4 OAI  All  bllllHOOlINL  IIAIAOl  fuM  bl  UA1A  lltOOL  I ION  I'HOlilUM 
bdilHoUl  INI  OAI  MU  I 

4 JNf’HI  uAIA  ANO  *#UMI  ’.UNION  NO  HAIlo'.  AHi  I'NINJfL  001  HlH| 

COWMON/UA  l/UlH#f'SU1  l ',(. . I’tMIX  . II’L  I 04  I > It I U I # Sf  N I /4  I . 4!  M ( ,'l|  | 

1 * UNI/.I.  I4IMI»'M)#UI(.’4I  .IbUMI.KNLXl.MOMtt  • M/LH4I#NI#I  #1  M41WI  um, 
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2Slo(24).HEX0t24).lTA».TAMa(4l.l0UT.IrR.IHM0R 
COMMON/OAO/IEXI24)  .S7W124)  H«'E!24) 

Cummon/ AMa/PAHe  . to*  . r w . dent  u*  i . k i m eubal  • v ah 

list  out  haw  input  DMA 

It-  I ( ItNBAL.EO.O)  .OH.  OU.EO.l)  >G0  To  *0 

WHITE  16. 17  I TWb,)Ob>PAMh.P0OX,TaIH,  (ThHhU  > < 1 = 1-*,) 

17  FurwaTUh  .'  HAW  INPUT  OATA  Fo*.  ENfit'GT  J.ALANCE  RUN*./ 

1 • TwU='.F6.1.'  DEGF  TU3-*  .F6,i . ' Dr$f  PAM(JS<  .H*.;).  * MiG  PROAS'. 

2Fu.2>'FH20  TAIR='.F6.3. ' MVV  • T*Nul=<  .FB.3>  • KV  TAM#*S« .( s.3»  > 

3MV  TAM03='.FS,3.«  MV  TAMt)4i*.F6.j.i  MV'l 
IF (Kl.tU.31wRlTEl6.61 

8 FoRMAT (•  Ft  No  T ‘U(MV J T4H(Mi;>  U'I  JOETsj  EUl-lV)  CAST 
1 Tc  TCS(MV)  M 

lMKl.NE.31*MlTe<6.9) 

9 FcHAATI*  PL  NO  TPLlNV)  TAP  (My  1 lUXIAVI  TFLulMV)  ullVOL 
IT  1 EC1MVI  CAST  TC  TCS(M;I  «) 

DO  2 N=l<24 
IF(N.GT,12>  GO  TO  3 

lMKl.EU.3!wHlfE(6.32IW.TFL  IN)  , TAR  (N)  »EH(N> , £SUj)  .N.TtbND 
lMKl«Ne,3)wHXTEl6.l51U*TPL(U> . TARUI)  * TCF"(M>  .LF'Mtw)  lEI'hil  .ESINl.W 

i.TCSlN) 

Go  TO  2 

3 1MK1.E0.3I  WHITE <6.39 1N.TPL1NI  .TAW  (0)  .(TWIN)  .ESIn) 

IM  Kl  *Nt  .3)  wHITEl6.90  IN.  TPL1HI . TAH  (fl) . ICFMI N)  .LFMI.NI  .fcUlli)  .C$1111 

2 Continue 

RETURN 

50  1MK1.NE.31  WHITE  (6.51  ITwB,  TCB.PAMbMAW'd  ) .TAWN12J  .TAl  Uli)  »fAMH(4 
1 1 . 1 A IH . VAH. PSTAT  >PBOX • sG. X2ER0 

51  FORMAT l 1H  .'  RAW  INPUT  OATA  FOH  STA'Dln  MIMUfcil  HUI.*/ 

1 • TwD=<  <F6.1> ' OEOF  Tlhs*  .Fa.  1 . * DEGF  PaFKF=  ' »1  6.2*  **HG* . 

2*  T»yBl  = *F6.3.  < MV  TAMB2S'  iFfc.3.  ' MV  TAMU3=*  ,lo,3, « Mv  TAMtlHi*  iFG, 
33,*  MV  * / • TAIH  = ' .F6.3.  < MV  Va  1H=»  ir7, 3,  ' *H2U  |>ST  A>  = * »F6.2.  * AHUCi 
4P(jo*  = 'P6.3.  *#H20  SPGR  = '.F6.3>  • X2fP0=»  .16.2.  'it  • 1 

lMtkl.E0.51  .ANU.  (IENHAL.Eu.fi  1 vHI TE(6.S4 1 
S4  For/aUIW  . * HAW  INPUT  OATA  F(.R  NO  THANS  SI  ANION  NUMSl'lt  HUN*  I 1 
IM  (kl.EU.3)  .ANU*  ( 1ENUAL  >Eu*  1 I ) V'HI  Tf  (6.  V ) 

52  Format < in  .»  haw  iupui  data  for  no  thans  i mekov  ha*  huh*) 

lMkl.EU.3l  WHITE(G.55|TwU,  TOB.PAMb.  T A 1 R,  V AR  , PSTA I »l>buX»  SGi  XKfc'HO 
Sb  FurMAM*  TwB='.F6,1«'  Ot&F  TUBS'  .FA.  1 > ' DEGF  FAMH-  • < f U . 2 * • V! «» * . 

1 /•  TAIH=*.F6.3.'  MV  VA1H=.  .t7,3.  "»W20  PS,  41  = * «Fu.2.  I/.II20 

2PuuX='F6.2.  <*H20  SPGH  = '.F6.3*'  X2F|'0=*  iF6,2. ' t ' 1 

53  1MK1.EU.31  WHITE  (6.33) 

33  Format i * pl  no  tplimvi  tah(Mv)  lhivulisi  lumvi  caki 

1 1C  TCS(MV)  STS  Sl/SIS  •) 

IMK1.NE.3I  WHITt(6.34) 

34  FORMAT  1 * PL  NO  TPL(My)  TAH(Mv)  |LFM(MV)  IFLOIHV)  lNIVOL 

IT)  tc  (MV ) CAST  TC  Tcb(Mv)  STS  Sl/STS  ST. STO  l,N 

MWlD/UM 

00  13  NF1.24 

1F(STH(NI .GT.0,)STHA0=STN(N)/STK(N) 

1MSTH(n).EU.0.ISTHAO=0. 

IMHEA(N)  .LT.HEXOIll  )G0  TO  70 

IF ( (NST .tU.O ) .OH. (Kl.EU.3) ) GO  10  &U 

UO  72  152. 1TAB 

IF ( RfcX (N) ,LT . HEXO 1 1 ) I GO  TO  71 

72  CONTINUE 
SroO=STOUTAB> 

GO  TO  73 

71  HATOOS<HEXU(l-HEXo(l-l>)/<HEXO(I)-HFyO<i-l>  > 
S[uO=STO(l-l)F($TQ(I)-STO<i-l) IwKAIOn 
GO  TO  73 
70  SI0U=ST0I1) 

73  STSTO=STN<N)/STOO 

IF  (iltE(N).GT,0.)FLNB=ALOG(UEC<  .)*l.)/UtElN) 

IF  (Bl.t(N)  .LL.O,  IFLNUFO.O 
IK (N.OT .12)00  TO  S9 

IF  (Kl.E0.3)*HITC<w.»C.i„  i PL  IN)  i TAR  (N I »EH (N)  .£S(N)  .N.TCS(II)  iSTH(ll) . 
lMRAO 

32  FORMAT  (UP  2F10,3.Fl0.4.Fl0.3>I10.F10.3.Fl0.bt3F-lU.i) 

IF  ( K 1 « NF  ,5) WHITE (6. 1 BIN. TPL 111)  .TAR INI » TCFM(N) . LFFI IN)  .LU(N)  .LSINI  .N 
I.TCSlN' .strini.sthao.ststo.flnu 
IS  FORMA’  ( ltO«4FlO,3.FU>.9»FU>.3« 1 10  »F  IP  .3.F1U . t>.  3F1U  .3 ) 

GO  TO  13 

69  IF  (X  , .EO.  3)  WHITE  (6. 39)  N.  TPL  (N) » TAR  (N)  *fcH(  N) . tS  (Ml » S IH  (I.) . STHaO 
39  F'OP  lA T(ll0.2FlU»3.F10»4.F  1U.3.20X.F1  (’ *5* 3F 10, 3) 

IK  iKl  «Nt . 3.' WH I TE 1 6. 40 1 N.  TPL  IN) . TAR  IN)  * TCFM1N)  .CFKIIl)  .LNIN) . tSUll  .S 

1 Tl  |N)  .SThAO.STSTO.FLNB 

SO  '3HMAT<UOf4F10.3.P'10.4iFiU<3.£QX»Flt>.b.3F10.5) 


GO  TO  13 
SB  CONTtNUt 

IF (N.UT. 12)00  TO  38 

IF  (Kl.EU.31wHlTEl6.32lN.  TPL  IN)  iTAR(N)  * tHIN)  • t$ 1 U)  (N.  TCSllll . STH(U)  . 

ISTHAO 

IF  (Kl.Nt.SlWHITEla.lBlW.TPElN)  (TAHUll.TCFMINI  .LFM(N)  »EM(HI  .ESIN)  ill 
l . ICSlNI .STHIN) .STHAO 
00  TO 

38  IF  ( K 1 .to.  3 1 wH  I Tt  (6. 39 ) U.  TPL  I *1 1 .TAH(Nl.tHlN) .(Sill)  .STHU.I  .SIllAO 

IF  ( Kl  «Nt . 3 ) WHITE  16. 40 ) i). TPUN) . TAM ( N| . T CFM( Fll  .LFMIUI  t LUND  .1  Sill)  .6 
1 T H INI >STHAU 
13  LUnTINUE 
RLTUHN 
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